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SUMMARY 

The paper sets out to give a comprehensive account of the effect of woven wire 
gauze on small disturbances to a uniform stream. The relevant aerodynamic 
characteristics of gauze are the pressure drop and deviation in direction which it 
imposes on a uniform stream incident at an arbitrary angle to the plane of the 
gauze. Measurements by Dryden and Schubauer suggest that these two charac- 
teristics may be uniquely related for all gauzes and Reynolds numbers, and enable 
the deviation to be determined when the pressure drop is known. 

The effect of placing a gauze at right angles to a uniform stream, on which is 
superimposed a small steady disturbance to the longitudinal velocity, is determined 
in terms of the gauze characteristics. 
the problem and the effect of gauze 


The limitation to a small disturbance linearizes 
is to reduce the intensity of the disturbance 
While leaving its (arbitrary) shape unchanged. Previous attempts to solve this 
problem have been made; it is shown that they all make special assumptions about 
the side-force exerted by the gauze and that Collar’s analysis is incorrect on another 
acount. The theoretical reduction in the disturbance is consistent with the few 
published measurements. There is one gauze—with a pressure-drop coefficient of 
about 2-8—which entirely removes all small steady longitudinal disturbances. 

The effect of gauze on a small arbitrary unsteady disturbance—i.e. on turbulence— 
is also determined, on the ass imption that the gauze wires produce no wake turbu- 


lence. The intensities of fluctuations parallel and normal to the plane of the gauze 


are found to be reduced by different amounts, and if the turbulence is isotropic far 
upstream from the gauze it is axially symmetrical downstream from the gauze. The 
Teduction in turbulent energy 


depends on the three-dimensional energy spectrum of 
the oncoming turbulence and is thus not unique, even when the oncoming turbulence 


[Quart. Journ. Mech. and Applied Math., Vol. II, Pt. 1 (1949)] 
5092.5 B 




















2 G. I. TAYLOR AND G. K. BATCHELOR 

is isotropic. Some estimates of the reduction in intensity obtained from an approxi- 
mation to the spectrum function in isotropic turbulence are presented. They ar 
found to be close to the reduction factors appropriate to a simplified type of tate: 
lence in which the longitudinal velocity does not vary in the direction of the stream. 
It is possible to deduce the effect of gauze on other characteristics of the turbul nce; 
one immediate result is that the longitudinal integral scale is reduced to zero by the 


gauze which removes steady longitudinal disturbances. 


INTRODUCTION 
Ir has long been known that when a turbulent stream of fluid passes 
through a wire gauze it becomes less turbulent, and that steady distur- 
bances are reduced in intensity. Several attempts have been made 
(Prandtl, 1; Collar, 2; Batchelor, 3; Dryden and Schubauer, 4) to give a 
theoretical description of these effects and recently some measurements 
(Dryden and Schubauer, 4) showing the effect of gauze on wind-tunnel 
turbulence have been published. In all theoretical discussions the gauze 
has been supposed to offer a resistance to the component of flow normal 


to its plane. If a pressure difference p,—p, is required to drive fluid of 


density p through a gauze at velocity U’, the drag coefficient k& is defined as 


and it has been assumed that / measures the only aerodynamic property 
of gauze which is effective in reducing turbulence. 

It was to be expected that different values of the ratio of the disturbance 
velocity in front to that behind the gauze might be obtained when different 
types of disturbance are considered, but two quite different expressions 
for this ratio in terms of k were obtained by Prandtl (1) on the one hand, 
and by Collar (2) and by one of us (3) on the other, for exactly the same 
disturbance. It will be shown later that the difference is due to the fact 
that these authors left out of consideration the force which the gauze 
exerts on the fluid in directions parallel to its surface. Prandtl and Collar, 
in fact, make different implicit assumptions about this force. It will be 
found that when allowance is made for this side-force Prandtl’s and 
Collar’s expressions are special cases of a more general formula. 

The work falls into three parts. In Part I a method is described by 
which the side-force as well as the drag of a gauze can be measured when 
a stream passes obliquely through it. Measurements made at the National 
Physical Laboratory by Simmons and Cowdrey (5) using this method are 
described and discussed. The results of more recent and more extensive 
measurements of the same kind by Dryden and Schubauer at the National 
Bureau of Standards have been communicated to the present authors and 


are reproduced in an appendix to this paper. Part II contains a solution 
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THE EFFECT OF WIRE GAUZE ON SMALL DISTURBANCES 3 


of Prandtl’s and Collar’s problem, viz. the effect of gauze on a small 
steady longitudinal disturbance, and a comparison with some experimental 
results. There is one particular value of k for which longitudinal distur- 
bances are entirely removed by the gauze. In Part III we discuss the 
effect of a gauze (whose mesh is so small that it creates no turbulence) on 
general unsteady disturbance, i.e. on stream turbulence, and expressions 
for the reduction of the mean-squares of the velocity components are 
derived. These reduction factors depend on the energy spectrum of the 
turbulence approaching the gauze. In the case of isotropic turbulence, 
previous measurements enable us to postulate a spectrum function which 
is likely to be approximately valid under certain conditions, and the effect 
f gauze can be calculated. A comparison with the experimental results 
of Dryden and Schubauer (4) is given, though owing to the nature of these 


experiments it is not possible to draw any definite conclusions. 


Part |] 
Weasurement of side-force and drag on a gauze sheet set 
obliquely to a uniform wind 
In the first instance it will be assumed that the flow is confined to one 


plane perpendicular to the gauze. This assumption implies that the gauze 


has aercdynamic properties which are symmetrical about the plane of flow. 
The flow through the gauze gives rise to a pressure difference p,—p, across 
its two sides. If l’ denotes the velocity of the stream approaching the 
gauze and @ the angle its direction makes with the normal to the plane 
of the gauze, the results of experiments in which p,—p,, U, and @ are 
measured can be ¢ xpresse d in the form 

Pi— Po = kg.4pU?. (1.1) 
The coefficient ky is identical with the k previously used when 6 = 0. 
Measurements show that is a function of 6 and U for any particular 
Yauze, 

The force parallel to the gauze gives rise to a change in the component 
of velocity parallel to its plane as the fluid passes through it. To measure 
this force it is sufficient to measure the change in direction of the stream 
on passing through the gauze. If ¢ is the angle which the stream makes 
with the normal to the gauze on its leeward side (see Fig. 1), the air 
emerges from the gauze with velocity U cos@sec¢ and the side-force per 
unit area of the gauze is 


ol nécosé Sil cos-0 sec d) = pl 2cos 6 sec d sin(@ —¢). 
If can be measured, the results of experiment can therefore be expressed 


In the torm . . 
in the fort i 2 cos see d sin(@—¢), (1.2) 
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G. I. TAYLOR AND G. K. BATCHELOR 
where F, is a force coefficient parallel to the plane of the gauze. In cases 
where steady or unsteady disturbances are producing transverse com- 
ponents of velocity which are small compared with the mean velocity, 
6 and ¢ are small, and (1.2) can be written 

F,/0 241 (d Q)!. (1.3) 


This form is useful because the experimental results show that 4/@ tends 
to a finite limit as @ tends to 0. The value of 4/@ for small values of @ 





c 


(= a, say) and the resistance coefficient k are the two parameters which 
will later be taken as specifying the effect of gauze on small disturbances. 

In general, if a gauze is set up at angle @ to a stream in a wind tunnel 
with rigid walls, the direction of flow will vary along the length of the 
gauze on both sides of it. For convenience in measurement it is best to 
try to arrange matters so that the velocity and direction of the stream, 
and consequently the pressure, are uniform on both sides of the gauze. 
One way in which this may be done is by turning the walls of the channel 
(shown as AC, BD in Fig. 1) about hinges at A and B till they make the 
appropriate angle ¢ with the normal to the gauze. When the walls AC, 
BD are set at this angle, it is possible for the velocity, direction, and 
pressure of the streams on the two sides to be uniform. In practice the 
experiment consists in determining the angle ¢ at which the walls must 
be set in order to make the variation along the length of the gauze as 
small as possible. 

This method was adopted by Simmons and Cowdrey (5) for measuring 
the values of kg and F, over a series of values of 6 and U for three gauzes, 
referred to by them as G,, G,, and G,, the geometrical characteristics of 
which are given in Table 1. The directions of flow up- and down-stream 


from the gauze lay in a vertical plane in each case. The measurements 
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TABLE | 
Geometrical characteristics of gauzes used by Simmons and 


Cowdrey, and Dryden and Schubauer 





Wire diameter Number of wires pei incl 
In 

<2 40 

52 (horiz.) 

- ( ’ 0067 ‘ 

- ‘ 60 (vert.) 
\ ¢ O75 40 
05 60 
ss 54 
S 005 4° 
7 24 









Gauze G, 





30 9° 40 50 60 
| 2. Measurement fd and @ made by Simmons and Cowdrey. 
f 4 for a given value of 6 were found to be insensitive to variations of U 


ver the range 8 to 34 ft. per sec. Their results for the three gauzes are 


shown with their permission in Fig. 2. It will be seen that in all three 


cases 4/8 is nearly constant over a large range of values of 6, viz. 


tor gauze G 66 0-48 (0 <= 6 < 20°), 
for gauze G 6/6 0-66 (0 <6 60°), 
tor gauze G, b/60 0-70 (0 <6 =< 60°). 


Dryden and Schubauer have recently adopted the simpler plan of 
eliminating the walls of the duct downstream from the gauze so that the 
stream discharges into free air. The jet emerges at an angle ¢ to the 
normal to the gauze and this angle can readily be observed. Dryden and 


Schubauer have been kind enough to send to us an account of their 
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experiments which is reproduced as an appendix to this paper. The 
significant feature of this extensive series of measurements is that //@ 
appears to be uniquely related to kg for the four gauzes and the various 
wind speeds used.+ An empirical relation which seems to fit the measure- 
ments sufficiently well over the range kg > 0-7 is 

F,/@ = 2—2-2(1-+-kp)? (1.4) 
which is shown in Fig. 5 in the Appendix. In our discussions of the effect 
of gauze on small disturbances, we shall be concerned only with small 
values of 0, for which (1.4) leads to (see (1.3) ) 

x = 1-1(1-+k)-. (1.5) 

It would be very useful if, as the evidence seems to suggest, this relation 
were valid for all wire gauzes and Reynolds numbers in common use, for 
the aerodynamic ‘smoothing’ properties of a gauze would then be com- 
pletely specified by the easily measured quantity k. 

It is probable that the above experiments were confined to the case 
of gauzes possessing geometrical symmetry about the plane of flow, i.e. 
one of the lines of wires was in the plane of flow while the other was 
perpendicular. We have no experimental evidence that if the gauze were 
rotated in its plane the downstream flow would remain in the same 
direction or even in the plane containing the upstream flow and the 
perpendicular to the gauze. However, for simplicity in the succeeding 
analysis, we assume that the gauze has aerodynamical properties which 
are rotationally symmetrical about the perpendicular to its plane. It will 
also be assumed that the same value of « applies at each point of the gauze 
in cases when the flow velocity is uniform neither in space nor time, and 
that the distance between the wires of the gauze is small compared with 
the scale of non-uniformities in the stream passing through it. 


Part II 
Effect of gauze on a small steady disturbance 

1. Solution of the problem 

The disturbance investigated by both Prandtl and Collar is a small 
steady spatial variation in the magnitude, but not the direction, of the 
velocity of a stream, and so affects all points along a streamline. Accord- 
ing to Prandtl (1) the ratio of the velocity increment (u,) to the mean 
velocity at any point downstream to that upstream on the same stream- 


line (1) is Us (2.1) 
—_ = =» “s 
Uy 1k 
+ A plot of J/sin@ against & shows that these quantities can equally well be represented 
by a single relation. 
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esult (2) for the same type of disturbance is 
Us 2—k 
: (2.2) 
) . 
u, 2+ 


of us has obtained the result (2.2) by a different method 
Collar considered the way in which the longitudinal 


d by the gauze without considering the lateral velocity 
icquired by the stream both upstream and downstream in 


d of the gauze in order to satisfy the condition of con- 


iwitudinal velocity varies. Batchelor took the induced 


nto account but made the assumption, now known to be 


that these lateral velocities are continuous across the gauze. 


ution of the steady disturbance problem can be given if 


ed that the disturbance velocity wu, is small compared with the 
el U, and tl the gauze behaves as in the experiments 
lescribed. so th it when 0 and d are small 
db 0, (2.3) 
is a constant fo1 particular gauze and for arbitrary small 
of the streamlines impinging on the gauze. 


wimuth 1] ungies 


der first the case when the initial disturbance velocity, i.e. the 


nt + +} 


mean speed far upstream, depends on the lateral co- 


und has nusoidal variation, so that 
u Uy COS py. (2. 1) 
of disturbance velocity near the gauze must also be 


wave-length 27/p. If the effect of viscosity be neglected 


xcept in so far as it is instrumental in giving rise to resistance at the 
gauze, the ition of moti can be written 
( ¢ — * 
l v V “Wy Q. (2.5) 
i cy 
here ub, stream function for the upstream disturbance velocity, 
represents distance in the downstream direction from the gauze, and 
€ |, 
¥ 29 
ci COX 
Neglecting ontaining squares of small quantities, (2.5) assumes the 
m 
v=) 0). (2.6) 
hat ; function of y only. Far upstream wu is given by (2.4), so 
it the fli ( the vauze civen by jute 


Uy pP sin PY. 

















G. I. TAYLOR AND G. K. BATCHELOR 
The solution of this equation which does not become infinite for large 
negative values of x is 
ub — isin yy +- Ae?” sin pr 2.8) 
Py p PI é sin py. (2.8) 
(The other solution A’e”* cos py is not needed in view of the symmetry of 
the flow about y = 0. Nor are the solutions representing higher harmonies 
needed since all the equations are linear and their coefficients would be 
found to be zero.) Similar considerations show that the stream function i, 
for the disturbance downstream of the gauze must be 
uso san M2 Sin py+- Be-* sin py, (2.9) 
Pp 
where uw, cos py is the disturbance velocity far downstream and A and B 
are constants to be determined. 
Now there are three conditions available for the determination of A, B. 
and w,/u,. The condition of continuity requires that w should be continuous 
across the gauze, with value uw, cos py, say, so that 


Ue Uy pa Us pb. (2.10 


The equation for the change in v on passing through the gauze is, from (2.3), 


/ / 


CW cus 
= X : at a 0 
COX CX 
i.e. B-+-xA 0, (2.11) 
YU, Us , 
and (2.10) becomes Us =. (2.12) 
~G 


Bernoulli's equation for an arbitrary streamline equates the difference in 
the total pressures far upstream and far downstream to the drop in static 
pressure at the gauze. Thus, to the first order of small quantities, 

P+ pUu, cos py—p,— pUu, cos py = ktp(U?+-2U ug cos py). (2.18 
But since the streamlines are initially and finally straight, the asymptotic 
static pressures p, and p, are independent of y, and p,—p, must have the 


stun bintit nce 
value kipU?. Hence Us —te = kus. (2.14) 


From (2.10), (2.11), and (2.14) it will be found that 





{ aa k | B Uy vk (2.15) 
. p\l x i} p \l+a- ky? 

and’? Us l X vA (2 16 
Uy | vtk i 


In the course of Collar’s work (2) he had some correspondence with one of us (G. I. T.), 
who pointed out that the difference between his formula and that of Prandtl was associated 
with the transverse force on the gauze. In a letter (9 Feb. 1939) Collar then proposed the 


formula (2.16) as a possible general expression containing both formulae. 
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Though it has been assumed that the disturbance is sinusoidal in y, the 
wave-length p does not appear in the result; and since in the first-order 
alysis small disturbances are superposable, it appears that any distur- 
ance of the Prandtl—Collar type will be exactly reproduced in shape far 


wnstream from the gauze but with the velocity altered in the ratio 


One particular case of a small steady disturbance to a uniform stream 


hich is not covered by the preceding analysis consists of a disturbance 

which the motion is purely lateral. A longitudinal vortex in which the 
iteral velocity is everywhere small compared with U is of this type. In 
his case there is no superposed irrotational flow near the gauze and the 
effect of the gauze is simply to reduce all lateral velocities to a fraction a 
f their value upstream of the gauze. The effect of the gauze is here 
ndependent of the resistance coefficient k, except in so far as « and k are 
related, because to the first order the magnitude of the velocity is U 

erywhere and the pressure-drop across the gauze is the same for all 
streamlines. It follows from the linearity of the problem that for a 
disturbance consisting of both longitudinal and lateral motion (the former 
f which must be steady, i.e. independent of x) the gauze alters all longi- 
tudinal velocities in the ratio (1+-~—ak)/(1+a«+4) and all lateral velocities 

the ratio a. The shape of the flow pattern far downstream will not be the 


same as that far upstream in this general case. 


2. Comparison with previous analysis 


It will be noticed that (2.16) reduces to Prandtl’s result that 
Us ] 
Uy l+k’ 
when a 0, and to the result obtained by Collar and by Batchelor, viz. 
tle » hk 
Uy 2+ k’ 
when |. In his proof Prandtl assumes that the flow downstream of 
the gauze is parallel to the mean flow, and that the pressure there is 
consequently uniform. This condition is fulfilled only when a 0, which 
is likely to be true for real gauzes only when k x. Batchelor’s proof, 


which is not unlike that used in the preceding section, involved explicitly 
the assumption that the lateral velocities are continuous at the gauze, 
. that a 1, which is also untrue, in general, for real gauzes with a 


hinite resistance coefficient 
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On the other hand, Collar’s analysis, which leads to the correct result 
for the particular case « = 1, appears to contain an error which we shall 


examine for its intrinsic interest.t He discusses the flow in a parallel- 


walled channel in which a narrow filament has velocity uw, in excess of 


that of the main stream. The filament passes through the gauze with 
velocity uz and finally settles down to velocity uw, far downstream. The 
rate of change of momentum of the filament is then equated to the excess 
drag at the gauze due to the fact that the filament is not moving with the 
same speed as the mean flow. In this way Collar derives an equation 


which, when associated with the equation 
1— Ue = hug (2.17) 


representing the effect of applving Bernoulli’s equation (see (2.14) ), vields 
the result Duy, (2.18 
which is identical with (2.12) in the case a L. 

This method of Collar’s depends on the assumption that the change in 
momentum of the stream outside the filament is zero. This assumption 
is not true.t The change in momentum of the stream outside the filament 
is, so far as the first-order term is concerned, exactly equal and opposite 
to that of the filament. If U is the mean velocity, and the stream has 
constant total width b, and if uw, is the variable disturbance-velocity far 
upstream, while w, is that far downstream, the rate of change in momentum 
of the stream is , 

p | [(U u,)? (U Uy)*| dy. 
By definition of U, 
b 
}u,dy=0 and | udy = 0, 
0 0 
so that the rate of increase in momentum is a small quantity of the 
second order when w,/U is small. The expression used by Collar is the 
part of the momentum integral lying within the filament. This is of 
the first order and therefore, when w,/U is small, is large compared with 
the true expression which is of the second order. The same considerations 
show that the difference between the resistance of the gauze when opposing 
a uniform stream, and that which occurs when the stream has a small 
Professor Collar informs us that he originally obtained his result by another approach, 
in which he used a single filament and adopted the assumptions of the well-known 
Rankine—Froude momentum theory of propellers: these assumptions are consistent wit! 
the results obtained by us for the particular case « 1. Subsequently he obtained experi: 
mental support for his formula, and realizing the objections to his original theoretical 


approach, published a different analysis which was wrongly used to obtain the same result. 


H. L. Dryden has also drawn attention to the error in private correspondence. 
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variable velocity, is also of the second order of small quantities. Thus the 
nomentum equation applied to a parallel-walled channel merely gives a 
ull result so far as the first-order terms are concerned. 

Though Collar appears to apply his momentum condition to the whole 
stream, by neglecting the change in momentum outside a stream tube he 
really only applies the condition to the fluid inside the tube, thus vir- 
tually making the assumption that the pressure on all parts of the wall 
f the tube upstream of the gauze is p, and on all parts downstream py, 
p, and p, being the mean pressures on the two sides of the gauze. The fact 
s that the pressure varies along the surface of the stream-tube. This 
variable pressure makes a contribution to the momentum equation for the 
stream-tube which Collar’s proof takes to be zero. The contribution would 
n fact be zero if and only if the distributions of difference in pressure from 
the mean on the two sides of the gauze were mirror images of one another 
vith a change in sign. When only first-order terms are considered, such 

listribution of pressure implies that the gauze produces a disturbance 


vhich is antisymmetrical on opposite sides, so that a velocity (U-+u,+4u, v) 


the point v,y) on the upstream side corresponds to a velocity 
U+u,—u,v) at the point r,y) on the downstream side; this is just the 
vpe of flow shown by (2.8 2.9), and (2.11) to exist in the case « : 
\t the gauze the w-component of velocity is the same on both sides and 
equal to | so that 

ly l l Us—(U)-= l Us 

This implies that Us +Us = 2Us, (2.19) 
which is identical with (2.18). Thus Collar is virtually assuming (2.19) to 
be true and his proof obtains the correct result (for « = 1) by virtue of 


the fact that (2.19) is in fact valid for a gauze which exerts no side force. 


3. The flow in the neighbourhood of the gauze 


We have seen that the assumption of antisymmetrical disturbance- 


velocities on opposite sides of the gauze includes the assumption that v 
s unchanged on passing through the gauze, or, in other words, that « = 1. 
this raises the query, How is the antisymmetry modified when «a + 1? 


lhe answer is to be found in the analysis which we have already given, 
is shown by the following considerations. 


From (2.6) it is seen that to the first order of small quantities the effect 


{the gauze on the upstream flow can be represented by superposing on 
the flow at infinity an irrotational motion specified by the potential d(x, y). 


Similarly the downstream flow can be represented by a combination of the 


ow tar downstream and an irrotational motion d'(x, y) in the neighbour- 
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hood of the gauze. If the transverse velocity is continuous across the 
gauze at x = 0, the appropriate relation between the two potential flows is 
d(x, ¥) d(—a, y). 


Likewise in the case when « ~ | the required change in transverse velocity 
across the gauze will be satisfied by 











db’ (2, y) vp(— 2, y). (2.20) 
The condition of continuity of the u-component of velocity at the gauze 
gives 

d(x, y) Cd’ (x, y) d(x, y) 
Us U4 : Uo - Us x} — : 
Cx r=0 CA r=0 Cx r=0 

in the general case. Hence 

x(U3— Uy) Us—Us, 2.21) 


which is identical with equation (2.12) already established, and together 
with the Bernoulli equation (2.14) yields 


Us 1lt+-a—ak 


—_ (2.16) 
1 Iltatk 


u 


as before. It will be noticed that 4(x, y) is determined by the distribution 


of u far upstream, for from (2.21) and (2.16) 


od(a, y) Ita 
: Us Uy R (2.22) 
Cx r=0 i Tar k 


These formulae are still applicable if wu, is a function of y and z, provided 
that the gauze has the property that the change in direction of the stream 
on passing through it is a function only of the angle between its plane and 
the direction of the incident stream. 


4. Comparison with experiment 
Measurements of the effect of gauze on a small steady disturbance are 
difficult to make for several reasons, not the least of which is the fact that 


wire gauzes tend to distort, producing non-uniformity in k. Two sets of 


measurements have been published, and although there is considerable 
scatter of the points, they provide some support for the validity of the 
result (2.16). Collar (2) measured the reduction in amplitude of a single 
disturbance for a number of gauzes, whose geometrical characteristics are 
given in Table 2, and his results are shown in Fig. 3. MacPhail (6) tested 
the effect of several gauzes (specified in Table 2) on a complex disturbance 
in a small water tunnel, and obtained the reduction in magnitude of the 
whole disturbance. His measuring points were at equal distances upstream 
and downstream from the gauze, and since with no gauze present the 
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TABLE 2 
Geometrical characteristics of gauzes used by 


Collar and by MacPhail 


miele Vumber of eres per inch 





2) Collar 
\x MacPhail 























Le °) 





3. Effect of gauze on a steady longitudinal disturbance. 


disturbance diminished by a factor 0-35 owing to viscous decay, we have 
assumed that his measurements of u,/u, must be divided by v0-35 to 
obtain the effect of the gauze alone, the assumption being that a negligible 


amount of viscous decay occurs downstream from the gauze. Both Collar’s 


ind MacPhail’s results show fair agreement in Fig. 3 with the theoretical 
ratio (2.16) if we make use of the empirical relation established by Dryden 
ind Schubauer, viz. a l-1(1-+-k)-? (k > 0-7). Also shown in the figure 
are the curves describing the reduction in the special cases a = 0 and 


‘= 1, neither of which fits the experimental points as well as the more 
general result 

The important practical point that small longitudinal disturbances are 
entirely removed by a gauze having a particular value of k appears to be 


supported by the experiments. This optimum value of k is given by 
k ] l/a, 
le. k = 2-76, if the gauze satisfies Dryden and Schubauer’s empirical 


relation between « and &. Lateral disturbances are not removed by such 


a gauze, of course, so that the choice of the most effective gauze for 
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the reduction of a general steady disturbance must involve a considera- 
tion of the relative amounts of energy in the lateral and longitudinal 
disturbances. 


Part II] 
Effect of gauze on turbulence 
1. Outline of method 

The method used in Part II of this paper to investigate the effect of a 
gauze or mesh on small steady disturbances in a uniform stream can be 
adapted to suit the case of unsteady disturbances. It will be necessary 
to assume that the turbulence far upstream of the gauze is homogeneous 
and has a root mean square velocity fluctuation which is small compared 
with the velocity of the uniform stream. As in Part II, the gauze will be 
assumed to be placed at right angles to the stream direction, and to be 
characterized by a resistance coefficient / and a side-force parameter « for 
incident streamlines of any (small) azimuthal angle. The treatment is 
essentially kinematical; that is to say, we take no account of change of the 
turbulence due to viscous or dynamical forces (other than those exerted 
by the gauze). The procedure is thus more accurate when the relative 
intensity of the turbulence is small, since the effect of the gauze depends 
only on the turbulence pattern and not on the speed of the uniform stream 
which carries it along. An attempt to solve the problem on roughly similar 
lines has previously been made by one of us (3). However, this attempt 
was in error inasmuch as it ignored the effect of side-force exerted by the 
gauze, and, moreover, the analysis was not wholly free from approxima- 
tion. The analysis given below overcomes the first and improves the 
second of these defects. 

With the assumption of small relative intensity of the turbulence, the 
effect of the gauze is linear, as in the case of a steady disturbance. Super- 
posed velocity patterns will therefore experience independent modifications 
due to the gauze; we take the hint and express the turbulent flow pattern 
as a Fourier integral of (spatially) periodic distributions of velocity. 
Similarly the velocities induced by the presence of the gauze can be 
expressed as a Fourier integral of periodic variations of velocity. We use 
the conditions that the lateral components of the velocity at each point 
of the gauze are diminished discontinuously to a fraction « of their value, 
that the local longitudinal velocity is continuous at the gauze, and that 
the local pressure-drop across the gauze is determined by the local longi- 
tudinal velocity and by the resistance coefficient k. These conditions are 
sufficient to relate the amplitudes of similar periodic disturbances in the 
turbulence far upstream and far downstream from the gauze. The integra- 
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ine the total intensity of turbulence far downstream from 
ires a knowledge of the three-dimensional energy spectrum 


turbulence, although an approximation to this function 


1 


vhen the upst ream turbulence is ise yt ropic. Since the energy 


tropic turbulence can vary (possibly within fairly narrow 


ie effect of a gauze also varies with the conditions. 


2. Representation of the velocity field 


id z-axes in the plane of the gauze, with the positive x-axis 


tream. Far upstream and far downstream the velocity 


the disturbance are w,, ¢,, w, and Uy, Vg, Ws, respectively. 
yauze the only change in the velocity field with time is a 
wnstream so that these six components are functions of 

and N if the disturbance velocity is sufficiently 


| with U’, vorticity is carried along by the uniform stream 


nge, except through the gauze. It follows that the velocity 


din th ghbourhood of the gauze has components 
j 
COM : 
wy- 1 for rE< 
{ C2 
l comp 
j 
COs, . 
te. = for r 0, 
- Cz 
Here 4, ire potential functions (of x, y, z, and t) which must tend 
mstat ues as x approaches —o and +-o respectively. 
We represe the distribution of disturbance velocity far upstream from 
gauze the form of a triple Fourier integral over all values of the 
wve-numbe mponent m, n: thus 
l, i | | lm, njetrmu-n) dldmdn, (3.1) 
re A, und A, m,—n) must be complex conjugates in 
er to make u, real. Likewise the components v, and w, are determined 
functic n) and C,(/,m,n), and the condition of incompressi- 
Ivy requ i, mB, nu; 0, (3.2) 
the sa the components us. Vo. w. of the disturbance velocity 
downst1 may be written as Fourier integrals like (3.1) and the 
resp 0 instorm I tions As, B,, C, are related by 
mB,-+-nC, 0. (3.3) 


m (3.1) is not rigorously possible when the turbulence 


| spatially homogeneous in the statistical sense, 
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for the functions A,, etc., are then infinite. However, we shall be interested 
in ratios (like A,/A,) only and the divergence of the functions may be 
safely ignored. The analysis remains unchanged if we think of the velocity 
field as large, but not infinite, and the difficulty does not then arise. 

It will be seen in the next sub-section that the equations connecting 
the disturbance velocities at different positions relative to the gauze are 
linear and homogeneous, as was found for the steady disturbance. There 
will therefore be no modulation of the Fourier components, and we can 
discuss one particular wave-number in isolation. The effect of gauze on 
a given field of turbulence is thus known when the effect on the triply 
periodic disturbance 

(3,01, W,) = (Ay, By, Cet +mvtn2) didmdn (3.4) 
is known. Far downstream from the gauze, this disturbance will be 
converted to 


(Uy, Va, Wy) = (Ag, By, C,)e*+mu+> dldmdn, (3.5) 


and the problem is to determine the functions A,, B,, and C, in terms of 
the known functions A,, B,, and C,. The relevant parts of the potential 
functions ¢, and ¢, have the same periodic character in ¢, y, and z, so that 
corresponding to (3.4) and (3.5) we put 


$d, = Pl, m, n)er2+*Ctmy+n) dldmdn, (3.6) 


dy P,(1,m, ne er (AC my+n) dldmdn, (3.7) 
where oc? m?+-n®; the dependence of ¢, and ¢, on x is obtained from 
Laplace’s equation and the boundary conditions at x x. and +x 
respectively. 

The characteristics of the turbulence which are of most immediate 
importance are the spatial means of the squares of the component fluctua- 


tions, e.g. uj. Now the average value of u? over a large cube of side 2D is 


D 
“rer . 
| uz dddydz 
8 D8 
D 


a [ [ [ [ [ A,(/,m,n)A,(l',m’, n’) > 


sin(/ ‘ SI ‘ Si 1 99’ - , 
in l')D sin(m m \D sin(n nD tdmdndl'dm'dn 


+0 m--m nn 
| Ve ( “ee ; 
D3 | | aut S), | | A,(L l, M@—m, N—n)> 


: sin LD sin MD sin ND ILAMAN | diduda. 
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ested where L +l’, M m-+m’, N n+n'. Thus when D is sufficiently 
Vv br irge L 
OCH us ae | 1,(/,m,n)A,(—l, —m, —n) dldmdn 
D8 

cti 

- A,(l,m,n) * dldmdn. (3.8) 
There D8 
' Hence if we define an energy spectral density F,(/,m,n) such that 
‘ ; 

rij F(l.m.n) dldmdn is the contribution to uf from wave-numbers in the 
Li} 

of 1 m.n) to (1 dl.m tm. 7 dn a then 

) 3 

0.4 F\(l,m,n) D3 A, = (3.9) 
1 

Likewise G,, H, and F,, G,, H, are proportional to B,|?, |C,|? and |A,/*, 
35 B, 2. C,\2 respectively. The mean energy of the turbulence per unit mass 
_ f the fluid far upstream from the gauze is 
Pnit : : ; : _" 
4] L(uf+eptuyz b ||| (4,4+-46,+4,) dldmdn, 
with a similar expression for the energy of the turbulence far downstream. 
3.7 3. Conditions at the gauze 
frol The conditions to be satisfied at a 0 enable us to relate the velocity 


fields on either side of the gauze. The first condition is that the lateral 
velocities on the downstream side of the gauze are a fraction « of their 
diate value on the upstream side. Thus 


ictu Cds Cd, 
97) 1s 2 : 1 
7” oY CY }o 


| C ds od, 
|e . *) ri w,— — * 
2 2 1 pia 
. Cx }9 
> 


nd for the single Fourier component described by (3.4), (3.5), (3.6), and 


Bb, imP, ( B, . imP,), (3.10) 
C,+inP, (C, +inP,). (3.11) 


The second condition is that the velocity normal to the gauze is the 


COs od, | 
{ 2 | = u.—-— 9 
2 : 1 : 
\ CA Cua 0 
ind for the typical Fourier component 


A,—oP, = A,+oFP,. (3.12) 


same on either side, i.e. 
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Finally, there is the condition that the instantaneous local pressure-drop 
across the gauze depends on the local velocity. Consistent with our basic 
assumptions, the equations of lateral motion can be written 


(: | af s| 
= (v. w) —-—{|—,—]P, 
ot Cx ploy cz 


so that if A represents the change in any quantity at the gauze, 


a(; +l = Te ur) _ f ’ . )ap = [ K(; m Jw, ws (3.13) 
ct Cx, p\cy cz Cy Cz 


There is no change in 7, w , V2, or w, following the mean motion, so that 
only the induced potential flow contributes to the left side. For the 
typical Fourier component, (3.13) leads to two identical equations of the 
form lF,—F,)—io(F,+F,) = ik(A,+oP)). (3.14) 

In addition to the four independent equations (3.10), (3.11), (3.12), and 
(3.14), there are the continuity equations (3.2) and (3.3). Thus, if A,, B,, 
and CO, are regarded as known, there are available five equations to deter- 
mine the functions A,, B,, C,, P,, P, specifying the disturbance far down- 
stream and in the neighbourhood of the gauze. 


4. Solution of the equations 


From (3.10) and (3.11), and the continuity equations (3.2) and (3.3), 
we find io*P,—1A, = a(io®P,—IA)). (3.15) 
P,, P,, and A, can now be found from this equation and (3.12) and (3.14). 


A, ofB—B—ik 


P. ‘ 3.16) 
. o 2B+i1(1+a+h) ; 

i a A, (B : vy B- B+-iak) ‘ (3.17 
~ o (B—t)[2B+7(1+a+h)] 

A, rete x—1)] (3.18 


(B—2)[28+-71(1+a+hk)] 
where f has been written for /(m?-+-n?)-*. Hence 


1.?8?+-(1+-a—ak)? > 


F,(l, m,n) A,|* _ (3.19) 
A,/? 462+ (1+a+k) j 


F,(l,m,n) 

I(B), say. 
Thus, as a result of the passage through the gauze, the contribution to the 
mean square of the turbulent fluctuation parallel to the main flow from 
Fourier components with wave-numbers in a small range about (/,m,”) 18 
reduced by the factor J(8). This factor increases monotonically from 
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{1 me vk , ) . ‘ . 
: at B 0 to a? at f x (for wire gauzes), and is always less 
| Y . j 
than unity. The case 8 = 0 (i.e. ] = 0) corresponds to a velocity distur- 


bance which has no variation in the #-direction, and the result (2.18) 
obtained in Part IT is recovered from (3.19) in this case. 

The reduction ratio for the total intensity of the longitudinal fluctua- 
tion, arising from all Fourier components, is thus given by 


( [ ( F;, ,m,n) dldmdn ( ( I 3) F, dldmdn 


es = ; 
i = E (3.20) 
|| | F\(, m,n) dldmdn || | F, dldmdn 
ind is entirely specified by the spectrum of the upstream turbulence 
und the gauze characteristics k and a. Numerical results will, of course, 
need some assumption about the upstream turbulence, but before pro- 
eeding further with (3.20) we shall consider the effect of the gauze on the 


lateral components of the turbulence. 


From equations (3.10) and (3.11) 


nb, mC, x( 72 B, mC’). 
Multiplying this equation by its complex conjugate and making use of the 
mtinuity equations (3.2) and (3.3), we find 
2 B,|? C,|?) [2|A,|? >| 0? | B,\*- C,|*) -[?|A,|?], (3.21) 


so that the energy of the lateral motion of the typical Fourier component 
is reduced in the ratio 
B,|?+- |C,|? A,|* - 
es = 1 —_[I(p)—a*]. 
B, ( 1 > Bb, ——_ C1 = 
In terms of the spectral densities this becomes 
G,+AH, F i 
t. = —+_|[I(B)—a?]. (3.22) 
G,+H, G,+H, 


The reduction ratio for the energy of the lateral motion of all Fourier 


> 


(ee) 


components is thus 


( ( ( ( Gs H,) dldmdn ( } ( B>( I x?) F, dldmdn 
} xu : - 


( ( [ (G,+H,) dldmdn 
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The expression (3.23) is negative for all values of 8 so that v is always less 
than a; v can be determined when p is known, and we must therefore 
return to (3.20). 


5. Special case of isotropic turbulence 

The essential need, if we are to make practical use of (3.20) and (3.24), 
is for information about the dependence of F, on 1, m, and n in the 
turbulence far upstream from the gauze. We need to know how the 
contribution to wu? made by a single Fourier component varies with its 
vector wave-number (/, m, vn). This three-dimensional spectrum function 
has not hitherto appeared in turbulence literature (7) and we know of no 
direct experimental data. Some guess-work will therefore be necessary, 
and in order to put ourselves on familiar ground we shall assume that the 
turbulence far upstream of the gauze is isotropic; this special case also 
has considerable practical significance. Even in this simple case there is 
no unique form for the energy spectrum, which will depend upon the 
origin of the isotropic turbulence in question. Prediction of the exact 
effect of the gauze would require a knowledge of the particular properties 
of the oncoming turbulence, so that (3.20) and (3.24) represent the limit 
of our exact analysis. ‘However, it seems probable from experimental 
evidence that there are definite limits to the variation of the structure 
of isotropic turbulence under different conditions, and we can seek an 
approximate expression for the effect of gauze on all types of isotropic 
turbulence. 

The postulate we have chosen to make in order to obtain a working 
estimate of the effect of the gauze on isotropic turbulence is that 


; - Co? 
F\(l,m,n) = ut — sso (3.25) 


where C and y are constants. This expression satisfies all the formal 
conditions imposed by isotropy and incompressibility of the fluid (7). 
Moreover, the contribution to uw} from periodic fluctuations with wave- 
number components in the direction of the x-axis which lie between / 
and /+-dl is 


x 


dl ( [ F,(/, m,n) dmdn “ae U (3.26) 
a 2((*+-y* 


Oo — ox 


This function of / is the ordinary one-dimensional longitudinal spectrum 
function used in the literature. and the special form (3.26) is the form 
found by Dryden (8) to be in fair agreement with experiment for several 
sets of conditions. The constants C and y can be chosen to give the 
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correct total energy and one length parameter of the turbulence such as 
the integral scale, but they will be found not to occur in the expressions 


for uw and v. 


The postulate (3.25) now enables the expression (3.20) for yu to be 
evaluated. Transforming m, n to spherical polar coordinates 7, 6, % 
such that ) ) ) Loud 
( , / Yr cos m- n- 5s o* r= sin2us, 

3.20 becomes 

rr tan cos* | ¥ vk)? sin? ¢ r= sin*ys es a 
————— es — 9" sin ds dibdi 
aa t COS“ | X kh)? sin" (7- il 





The integrations with respect to 7 can be carried out immediately, giving 


3 [ 4a? cos*4s+(1 k)?sints . . ee 
aan —____sin3s dub. (3.24) 
2} 4c0s%+(1+a+k)?sin®s 
0 
Changing the variable of integration to y (= cosy) and writing 
(] \ vk)? . (1 yk)? 
(1+-a—ak)?—402’ (l+a+kP—4’ 
We have 
) | vA )2 Ly“ oe &- > , 
: —- i] lay 
Zz | k) | “ sa) x ry 
« ij 
| ) dye | \ k)? tr-. 2__ ¢2 n—I1 
nae $(1—7?)(1+ +>. log+ 
| o_4 | vt+k)?—4 “ 2n n+1 


(3.28) 
This can be evaluated directly when the values of k and « for the gauze 
we known. The reduction ratio for the lateral intensity is then given by 
3.24), which becomes, when the upstream turbulence is isotropic, 


1 


y y*+ 4) (1 x—ak)?—(1-+-a-t k)*w}. (3.29) 
The effect of gauze on the total energy of the turbulence is to reduce it 


in the ratio 1(u-+- 2p). 


al 


6. Numerical results 
It has already been pointed out that the measurements of « and k for 
woven round-wire gauze made by Dryden and Schubauer (see Appendix) 


indicate that these quantities are uniquely related, and that, for k not less 
than about 0-7, 


‘ 1-1(1+-k)-3. (3.30) 
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On the assumption that this relation holds accurately, the turbulence-reduc- 
ing properties of a gauze are specified by k alone. The reductions in longitu- 
dinal and lateral intensities of the type of isotropic turbulence specified by 
the three-dimensional energy spectrum function (3.25) have been worked 
out for several values of k with the aid of (3.28), (3.29), and (3.30) and are 
shown in Table 3 and Fig. 4. The variations of » and v with k have many 
features analogous to those for the steady disturbances discussed in Part IT. 

It is worth noticing that the theoretical reduction factors are not very 





‘ | y= reduction in energy of longitudinal fluctuations 
| s ” ” ” ~  |aEeral * 























Fic. 4. Effect of gauze on isotropic turbulence. 


different from those valid for a simplified type of upstream turbulence in 
which the longitudinal motion has no variation in the direction of the 
stream. It was seen in Part II of this paper that the effect of the gauze 
on such a disturbance is to reduce the longitudinal and lateral disturbances 
independently, the reduction ratios being 
pe {(l1+-a—ak)/(1 x k)}?, v x", (3.31) 
These expressions are nowhere different from the theoretical values by 
more than about 0-07, while the corresponding value of the energy reduc- 
tion factor 4(u-+2v) does not depart from theoretical value by more than 
0-03. Numerical values of the expressions (3.31) are included in Table 3. 
The closeness of the figures is probably to be attributed to the fact that 
the longitudinal intensity of real turbulence is already considerably reduced 
by the time it reaches the gauze, and the longitudinal disturbance exercises 
comparatively little influence on the lateral disturbance in the neighbour- 
hood of the gauze. 
The only published measurements of the effect of wire gauze on turbu- 
lence are those made by Dryden and Schubauer (4). They give the values 
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TABLE 3 


Effect of gauze on isotropic turbulence 











I 2 3 4 5 9 15 

5 $03 | 0°302 0*242 | 0°202 | Ov12I1 0076 

( I tud | intensit } 0:046 ‘020 O’OIg | 0°022 | 0°031 0°031 
I 543 0°345 | O°251 o199 | 0-166 | ovror 0064 
I 7 245 "174 0°139 | o118 | 0°078 | 0°053 
7 25 5 | o'690 0954 0-955 

{ } I °) $| O°007 | OOI5 | 0°031 0'031 
| { } $47 72 201 104 o*I40 o’ogI e I 
) ‘ = > 0167 orl 0062 





fx and »v for four different values of k, as shown in Fig. 4. Their con- 


lusion from the measurements is that the energy of turbulence is 


reduced in the ratio 1/(1-+-/) and it is implied that the energy contributions 
from the component fluctuations are each reduced in this ratio. This is 


pposed to our own result that the lateral and longitudinal intensities of 


turbulence have different values behind the gauze, and indeed the experi- 


nental points shown in Fig. 4 are not inconsistent with our result in this 


respect. The measurements of v are in quite good agreement with the 
theoretical values, although the measured reduction in longitudinal fluc- 
tuation is much less than the theory predicts. The conditions under which 


the measurements were made suggest two possible reasons for this 


liscrepancy. The root mean squares of the three components of turbulence 


were measured in the working-section of a wind-tunnel, both in the absence 


ind in the presence of a gauze of known resistance in the settling-chamber 


of the tunnel. The effect of a contraction is known to be greater on the 


ngitudinal fluctuation and it is not unlikely that this effect is different 


for the different types of turbulence passing through it. In the second 


place, the considerable distance between the gauze and the measuring 
int would allow appreciable decay of the turbulence, and since there 
sa general tendency to isotropy during decay, the longitudinal intensity 
would be relatively increased. It is Dryden and Schubauer’s intention 
to make further measurements of the effect of gauze on isotropic turbu- 


ence in a manner which avoids these difficulties. 


7. Nature of the turbulence downstream from the gauze 
The preceding analysis has shown that the gauze reduces the lateral and 
ngitudinal intensities by different amounts. However, rotational sym- 


etry about the x-axis of the turbulence is preserved through the gauze so 
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that if the turbulence far ahead of the gauze is isotropic it will be axisym- 
metric downstream from the gauze (3). One consequence of this is that 
the usual decay laws known to be valid for isotropic turbulence will not 
necessarily apply accurately to the turbulence downstream from the gauze. 
The three-dimensional energy spectrum functions for the downstream 
turbulence follow immediately from (3.19) and (3.22). If A,(l, m,n) is the 
density in wave-number space of contributions to the mean square of 
turbulent velocities parallel to the x-axis in the turbulence far ahead of 
the gauze, the spectral densities of fluctuations parallel to the (x, y, )-axes 
far downstream are F,, G,, and H,, where 
F(1, m,n) T(B)F, (1, m,n), (3.32) 
G,(1,m,n)+ H,(l, m,n) v*[ G, (1, m, n)+- Ay (1, m, n)|+-B?(I—a?)F,(1, m,n) 
x Fy(m, n,l)+-F,(n, 1, m)| B?(1—o?)F,(1, m,n). 
(3.33 
Most of the useful statistical characteristics of the turbulence can be 
obtained from these spectrum functions.+ 
A quantity of great practical significance is the spectrum of the variation 
of the longitudinal fluctuation along a line in the direction of the x-axis; 
far downstream it is given by 
a 
| | F,(1, m,n) 


Dn—ox 


4.2]? + (1+ a—ak)*p? 
a) PY dmdn. (3.34) 
5) 9 9 
12+ (1+a+k)2p 
Given a knowledge of F,(/,m,n), it is then possible to predict the effect 
of the gauze on any of the parameters describing the longitudinal correla- 
tion between longitudinal velocities. For instance, the square of the radius 
of curvature at the origin of the correlation function is changed according 
to the ratio m 
||| Fy, m, n)L(B)P dldmdn 
ri. 2 fe sd 
\ ——— - e (: . 
“k=0) 


» [ | [ F,(1,m, n)l? dldmdn 


9 9 


Since J(8) has its minimum value when / is small, comparison with (3.2) 
indicates that A? is reduced on passing through the gauze. The paramete! 
\ has not the same direct connexion with energy dissipation as in isotropi 
turbulence, so that it is not possible to deduce that the fractional rate of 
decay of energy is increased. Indeed, forming the spectrum of the longt- 
tudinal fluctuation along a line in the direction of the y-axis shows that 


+ In order to determine G, and H, separately we should have to consider the relation 
between the phases of different Fourier components in the upstream turbulence. 
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the change in the corresponding parameter (A’, say) is given by 


( ( ( F (l,m, n)I(B)m* dldmdn 
(5 ; - ; (3.36) 


k=0 2 2 


wu 1|| Alm, nym? dldmdn 


nd the same reasoning indicates that A’ is increased by the gauze. 

An interesting result is that the change in the ordinary longitudinal 
integral scale L does not depend on the function F,(/,m,n) and is thus 
identical for all types of turbulence. ZL is proportional to the value of the 


5) 


energy spectrum function (3.34) at l 0, divided by u?, so that 
L, uy (l+a—ak)? (1+a—ak)? 
Li.=0 yz (l+a+h) p(l+a- hk)?" 


Values of this factor are shown in Table 3; it is unity when k is zero or 


(3.37) 


very large and has a minimum of zero at (l+a—ak) = 0, i.e. when 

2-76, according to Dryden and Schubauer’s experimental relation 
between a and k. This vanishing of the scale is a direct consequence of 
the result established in Part II that a disturbance to the longitudinal 
velocity which has no variation in the x-direction is entirely removed by 


1 gauze for which (1+-a—ak) = 0. Evidently the effect of gauze on the 
longitudinal double-velocity correlation is to make it negative over a 


considerable range of the distance between the two points. 

Since the three-dimensional energy spectrum downstream from a gauze 
can be predicted when the upstream spectrum is known, it follows that the 
effect of a succession of gauzes of the same or different resistance can be pre- 
dicted. Equation (3.19) holds for each gauze, so that the reduction in the 
longitudinal intensity of (initially) isotropic turbulence owing to 7 succes- 


sive gauzes having the same resistance is 
|} | L2°(B)F(l, m,n) dldmdn 


( ( ( F, (1, m,n) dldmdn 


Equation (3.22) shows that the reduction in lateral intensity due to r 


successive gauzes is related to that due to r—1 gauzes by 


[| [ BI" [—o?)F (l,m, n) dldmdn 
l ) ag | 1 ee ——————— 


‘| F,(l, m,n) dldmdn 


x 


v7*v,-y—h(l+athk)p,+}(1+a—ak)*u,_). (3.39) 
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Numerical values of yw, and v, for r = 2, 3, etc., can be calculated on 
the assumption that (3.25) describes the spectrum function F,(/,m,n) and 
they indicate, in general, that each added gauze (of the same k-value) 
produces a slightly smaller additional reduction in the intensity of turbu- 
lence. As the number of gauzes is increased, the turbulence tends to 
become a wholly lateral motion and the value of v,/v,_, tends to a?. The 
value of p,/u,_, increases with 7, but its value is not very important since 
uw is very small for r > 1 and k > 1. The energy reduction factor for each 
gauze in addition to the first gauze will thus be ga? approximately, 
However, these theoretical results should be used with caution since we 
have taken no account of the dynamical tendency to isotropy in turbulence 
whose intensity is not vanishingly small. There will undoubtedly be some 
transfer of energy from the lateral to the longitudinal motion in the region 
downstream from each gauze, and this will be greater the greater the 
departure from isotropy. 


APPENDIX 


NATIONAL BUREAU OF STANDARDS MEASUREMENTS OF 
LATERAL FORCE ON GAUZES OF ROUND WIRES 


(Communicated by H. L. DRYDEN and G. B. SCHUBAUER) 


A preliminary version of the preceding paper was made available to H. L. Dryden 
and G. B. Schubauer of the National Bureau of Standards, who had published som« 
information on the effect of damping screens in reducing wind-tunnel turbulenc 
(Journal of the Aeronautical Sciences, 14, No. 4, April 1947). Preparations wer 
being made to proceed with additional and more direct measurements of the effect 
of damping screens on isotropic turbulence at the National Bureau of Standards. 
In view of the apparent importance of knowledge of the lateral force on such screens, 
W. Spangenberg of the Bureau staff undertook the measurement of the lateral fore: 
on screens made of round wire. 

The method adopted was based on the formulae given some years ago by Taylor 
for the relation between the deflexion of an air stream impinging on a screen at 
an angle and the lateral force. The measurements were made at the end of a 12-in.- 
square duct discharging into free air, the end of the duct being cut off successively 
at various angles @ and various screens being placed over the end of the duct. 
Precautions were taken to secure a uniform velocity distribution upstream from 
the screen with thin boundary layers. The direction of flow was measured by fin 
silk fibres at ten stations about 1 in. apart in the central region of the stream and 
along a line perpendicular to the flow. In some cases measurements were mad 
upstream of the screen as well as downstream. The angles were read to the nearest 
+ degree by a protractor, using a mirror to avoid parallax. The mean deviation fo! 
ten stations varied from 0-2 to 0-6 degree, differing from screen to screen, but fairly 
consistent for each screen. 

The results of observations of ¢ and kg are shown in summary form in Table 4. 
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d o1 TABLE 4 
) and Summary of values of bd and kp from faired curves 
alu ctieaeaeeii 
54-mesh 40-mesh 24-mesh 
urbu- “ 
1s { d " d ko d ko 
ris { — - 
Th 228 6:9 | 1:60 | 8-4 | 080 
7 1°85 7° 1°39 8-5 | 06g 
since 7 1°65 7°! 1°25 o5 0°04 
Cacti 2°31 9°5 | 165 | 128 ) 
itely 1°8 10°! 1°43 13°70 | 0°76 
i } I *O6 1o°s I*3I 13°00 o*72 
f Wwe - — —— 
lene 1 1-92 | 205 | 142 | 246 | 0°75 
1°5( 21°7 | 1:20 | 25°5 | of 
some 4 | 20°5 1°36 22:2 10g | 25°9 | 0°58 
C2101 2 I 32°4 1°05 37° | O51 
r the 1 | 108 | 33°4 84 | 39°5 | 0°44 
2 19 35°3 74 40°0 O41 
F,/6 | | from the formula 
i 2 sin(@—d)cos 6 
Acosd 
s ce he relatio1 
Po kg. 4pU?. 
Ww —— - 7 
dar o 
fi C - 
Pe s 
rn fe 
c > a . ° “6 a 
z= io a o 6=10 \ Dryden & 
: A e @=15°30°45°f Schubauer 
9. ay a Simmons 
et) | ee a, ee acer Fo _ ake 
- y 6 Btkg 
4 mesitaiain F =? ~-2s 
I 6 ~ vVitke 
\ = | 
: ] 2 ke 3 4 5 
Fic. 5. I rimental values of F,/@ against Kg. 
al 
eps Ther variation of ¢ with the speed (and consequently with kg) for a 
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here is a fairly definite relation between F;/@ and kg 


screen as shown in Fig. 5. 





Suppose the 
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screen to have n wires per foot. Let D be the drag of a foot length of the wire 


flow at the screen is assumed to be equal to 6 and the velocity at the sereen is U, then 


Fy. 3pU? nDsin8@, 


kg.3pU? nDeos6+-nD. 
ko sin@ 
Henc: Fy A" 
n l cos6@ 


But the direction of flow at the screen will probably be intermediate between 6 and 
dé, analogous to the induced angle of attack on a cascade of airfoils producing lift 
As an approximation, assume that the induced deflexion at the screen is equal to 
one-half the final deflexion (9 —¢); i.e. replace 6 in the first-order theory by 6 1(0—4). 
Thus, designating 6—¢ by 26, we have the pair of equations 

2cos@ 


k= — sin 26, 
6 cos(@— 20) 


2 kg sin(@—8 
Ky . 
] cos(G—o) 


On elimination of 6, we obtain a relation between Fy and kg depending slightly on 6. 


Typical values are given in Table 5. 


TABLE 5 
Theoretical relation between Fy and kg 











d Fy F,/0 ko 

6 Is 5 37 I 237 3°85 
7 271 I'035 | 2°30 

204 790 1°94 

II ’ 524 1°21 

13 265 57 

0 = 20 Io I I*I49 z41 
I4 4 49 | 2°53 

IS | 37 723 | 1°78 

22 26 497 I'l3 

at 124 257 54 

6 45 I 771 8x 2°97 
17 | 884 | 2°50 

21 616 754 2°05 

25 534 ‘079 1°69 

: $4° *505 1°33 

33 "351 446 *99 

37 247 314 66 

}1 131 ‘166 "33 





If 8 and ¢ are small, the two equations reduce to 


Fy = 4k(0—8), 

M k 
whence fo =... 

6 8+kg 


forming one side of the square mesh and suppose it to be independent of the angle 
between the axis of the cylindrical wire and the wind direction. If the direction of 
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This simple equation fits th xperimental data very closely for kg 1-4. The 
leparture at larger values of kg is not surprising in view of the probable interactions 


tween the individual wires. 


Details of these measurements will be published by the staff of the National 


jureau of Standards in conjunction with further data on the effect of screens on 
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THEORY OF A LOOP REVOLVING IN AIR, WITH 
OBSERVATIONS ON THE SKIN-FRICTION 


By C. C. L. GREGORY 


(University of London Observatory, Mill Hill Park, N.W.7) 
[Received 22 December 1947] 


SUMMARY 

An investigation is made of the dynamics of a flexible loop, driven by a small 
pulley at a constant speed, for the case when the speed is sufficient to cause the loop 
to become air-borne. It does this when the air-friction is greater than the weight 
of the loop. The shape and position of the loop may be calculated for values of th: 
ratio, 2a, of air-friction to weight and the angle, d*, that the tangent to the ascending 
portion of the loop at the pulley makes with the horizontal. This angle may bi 
controlled by means of a second pulley so arranged as to press the loop against the 
driven pulley. If « 1, theory predicts a point of the loop at which its radius of 
curvature becomes zero. For values of « between } and 1 the motion is stable, and 
the precise value of « may be found from the photographed shape of the loop. As 
the result of some trials, values of a have been obtained which have given a valu 


of the air-friction drag coefficient of the right order of magnitude. 


Introduction 

It has long been known to engineers that a belt when driven at a high speed 
will occasionally be seen to stand up in the air in the form of a loop. This 
may happen when a shaft breaks, or when a belt slips off a driven pulley 
at the end of a shaft. The dynamics of such a revolving loop does not 
appear to have received previous consideration, nor does the more general 
problem of an endless flexible cord subjected to external forces which do 
not all pass through a fixed point. 

The author’s attention was drawn to this problem at the Admiralty 
when examining a proposal to use such a loop as a means of defence 
against aircraft. At that time no very precise values of the air-friction at 
a moving surface were available, and it occurred to the author that the 
revolving loop, found to be valueless for the purpose proposed, might 
prove useful as a means of obtaining this information. Accordingly, the 
theory of the revolving loop has now been rewritten. Some tests of the 
suitability of the loop for measurement of surface air-friction were made 
at the University of London Observatory in February 1947 when astro- 
nomical observations were interrupted by bad weather. 


Part I. Theory of the revolving loop 
Let the are APO (Fig. 1) be the locus of any point, P, of a uniform, 
flexible, inextensible string of weight mg per unit length, which is moving 


(Quart. Journ. Mech. and Applied Math., Vol. II, Pt. 1 (1949)] 
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THEORY OF A LOOP REVOLVING IN AIR 31 
with constant tangential velocity V. Suppose the string is subject to the 
following external forces: 
(i) its weight acting vertically downwards, 
(ii) a constant tangential retarding force 2amg 4s on every element of 
length 8s, 
(iii) a constant tension, 7), at O, where the tangent to the string is 
horizontal, 
iv) a constant tension, 7'*, at A, the tangent at A making an angle ¢* 
with the horizontal. 














Fic. 1. 


For uniform motion the forces (i) to (iv) must balance the centripetal 
force (mV2 5s)/R, where R is the radius of curvature at P. If 7 is the 
tension at P, and the tangent at P makes an angle ¢ with the horizontal, 


then 


éT mg(2x-+-sin ¢) ds, 
sas (¥* ’ 
T’' d¢ mg\ — cos} ds, 
gh 


where the motion is in the direction A PO, and s is measured from O along 
the curve of the string to P. The differential equations to be solved are 


thus 
dT : 
: mg(2a-+sin ¢), (a) 
mV2—T mg cos ¢, (b) (1) 
where R = (c) 
aod 
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From 1(6) and l(c), 


iT iR ‘ 

i mg — cos¢+ mg sin d, 

ds ds : 

dR ; 
so that, by l(a), * = cos¢@ = 2(a-+sind), 
iR 

or, a = 2(asec¢-+-tan ¢) dd. 
Hence, we obtain k = k{tan({7+-4¢)}°* sec’, 


where / is a constant of integration, and ee 1(b), 
mV?—T = myk{tan(ia7+-4¢)}?sec d. 
Consideration of the conditions at ¢ = 0 leads finally to the expression 
mV?—T = (mV?—T,){tan(ia7+-4¢)}?% sec ¢d. (2 
Since change of sign of « corresponds to ane of the direction of J, 
the above equation, when written in the form 
mV?—T' = (mV?—T,){tan(47-+- $¢)}-**sec ¢, 
represents, for negative values of ¢, the ‘eiiietii branch of the curve 
from O to the po C,.where the tangent is vertical. 
Writing ¢ for tan(47-+-3¢) and K for (V?/g—T,/mqg), (2) becomes 
mV2—T = 4Kmg(t?**1+- 724-1), (3) 


Since T' > J at C, it follows that 


oa > 1. (4) 
We also have, from (1) 
R = K{tan(47+-4¢)}**sec*d = }K (t?%-2-+ 2424+ 72072), (5) 
and since ds/dt = Rdd/dt = 4.K(t?*-?-+-#%), 
it follows that s = $K[t?°/(2a—1)+#%+1/(2a+1)]. (6) 


Taking axes through C, x to the left and y downwards, the coordinates 
of P are 

2 = 4KF4*/a, y = $K[P%1/(204+1)—P*-1/(2a0- 1)]. (1) 

We consider now the equations of motion of an element 5s’ of a similar 

string at the point P’ and moving with the same velocity, V, the tangent 

at C being vertical and the tension at C being mV? (Fig. 2). These 


equations are 


87" —mg(2a—sin ¢’) ds 


T” 54’ - ma oe + cos $'] 5s’ 
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Using a similar method to that given above we obtain 
T’—mv? 4 K’mg(t'2*+1 + ¢/24-1) (9) 
where ¢’ is the angle which the tangent at P’ makes with the horizontal, 
K’ = (T*'—mV2)/4mg(t#"2%+1-+ #2071), 
t’ = tan(i7—}¢’), t*’ = tan(j7r—3¢*’). 
T*’ and ¢4*’ are, respectively, the string tension and the angle of the 
tangent at a point, B, at which the constant driving force, 7'*’, is applied 








to the string. Accordingly, P’ lies on a curve similar to the portion, OC, 
if the curve in Fig. 1, except that it is inverted and the scale is increased 
by the factor K’/K. 

We may now suppose that the two portions of string form one con- 
tinuous loop COPA P’C, by joining the curves P and P’ at C and making 
B coincide with A. [In practice this may be done, approximately, by using 

large loop constrained to bend around a small driving pulley at A.| 
Equating the w-coordinates of A and B, we have 

K’/K = (¢*/«*"}**, (10) 
where t* is the value of t at A. Equating the y-coordinates of A and B, 
ind using (10), we obtain 


t*t* (2a+1)/(2a—1); (11) 


(12) 









ind, from (6 S arc COPA 1K 
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Equations (10) to (13) give 


S = 8S’, or arcCOPA = arcCP'A. (14 
From (9), 


T*’—mV? = 4K'mg(t*' +t*’-1)t*’24 
= $Kmg| (2a—1)t*?*+1/(2a+1)+(2a+ 1)t*224 (2~—1)| 
and from (3), mV2—T* = 4Kmg(t* +t*-1)t*2%, 
whence a tamgs, (15) 
and also 


T*' + 7* = 2mV2—K[t*/(2.+1)—t*-1/(2a—1) tg 





aie, (16 
2mV2—2y*mgq, } 
where x*,y* are the coordinatés of A. Accordingly, 
T*’ = mV?+-2amgS—mgy* (17 
and f lag mV?—2amgS—mgy*. (18 
Also, if x’, y’ be the coordinates of P’, 
x’ = 4K't'24/q, (19 
, 1 ¢ 2x—1 f'2a+1 
y = 3K — (20 
7 2a 2a+1 
The area of the loop is given by 
A A 
A y’ dx’ — | y dx 
O O 
gee ger te 
r e f4a-2 eff 4a P 4a . fa—2 
Lk’ dt ~— | dt|—_ 1k? | t dt t dt| 
J 2a—1 J 20 2x 2a—1 | 
0 0 0 0 7. 
=| t ge-1 t 
4 L(2a+1)(4a—1)  (2a—1)(4a+1)  (2a+1)(40+1) 
t* 1 | 
(2x—1)(4a—1) | 
2a S2* 
(21 
16a0?— ] 


The left-handed couple about A, required to balance the weight of 4 
and S’ together (considered as a rigid, stationary loop held at the point A) 


is given by 
rm 1 A 


G mg} | (~x*—x)ds+ | (a*—2’) ds’ 
Oo 0 
1 A 
2mgx*S—mg | x«ds—mg | x’ ds’. 


Oo O 
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Since 
Kt /2a, a K't'**/2a, ds LK (2%+-029-*)dt, ds’ 4K’ (t'2*+-¢'29-2) dt, 


t follows from (21) that 
160° , 
G aa} mg Sa* tamgaA. (22) 
Here G is the left-handed couple about A, arising from air skin-friction 
when the loop is rotating with constant speed V. Air-friction will exactly 
balance gravity, and the loop becomes airborne provided «a > $, as 
required by (4) 
In order to determine whether the radius of curvature, R, of the loop 


ilwavs remains finite, we have, from (5) 


IRidt EK (a | ¢2%+1 4 Qayt?—-1 + (a — 1) #7 " 
] \ 
0 when ? 0 (when 2a a), l, or ; 
Ito 


i 
| \ 


respectively. The corresponding values of R, given by 


rives values of é from —30° to —90° for values of « from } to 1 


| K(1 x)*—-1(1 + ax)-2-1, 
we finite for a lL. For lL. Ruin = 0 at C, where —¢ = 3a. [In an 
ctual case some distortion of the loop occurs when a > 1, since no 
material loop is infinitely flexible; the rotational effect of air-friction, in 


this case, causes the top of the loop to bend back somewhat. | 

From (21) it will be seen that, when V (and consequently a) is constant, 
the area of the loop, A, varies as x*; the locus of A when C is supposed 
fixed may be found as follows. 

Let the two tangents at A meet the vertical through C in £ and D 


Fig 3 the 


AK v* sec b* Kt*2%(¢* +-#*-1) /4a 
mV?— T'*)/2amg (S+y*)/2a, 
from (18 AD = x* sec d*’ 
K 120] = l , ~ Lis 7 
| 2a—] 2a+-1 
Thus 1H+4AD Kt au (20-4 1 )-1¢* +- (2a— 1)—"¢* | 
2 constant, 
1D S—y*/20 
CE v* tan d* —y* S/2a. 


v* tan d*’ 
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hence the locus of A is the ellipse 


‘ai - .. - 
y** + 28 = §* 
4y*— ] 
D and E being the foci, C the centre, and S and S(4a?—1)*/2a the 
semi-axes. 





? (23) 








Part II. The airborne loop in practice, and the determination of 

skin-friction at various speeds 

The determination of air-resistance to a surface, when in uniform motion 
relative to still air parallel to the surface, has always proved difficult from 
an experimental point of view on account of end-effects of the material 
body which provides the moving surface. An endless cord or belt of 
uniform cross-section driven by a pair of wheels mounted on fixed parallel 
axes, and made to revolve at a constant known speed would, indeed, 
overcome the difficulty of ‘end effects’, but other difficulties would arise. 
Only at very low speeds would such a driven belt be free from ‘flutter’ 
unless it were stretched very tight, in which case the presence of eddy 
currents might be concealed. Another difficulty with such an arrangement 
would be the accurate determination of the rather small additional force 
required to drive the belt against air resistance, over and above that 
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required to overcome the friction and air resistance of the driving wheels 
and other moving parts of the apparatus. 

When a belt is driven by a single small wheel or pulley, however, these 
difficulties do not arise, the motion of the ‘loop’ being observably steady 
so long as the speed is such that the value of « lies between } and 1. The 
precise value of a can be found from the curvature of the loop in the 








Fic. 4. 


vicinity of the point C (Fig. 1), no measurement of force being required. 
Both the shape of the loop (giving «) and its rotational speed, V, may be 
determined by a cinematograph camera working at a known speed, a short 
length of the (white) loop being blackened so as to be visible. 

From the foregoing theoretical discussion the necessary and sufficient 
conditions for the maintenance of an air-borne loop in a steady state can 
readily be seen. 

Let a (Fig. 4) be a cylindrical pulley wheel driven at a constant peripheral 
speed, V’, by an electric motor, say. Around a is placed the loop of 
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length 2S and weight 2Smg. A loose pulley, b, on an axle parallel to that 
of a, keeps the loop in contact with a at A. 

Let QA be the common tangent of a and b at A, and let QA make an 
angle ¢* with the plane of the horizon. On first starting the motor, the 
part of the loop in the vicinity of A is projected along this tangent QA 
with velocity V, and, if air resistance is neglected, would attain a vertical 
height equal to V?/2g. Taking air-resistance into account, V must be 
sufficiently great to make 7'* greater than zero in (18). 

Once the loop is airborne, the whole family of possible curves given by 


(23) for all positions of A in Fig. 3, and subject to the restriction d* > 90°, 
may be obtained by rotating the line of centres of a and b, and thus alter- 
ing d*. However large V may be, the loop will fall as soon as 4* = 90°, 


as may be seen from (2). 

Should it be possible to measure with accuracy the constant couple 
which must be applied to a to maintain the loop at any (sufficient) speed, 
v would at once be obtained from (15), the mass (2Smg) of the loop being 
readily measured. This value of « could then be compared with that 
derived from the radius of curvature of the loop at C, as mentioned in the 
second paragraph. 

It is worth noting that measurements of surface air-friction over a con- 
siderable range of speed can be made by employing loops of very different 
cross-section, V being greatest (for «a = 0-6, say) when the cross-section is 
circular. If the cross-section be rectangular, of thickness a and breadth } 
much greater than a, then 

vmg = cVi(a+b), 
where c is assumed constant. 

If y be the radius of the circular section of a loop of similar material 
such that ar? = ab, the ratio of the speed, ), which it must have to give 


the same a, to JV, is given by 


(zb/a)! (approx.) 


tr * 


. (a+b) 


With b/a = 800, the ratio of the speeds is about seven. 

A few trials were made at the University of London Observatory in 
February 1947. Loops from 7 to 20 feet in length were maintained in the 
air by an electric motor driving pulleys of 0-298 and 0-255 ft. diameter 
respectively at about 50 r.p.s. The loops were made from strips of chrono- 
graph paper 0-0833 ft. wide and 0-00031 ft. thick, m being 0-0012 Ib. /ft. 
. was found to vary from about 0-5 to 0-7 with speeds of about 40 and 46 
ft./sec. and was obtained from photographs of the loop by measuring Fyjn- 
It was found that the loop could be kept steady for any length of time so 
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o that mg as it was subjected to a slight draught of air perpendicularly to its 


plane. This was probably necessary in order to prevent circulatory air 


— wrents forming and so reducing the air-speed of the loop. When driven 
» a t constant speed in still air the loop was observed slowly to drop a little 
it OA nd then more quickly recover its original height, this ‘hunting’ continuing 
t 
ti na somewhat irregular manner. 
ist. by Regarding the accepted value of the skin-friction drag coefficient K,, 
the value suggested by M.A.P. for a steel wire is between 0-003 and 0-0015. 
en by \ssuming K, pV*(a-+-b) xm, 
() here p 0-00238. a 0-00031,. 6 0-0833, and m 0-:0012. we deduce 
ute] that , "9 
‘ kK . 0-166V2, 
{) 
hence, i) fora 0-5 and V? 1600, we obtain K, = 0-0019, 


2100, we obtain K, 0-0020. 




















THE PLASTIC YIELDING OF NOTCHED BARS UNDER 
TENSION 


By R. HILL (Cavendish Laboratory, Cambridge) 


[Received 3 February 1948] 


SUMMARY 


The state of stress in the core of a notched bar is analysed at the moment when 
pronounced plastic yielding begins. The theory is two-dimensional and an idea 
plastic-rigid material is assumed. Following a general analysis of the problem, tli 
magnitude of the constraint factor is calculated for a deep notch with a semi 
circular root. The correct approach to problems of plane plastic strain is discussed, 
and is further illustrated by a re-examination of the classical work of Prandtl on 
indentation (which is closely related in principle with the notched bar problem). 


1. General considerations 
THE paper is concerned with the state of stress in the core of a notched 
bar under an applied load sufficient to cause plastic flow. It is known 
that the fairly well-defined bend in the load-extension relation for a 
notched bar corresponds to a mean stress in the core that is greater than 
the yield stress in a tensile test of an unnotched specimen. The investiga- 
tion will be directed particularly to a calculation of the magnitude of this 
effect in relation to the shape of the notch. The problem will be treated 
as one of plane strain, so that the specimen is a wide rectangular block 
symmetrically notched through the thickness on two sides. The notch is 
shown in Fig. 1, the plane of the paper representing a plane of flow. Let 
it be supposed that the applied load at the ends of the specimen is staticall) 
equivalent to a tension along the central longitudinal axis. It is also 
supposed that a distribution of stress, suitable to ensure plane strain, is 
applied to the two lateral surfaces of the block which are parallel to the 
planes of flow. The specimen is assumed so long that the state of stress 
in the neighbourhood of the notch is, to any desired approximation 
independent of the precise distribution of the end load. It is further 
assumed that the notch is of sufficient depth to ensure that the plasti 
region is localized in and around the core (this will be discussed more 
carefully later). 

If the curvature of the notch is greatest at the root, the yield limit is 
first attained at that point. As the applied end-load is further increased. 
the plastic region can be expected to spread in the manner shown in 
Fig. 1. If the constant maximum shear-stress criterion of Tresca is adopted 
(the work-hardening being supposed zero), the problem of calculating the 
[Quart. Journ. Mech. and Applied Math., Vol. II, Pt. 1 (1949)] 
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stress components in the plane is statically determined (the stress normal 
tothe plane is the intermediate stress). While part of the core is still elastic, 
the strains are small enough for the changes in the external surfaces to 


e neglected. The stress boundary conditions, known a priori, are then 


sufficient to determine the distribution of stress (in the plane) without 














Plastic 


regions 

















It Che notched tensile specimen. A ty pical plane of flow. 


reference to the velocities of flow.+ Since the elastic compatibility equa- 
tion does not involve elastic constants (there being no body forces), the 
entire distribution (in the plane) is independent of the elastic constants. 
In the elastic region the stress normal to the planes of motion is equal 
to vXthe sum of the other two principal stresses (where v is Poisson’s 
ratio); in the plastic region it can only be found by following the history 
of the deformation, using Reuss’s equations.{ The plastic-elastic boundary 
is a curve along which the maximum shear stress is constant (not a slip- 
line). It seems that the stress can only be universally calculated by 
numerical methods. The relaxation technique has been applied by South- 
vell and Allen (14) in determining the plastic and elastic stress in a 

Strictly, the velocity distribution at each moment should be examined to see whether 

® plastic work done on each element is positive ; otherwise unloading takes place. 
These 1 st conveniently be found in ref. (3). 
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notched bar with a shallow notch, but no calculation appears to have been 


within 

made for the deep notch considered here. aon 
In the absence of such a calculation we can proceed further by intro- appli 
ducing the assumption of a plastic-rigid material, which is rigid up toa —_ 
well-defined yield point and does not work-harden. This does not affect 19). 


the stress distribution in the plane, nor the shape of the plastic region, 
corresponding to a given end-load. In the plastic region the stress normal 
to the plane is now equal to the mean of the principal stresses in the 
plane; hence there is a discontinuity in the stress normal to the plane 
(unless » = 4) along the line of separation between rigid material and 
material undergoing deformation. So long as part of the core is not plastic 
no displacement of the ends of the specimen is possible. However, even 
when the plastic region has spread right across the core it still does not 
necessarily follow that deformation is possible, since the incipiently plasti: 
material may be rigidly constrained by the adjacent bulk of non-plastic 
material. To state this more precisely it is necessary to introduce the 
equations of Saint-Venant governing the flow velocities in the plastic 
region. These are simply (i) that the volume change is zero, and (ii) that 
the principal axes of the rate of strain coincide with the principal axes of 
stress. The characteristics of the velocity equations are the slip-lines 
(curves whose directions coincide with the directions of maximum sheat 
stress). A little consideration shows that no extension of the specimen 
is therefore possible until the plastic region has spread sufficiently far to 
include the entire slip-lines OS, OS’ from the centre O to the free surface 





(Fig. 2).¢ This is the mathematical expression of the reason normally 
offered to account for the raising of the apparent yield stress by a notch: 
namely, the constraint on the possible deformation of plastic material 
exerted by the adjacent bulk of elastic material. At the moment when 
extension first occurs, the plastic boundary must touch OS, OS’ as in (a 


ire 

or must pass either through S and S’ as in (b), or through O as in (c). It wh 
does not seem possible to say which configuration actually occurs without -_ 
the full solution; however, this is not necessary for the present purpose.2 wd 
Since both velocity components are known on the plastic boundary in * 
terms of the (unit) speed of extension, the velocity solution is uniquely Sy 
determined in regions SOB, S’OB’ by virtue of the properties of charac- 63 
teristics. Hence the material within SOB, S’OB’, though plastic, moves i 
as a rigid body attached to the non-plastic material. The velocity solution ¥ 
The initial assumption of a sufficiently deep notch is used here; with shallow notches ay 

the plastic region will be quite different. J; th 
+t It appears from calculations of Southwell and Allen that (b) is probably the actual by 


configuration (private communication). The calculations have so far been carried out fo! 
shallow notches only, and for an intermediate range of loads. 
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rithin SOS’ is uniquely determined by the known velocity components 
rmal to the slip-lines SO, S’O. It can be found, if required, by the 
pplication of Riemann’s method (see Geiringer, 1). Across SO, S’O there 
sa discontinuity in the tangential component of velocity (constant saltus 


Such discontinuities in tangential velocity across slip-lines are of 




















(b) (c) 


Fic. 2. The three possible modes of development of the plastic region. 


frequent occurrence in plane-strain deformation of a plastic-rigid material] 
vhich does not work-harden. They represent an extreme case of the 
nalytie possibility of arbitrarily prescribing the normal components of 
velocity along two intersecting slip-lines of a given field. The discontinuity, 
be it noted, is only momentary in each element crossing the slip-line. 
Such an element abruptly changes its direction of motion, and undergoes 
\ finite instantaneous change of strain. In a real metal the discontinuity 
vould be diffused, partly by elastic strain-increments and partly by work- 
lardening, into a region where the rate of shear is large. The closer the 
ipproximation effected by the plastic-rigid material, the narrower will 
this region be. Also, in a real metal. regions BOS and B’OS’ would not 


be } 


gid, but would be undergoing plastic strains of the same order as 


elastic strains 
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We now examine how the load can be calculated at the moment when 
extension begins. Since the slip-lines are also characteristics of the stress 
equations in the plastic region (Hencky, 2), the slip-line field and stregges 
within SOS’ are uniquely determined by the shape of the stress-free 
surface SS’.+ Furthermore, they are independent of the notch depth 
provided this is sufficiently large; the plastic boundary will of course vary 
with notch depth. The load and distribution of stress across the narrowest 
section can thus be found at the beginning of the extension without a 
precise knowledge of the shape of the plastic region or of the elastic stress 
distribution. It is clear that this fortunate possibility is due to the coin- 
cidence in pairs of the characteristics for the stresses and velocities. The 
possibility would not necessarily occur for a material yielding according 
to Mohr’s criterion (9). The four characteristics are then distinct. The 
stress characteristics are in the directions corresponding to the points of 
contact of a Mohr circle of stress with the yield envelope (Mandel, 7 and 
8); the velocity characteristics are still the slip-lines. It must be empha- 
sized, too, that the assumed manner of the spreading of the plastic region 
for a deep notch, while almost certainly correct, has yet to be confirmed 
by a full solution. Similar methods cannot be used for shallow notches, 
since the general shape of the plastic region can scarcely be surmised in 
advance. For this reason the present considerations are limited to notches 
which are deep in comparison with the core diameter. 

A ealculation of the variation of the load with increasing extension is 
outside the scope of the present investigation. It is merely necessary to 
remark that, since the core width decreases during the deformation, the 
load needed to maintain flow depends only on the changing contour of 
the notch. Hence the load-extension curve for the plastic-rigid material 
has a sharp corner at the beginning of extension. In a real metal with a 
sharp yield point and a small rate of hardening, this ‘corner’ is rounded 
off by the influence of elastic strains, but is sufficiently well defined to 
allow an approximate estimate of the constraint factor. Plastic strains 
of elastic order occur, of course, before the corner is reached, following 
the initiation of plasticity at the root of the notch for a much smaller 
load. When the metal work-hardens rapidly, any basis of comparison 
between notched and unnotched specimens must be somewhat arbitrary 
in view of the completely rounded load-extension curves and the non- 
uniformity of the deformation. However, the constraining effect of the 
notch is still apparent from the general raising of the load-extension curve. 

In the above discussion the problem was simplified by the use of Tresca’s 


The stresses are also uniquely determined in terms of the notch contour within a part, 
at least, of regions BOS, B’OS’. 
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vield criterion. If the Huber-Mises criterion is used, certain differences 
must be noted. For the plastic-rigid material the yield criterion reduces 
to Tresca’s with a modified yield stress 2Y/v3, since the stress normal to 
the plane is equal to the mean of the other two principal stresses. The 
jip-line field within SOS’ is therefore unchanged, and all stresses there 
we simply multiplied by the factor 2/V¥3. The plastic boundary is, how- 
ever, different unless v = 4, since the yield criterion now involves the 
stress normal to the plane and this is discontinuous across the plastic 
yundary. Hence the stress component acting normally to a plane element 
perpendicular to the boundary is also discontinuous unless v = 4. For the 
same reason the stress distribution in the plane before extension begins 
snot everywhere independent of elastic constants as before, but depends 
n the value of v. With a general plastic-elastic material the stress problem 
sno longer statically determined, since the stress normal to the plane in 
the plastic region depends on the previous strain-history. The stress 
listribution in the plane should nevertheless be closely approximated by 


that for the plastic-rigid material. 


2. Notch with a semicircular root 

Suppose the notch is parallel-sided with a semicircular root, centre C 
ind radius r (Fig. 3, showing half of specimen only). The slip-lines in the 
field SOS’ defined by the stress-free semicircular surface are logarithmic 
spirals; the trajectories of principal stress are the normals to the surface 


ind concentric circular ares. If 2R is the width of the narrowest section, 


the angle SCS’ = 26 (Fig. 3a) is given by the relation 
R , 
e’—1 (0<0< $n). (1) 
- 


This slip line field is valid in the range 

i) ei7— ] 3°81 approx. (2) 

- 

For larger values of R/r the field is shown in Fig. 3b. The regions SQR 
ind S’Q’R’, containing straight slip-lines, are uniquely defined (in type, 
ut not extent) by the parallel free sides of the notch. Region RT'R’ is 
lefined, as before, by the semicircular surface, and CT’ = e!7r. Regions 
PQRT and P’Q'R'T are then uniquely determined by the base slip-lines 
RQ, RT. and R’'Q’, R'T. The normals to TR, TR’, respectively, constitute 
one family of slip-lines, the other being the orthogonal trajectories of these 
normals. Furthermore, 7’P, QR, TP’, QR’ are all equal. Finally, region 
'PP'O with straight slip-lines is uniquely defined by the base lines 7’P 
ind TP’. The rigid part of the plastic region lying outside SOS’, at 


present unknown, is not indicated in the figures. 
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The stress can now be calculated by integration of Hencky’s well-known 
equations 


dp+2kdé = onan a-line, | 
. (3) 
dp—2kdé = 0 ona f-line, | 
where p is the mean compressive stress in the plane (i.e. —4(o,+-0,), where 


bai 


<— 2 — 

















+. 3. The slip-line field for a circular notch. 


o, and o, are the principal stresses in the plane), k is the maximum sheat 
stress, and ¢ is the anticlockwise angular orientation of a slip-line to some 
fixed direction. The « and f families of slip-lines are distinguished (Fig. 3 
by the sense of the shear stress acting across them. Of most interest for 
the interpretation of notched-bar tests is the distribution of the axial 
stress o across the narrowest section. On the free surface of the plastic 
region p hk, and so, in the configuration of Fig. 3a, o is distributed 
according to the formula 

oO 

2k 


where x is distance measured from the root. Hence o rises steadily from the 


ee log(1 J (0 < z= 3), (4 


value 2k at the surface. Plastic flow has thus dissipated the stress concentra 
tion present when the material is still elastic. The transverse stress in the 


neck is always o—2k, and is a tension. By integration the load is 


L 2Re ueR(L 4 p)loe( 1 + ) (= <e*—1), (5) 


1 


/ 


r } 








where 
const! 
from — 
distri 
and e€ 
sufficl 


it has 


Th 
is th 
ment 
point 
facte 
Kun 
Whe 
tang 


obvi 


al 











-knowr 





sneal 


Some 











PLASTIC YIELDING OF NOTCHED BARS UNDER TENSION 47 
chere ¢ is the mean value of the axial stress, and @/2k is, by definition, the 


nstraint factor measured in notched-bar tests. This factor rises steadily 
a] 


fom the value unity with increasing R/r. In Fig. 36 the axial stress is 


listributed up to 7’ according to equation (4). Along 7'O, o is constant 


nd equal to 2h4(1+-37). It is clear that this will always be the stress 


sufficiently far from the root, whatever the shape of the notch, provided 


t has parallel sides and is not re-entrant. The constraint factor is 

ed fh hee A (e4™ —1—47). (6) 
2k ° R _ 
The maximum constraint factor obtainable with a parallel-sided notch 
s thus 1+-}7, or about 2-57. This is close to the value found experi- 
entally for cylindrical notched bars, when the material has a sharp yield 
int. Thus Orowan, Nye, and Cairns (11) report a maximum constraint 
factor 2-6 for annealed mild steel. The significance of this in relation to 
Kuntze’s so-called ‘cohesive strength’ has been discussed by Orowan (10). 
When the notch is not parallel-sided but wedge-shaped (the sides being 
tangential to the circular root), the modification of the above solution is 
bvious. If 2a is the wedge angle, the maximum constraint factor is 
l7—x, corresponding to vanishingly small root radius. 

Certain properties of the velocity solution are worth noting for the 
ld of Fig. 3b. If u, v are the velocity components along the a, f slip- 

nes respectively, Geiringer’s equations (1) are 
du—vdd = 0 along an a slip-line, | (7) 

dv+u dd 0 along a f slip-line. } 

the speed of the extension is U’. the material within OPP’T moves as 


rigid body with speed U in the direction 70. In the right-hand half 
f the field (say) the velocity component wu is constant on each a-line 
thin P’Q’ R'T. Hence since v is constant on 7'P’ (and equal to — LU’ v2), 
is also constant on each a-line. In particular v is equal to —3U’/v2 on 
RR’. The material within Q’ R’S’ moves, therefore, as a rigid body with 
component speeds U’, 2U’, parallel and perpendicular to the side S’R’. 
rhe effect of the plastic deformation is thus towards the formation of a 
ey-hole it the bottom of the notch. The velocity within TRR’ can be 
letermined by applying Riemann’s formula. 

Van Iterson (6) has inde pendently proposed various slip-line fields in 
the neighbourhood of wedge-shaped notches when the root radius is 
vanishingly small, but without discriminating between deep and shallow 


otches. He did not discuss the inter-relation of the solutions in the plastic 


ind elastic regions, and was evidently unaware of the need to consider 


le associated velocity solution in deciding when extension is first possible. 
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3. Relation to the problem of indentation with a flat punch 
When the strains are so small that changes in external surfaces can be 
neglected, the problem of the indenting of a semi-infinite mass by a flat 
rigid punch is clearly identical (apart from changes of sign) with that of 
a notched tensile specimen when the notch is infinitely deep and sharp, 
The boundary conditions on the punch face RR’ (Fig. 4), namely, zero 











(2) (b) 


Fic. 4. Comparison of present solution (a) with Prandtl’s solution (6) 


for the indentation problem. 


shear stress and constant normal velocity, are just those on the narrowest 
section of the core in the notched bar. Again, the surface of the mass not 
touching the punch is stress-free, as are the sides of the notch. This 
indentation problem, now classical in the development of plastic theory 
was first considered in 1920 in a highly original paper by Prandtl (12). 
Although Prandtl’s work preceded by three years the general equations 
of the slip-line field obtained by Hencky, he recognized the hyperbolic 
character of the stress equations and derived the two special fields needed 
in the solution. However, Prandtl seems not to have fully appreciated 
the nature of the problem and the underlying assumptions, with the result 
that his solution is not entirely accurate in detail. It has not been corrected 
by later writers. It is necessary, therefore, to re-examine the problem 
carefully in the light of the previous discussion. 

A smooth flat punch is pressed into the plane surface of a semi-infinite 
mass of plastic-rigid material (Fig. 4a). Indentation cannot begin until 
the plastic region has spread sufficiently far to include the entire slip- 
lines OPQS and OP’Q'S’ from the centre O of the punch face.t The 
curved plastic boundary SS’ can be expected to be of the form indicated 


in Fig. 4, which seems the most likely of the three possibilities of Fig. 2. 


At any earlier stage the slip-lines through the extreme points S and S’ of the incipient!) 
plastic surface do not meet on the punch face. 





It wo 


The si 


of the 
half o 
straig! 
the cc 
is unl 
tinuit 
still h 
with | 
r 
vield 
load, 
load- 
tion 
Wher 
vield 
ways 
of th 
0 §] 
struc 
the ¢ 
relat 
OL SU 
W 
mad 
toh 
only 
elas 
dete 
shor 
the 
reas 


equ 


Ift 
ind 
fg 


mo 





can be 
7 a flat 
that of 


shary 


Vv. Ze] 


itio S 


ecte 














PLASTIC YIELDING OF NOTCHED BARS UNDER TENSION 49 


t would be interesting to have this checked by detailed calculations. 
fhe slip-line field within OPQSR is uniquely determined by the shape 
fthe free surface and the stress conditions along it. In the left-hand 
alf of the diagram, for example, the slip-lines in OPR and QRS are 


straight lines meeting the surface in 45°, while in PQR they are radii from 
the corner and circular ares. The pressure at the beginning of indentation 


s uniform and equal to 2/(1 7). There is a tangential velocity discon- 


tinuity across OPQS (and OP’Q'S’); the plastic material on one side is 


till held rigidly by the non-plastic material, while on the other it moves 


ith constant speed v2 l’ along the slip-lines of the same family as OPQS 
speed of the punch). In a real plastic-elastic material with a sharp 
eld point, irrecoverable indentation occurs at the first application of 
id, due to plastic flow around the corners of the punch. However, the 
,d-penetration curve should have a fairly well-defined bend, the penetra- 
m prior to this being restricted by the constraint of elastic material. 
When the material work-hardens rapidly and does not possess a sharp 
ield point, it is observed that the displaced material does not flow side- 
ys towards a raised lip, but is accommodated by the compressibility 
{the whole mass. It seems that rapid hardening causes the plastic region 
spread downwards rather than sideways. A balance is, as it were, 
struck between (i) flow out to the surface with hardening accompanying 
the associated severe strains, and (ii) a downward displacement requiring 
elatively small plastic strains and hardening, but against the constraint 
f surrounding elastic material. 
We can now compare the present solution with Prandtl’s. Prandtl 
ule no mention of the rigid part of the plastic region, and seems not 
have realized the significance of the fact that the plastic boundary can 
nly be determined by a calculation of the stresses in both plastic and 
istic regions. He was also apparently not aware of the condition 
letermining the start of indentation: namely, that the plastic region 
should just include the slip-lines OPQS, OP’Q'S’ from the centre O of 
the punch face to the free surface. Nor was he apparently aware of the 
eason for this, namely, that the slip-lines are characteristics of the velocity 
lations. In Prandtl’s solution (Fig. 44) the extreme slip-lines CBA R’, 
B'AR pass, not through the centre of the face, but through the corners. 
ithe present argument be correct, plastic flow has not spread so far before 
indentation begins; the velocity is. in fact. indeterminate in the con- 
guration of Fig. 46. Prandtl arbitrarily assumed the region RAR’ to 
ve downwards as a rigid body attached to the punch. The indenta- 
i pressure corresponding to this slip-line field is, however, still 24(1+ $7). 
these objections to Prandtl’s solution have not been noted by other 
E 
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writers. It should be mentioned that it may be difficult to distinguish jp 
an experiment the moment corresponding to the present solution, not 
only because of the influence of elastic strains, but also because thy 
numerical work of Southwell and Allen for shallow notches shows that 
the plastic region develops very rapidly at this stage. 

Prandtl also considered the problem when yielding depends on hydro 
static pressure, and took an arbitrary yield envelope in accordance with 
Mohr’s suggestion. His solution for this case is open to a further objection 
Following Mohr’s formulation of his yield criterion (9), Prandtl assumed 
that, after yielding, the deformation consisted of simple shears in th 
direction of the so-called surfaces of sliding. These correspond to th 
points of contact of a Mohr circle with the envelope and are, as has been 
mentioned, the characteristic directions for the stress equations. In view 
of Mohr’s formulation it was natural for Prandtl to take the plasti 
boundary and the streamlines to be surfaces of sliding. However, for an 
isotropic material the principal axes of stress and strain-rate must coincid 
(elastic strains being neglected). This condition is sufficient to determine 
the velocity solution for an incompressible material in plane strain, once 
the stress distribution is known. Mohr’s deformation hypothesis cannot 
thus be universally valid. In the present problem the plastic region when 
indentation begins must, as before, have spread sufficiently far to include 
the two slip-lines from the face centre to the free surfaces, since the slip- 
lines are still the characteristics of the velocity equations. The region 
where plastic deformation is taking place is bounded by these two slip- 
lines and not, as Prandtl assumed. by the stress characteristics. Th 
fortunate situation whereby the indentation pressure can be surmised 
without a full solution no longer exists. The part of the plastic region 
uniquely defined by the free surface does not, necessarily, now extend 
over the entire face of the punch. The indentation pressure cannot then 
be found without a calculation of the stress distribution in the plastic and 
elastic regions. For this reason also, Prandtl’s solution for a general 
material cannot be correct. It has been applied primarily in soil mechanics 
where the ubiquitous surface of sliding is substituted for a general law o! 
deformation. 


4. Conclusions 

A main aim of the present paper has been to clarify the approach to 
plane strain problems in plasticity, and to formulate the assumptions 
explicitly. Not only the first writers on the subject, but also the most 
recent writers (e.g. Sokolovsky, 13, and van Iterson, 6) have suggested 
slip-line fields for various problems quite arbitrarily and with no adequate 
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PLASTI( 


discussion of their plausibility. In some cases elementary considerations 
suffice to prove them wrong, but in others only detailed numerical calcula- 
tions can decide. In many problems, it is true, the slip-line field in the 
region of plastically deforming material may be surmised without a full 
solution for the stresses in both elastic and plastic regions. This is particu- 
larly so in steady-state problems,t even more than in non-steady-state 
problems such as the one considered here. In a steady-state problem it 
has only to be assumed that the rigid material outside the region of 
plastically deforming material is capable of supporting the applied stresses 
nits perimeter, corresponding to the proposed slip-line field (which must 
f course be compatible with the velocity boundary conditions). General 
considerations may suffice to show whether the assumption is likely to 
be correct. In a non-steady-state problem it has to be supposed that the 
proposed slip-line field is actually developed during the prior loading: 
this automatically includes the assumption that the still-rigid material 

n sustain the calculated forces on its perimeter at the moment under 
consideration. The difference in degree of the assumptions corresponds, 
n fact, to the extra independent variable characterizing non-steady-state 
problems: that of time, or progress of the deformation. 

It is easy to suggest apparently plausible slip-line fields for a notched 
bar under a bending couple, or for a tensile specimen with a deep rectan- 
sular notch, or even for shallow notches of arbitrary shape, but the 
ssumption that such fields will actually be developed is too risky. This 
s particularly so when, as the work of Southwell and Allen shows in the 
case of a shallow wedge-shaped notch, the actual plastic region develops 
ina direction which could hardly have been surmised in advance. For 
many problems, then, detailed numerical calculations constitute the only 
safle basis. The successful application of relaxation methods to the simpler 


types of boundary-value problems is a valuable beginning. 
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SUMMARY 


In this paper we determine by complex variable methods the distribution of stress 
the neighbourhood of a two-dimensional Griffith crack when the pressure varies 
g the crack. The analysis is extended to deal with problems involving cracks 
aeolotropic materials possessing two directions of elastic symn v. The paper 
neludes with an investigation of the stress distribution in the n bourhood of 
two collinear cracks of equal length, and formulae are found giving the shape of 
cracks and the critical tensile stress normal to the cracks which will produce 
ture. It is shown that the influence of one crack on the other is very small 

rovided that the distance between the cracks exceeds the length of each crack. 
The methods of this paper can also be used to solve the problem of the indentation 
plane boundary of an isotropic or aeolotropic material by a single punch 

Tany shape » Ol by a flat-ended double punch. 


1, THE theory of cracks in a two-dimensional elastic material was first 

veloped by Griffith (1 Recently Sneddon (2) gave an alternative 
treatment of a crack opened by a uniform internal pressure, by using 
Westergaard’s complex stress function (3). In a recent paper (4) Sneddon 
ind Elliot show how to solve the problem when the pressure varies along 
he erack, by an elegant but rather complicated method involving the 
solution of a pair of dual integral equations, and the authors state that 
he complex stress function corresponding to a variable internal pressure 
s apparently unknown. 

In the first part of this paper we give a simpler solution by complex 
variable methods and obtain stress functions corresponding to a variable 
listribution of pressure. In the second part we show how the analysis 
may be extended to deal with problems involving cracks in aeolotropic 
materials possessing two directions of elastic symmetry. In the last part 
€ investigate the stress distribution in the neighbourhood of two collinear 
cracks of equal length. the case of an infinite number of cracks having 
been solved by Westergaard (3). Formulae are found giving the shape 
t the cracks and the critical tensile stress normal to the cracks which 
vill produce rupture. It is shown that the influence of one crack on the 
ther is very small, provided that the distance between the cracks exceeds 
the length of each crack. The theory is illustrated by drawing curves of 


[Quart. Journ. Mech. and Applied Math., Vol. II, Pt. 1 (1949)] 
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maximum shearing stress in the neighbourhood of the cracks when they 
are separated by a distance equal to two-ninths of the length of each 
crack. It is shown that the presence of one crack considerably increases 
the stresses near the other only in the immediate vicinity of the ends of 
the crack. 


2. In this section we consider the effect of a single plane crack in ay 
isotropic material which is in a state of plane strain. The crack may be 
represented by a cut from z c to 2 c, where c is real. The trans- 


formation ae as 
2=ccosl (¢ = é+in) (1 


transforms the region outside the cut into the semi-infinite strip » > 0 


0 


uv 


27 in the ¢-plane. the line » = 0 corresponding to the boundary 
of the cut and y — ~ corresponding to the point at infinity in the z-plane. 
The transformation inverse to (1) is 
git fz—(22—c?)}}/e, 

where the sign of the square root is suitably chosen. 

When an isotropic solid is in a state of plane strain the distribution 
of stress has been represented in compact form with the help of complex 
variable analysis by. Muschelisvili (5) and others. Among the possible 
equivalent representations. for the present paper we write the displace- 
ment and stresses as the real parts of 


2uu = (3—40)9(z)—f’(z)—2Zg'(z), 


2uv = —i(3—40)g(z)—i{ f’(z)+29'(z)}, 

vx f(z) —2g" (z)+29'(z), (3 
yy = f" (2) +29" (2) +29’ (2), 

ay if” (z)—i2g" (2), 


where o is Poisson’s ratio, » is the modulus of rigidity, and f and g are 
chosen to give single-valued stresses and displacements which vanish at 
infinity. The corresponding formulae in generalized plane stress are found 
by replacing o by o (1+). keeping p» unaltered. 
The integrated effect of the stress components (x, ny) round part o! 
a curve 7 = constant which surrounds the crack is given by the change 
in value of 
| (ny+inx) ds [ f' (2) +29'(z)+-9(2)] = — | (yn +1&n) di, (4 
0 
as z traces out the curve. In particular, this relation holds on the crack 
n = 0. 


We shall consider the case when there is no shearing force along the 
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to apply when shearing force is present. 


system. 
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En 0 when z is real, 


We shall 
These conditions can be satisfied by taking 


f'(z) = g(z)—2q'(z). 


fhe stresses and displacements are then given by the real parts of 


Qu 2— 40)g(z)+(z—Z)g'(z). 
Que 1i(1 —o)g(z)+i(z—2Z)g’ (2). 
C2 2q' ( Z)g"(z), 
yy = 29'(z)—(z—Z)g"(z), 
cy i zg" (z). 


e maximum shearing stress 7 may be found from the relation 
T 1 ¢x—yy+ 2ixy|, 


l.on using (6), this becomes 


ar 4y7g"(z).g”(Z). 


m the crack when 0, equation (4) reduces to 


/ 


c | msin€ dC. 


also restrict 
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but the analysis may easily be extended 


our 


ittention to those distributions of pressure which form a self-equilibrating 


(5) 


(6) 


(8) 


(9) 


From equation (9) we find suitable complex stress functions giving single- 


ilued stresses and displacements which vanish at infinity. Corresponding 


1 uniform distribution of pressure 77 Po, we obtain 


q(z) } p(z— (z*—c*)}}, (10) 
hich agrees with Westergaard’s solution used by Sneddon. With 
' p, cos € the corresponding solution is 

) ” ° 1)9 
q(z Pi fe (22 —¢2)))2, (11) 
SC 
nd with 7 Y) P, Cos nc n 2) the solution is 
(2) cp | | (2 ~2 c*) \" J ] (z (2? ayy (12) 
} n 1 | ( | n- 1| c | 
When 47 q, sin n€ (n 2) the stress function is 
" icq L | sre)" 1 {z—(z?—c?)!\"+1 (13) 
ale = 2 2 
4 “ees ( | n+1\ ( | 
The value 77 q,sin ¢ is not considered as the corresponding stress 
ystem is not self-equilibrating. By superposition of the solutions (10)- 


. we obtain the resulting displacement and stresses associated with 
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a general variable distribution of normal pressure along the length of th 
crack. The energy of the crack and the criterion for rupture may }y 
calculated as in Sneddon’s paper (2). 


3. Generalized plane stress problems in aeolotropic materials with two 
directions of elastic symmetry have been considered by Green and others 
We shall use the notation of one of Green’s papers (6), suitably modified fo 
application to problems involving plane strain; and we shall consider onl 
the case when the crack is parallel to one of the directions of symmetry, 
where c is real 


The crack is represented by a cut from z = ¢ to z C, 


Then if z, = w+iA,y, 2% = a+r, y. the displacements and stresses corre- 


sponding to plane strain are given by the real parts of 


u (aA +b) f’(z,) —(aA3+)q’ (22). 
v i(bA, +dAz1) f'(z,)—t(bA,+-ddaz})g' (22), 
LX Az f"(z,)—AZg" (z2). (14 


yy =f" (ez) +9" (22). 
xy iA, f "(2,) 


where f and g are chosen to give single-valued displacements and stresses 





As g' ( 2), 


which vanish at infinity. The constants a, b, d, A,, A, are related to the 


elastic constants s;; by the equations 








9 
Ss Bas SJ< 
@ — °11°337*13° 
> 9 
8 ' 2 
“33 9 S Sac S5< 
dF AS __ Mites 3 
813 803 — 810 Sg $11 533 — $13 . 
b — °13°237 *12%33° (15 
> » ) 
*33 we. 33 See “913 523 T “512 S33 
9 2 a aa 2 
d — °22°33— *23° Ni AR S22 533 98 
°33 


For generalized plane stress the corresponding results may be obtained 
from Green’s paper (6). 


We shall restrict ourselves to the case when there is no shearing force 


along the crack, and the applied forces are self-equilibrating. Thes¢ 
conditions can then be satisfied by taking 
A f(z) +Asg (2) (0). (16 


The displacements and stresses then become the real parts of 


u (ad3+b) f’(z,)-+ (aA3-+b)A, Az f"(29), 
v i(bA, +d Az?) f’ (2) + i(bAg+-dAz2)A, Azf’ (22). 


Ee rE f(z) +A, As f” (22), (li 
yy = f"(a)—A Az f’ (22). 
2 ir, f" (2) Hid, f’ (22). 
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Suitable forms for the function f(z) may be deduced from the corre- 

sponding stress functions for isotropic material. In the case when the 
rmal pressure along the crack is uniform, the appropriate stress function 

3 given by 
. pas ( ‘ 9 ) . 

2_ {> __(22__¢2)hh. (18) 


[he displacements and stresses are then obtained from equation (17). 


In particular, the displacement at points on the crack is given by 


, pa(™ l *2) xr 


A, As 


(19) 


showing that, as in the isotropic case, the shape of the crack is elliptical. 
[he potential energy of the crack may be shown to be 
[Ay +A, 
i jdp*e2| co ’ (20) 
: A, A> 
By using this expression and an argument similar to that of Sneddon (2), 
»can show that the crack will extend when the normal stress exceeds 
the critical value p,, where 
( 47A,A, | (21) 
f ’ - 
Ay \aed(A, T A,)} 
nd T is the surface tension acting round the crack. Our formula (21) 


educes to Sneddon’s formula (2.3.4) when the material becomes isotropic. 


4. In this section we find the distribution of stress in the neighbourhood 
{two equal collinear cracks in an isotropic material, when a uniform 
pressure acts normally along each crack and there is no shearing force. 


[he complex stress function found in § 2 for a single crack under uniform 


ressure could also have been obtained by considering the following 


undary problem for the half-plane y > 0. When z is real, we must 
tisty 
I yy Pp { | 
} () (\a ( 
b) xy 0 for alla 


e stresses and displacements given by (6) satisfy condition (6b); and 
dition (a leads to 
Real part ot 2q'(z) p (|r ; Cc), 
21(1 oC 


g(z)| 0 (jz| >). 


Real part of 
1 | 


Rn 
D 


it these conditions correspond to the hydrodynamical problem of the 


inform motion of a flat plate of width 2c moving normal to its plane in 
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a fluid at rest at infinity, and the complex potential function which 
represents this motion is well known. 

In the same way, the problem of two collinear cracks on the x-axis fron 
z= ktoz= 1 and from z k to z 1, is equivalent to the hydro- 
dynamical problem of the uniform normal motion of two equal collinear 
flat plates through a fluid at rest at infinity. The appropriate complex 
potential function is given in Durand (7); and we deduce that a suitabk 
expression for g’(z) is given by 

9 (2) = $7] = — “ a =1| (k <1), 22 
V{(z?—1)(2?—k*)} 
where A? = E’/K’, and K’, E’ are respectively the complete elliptic inte- 
grals of the first and second kind associated with k’, the modulus comple- 
mentary to k. Our choice of coordinates for the ends of the cracks involves 
no real loss of generality. 
From (22) it follows that a suitable form for g(z) is 
g(z) = 3p 


1 


= ——. dg — 4.2. 93 
(e—1(e—ey 2? 


The substitutions z = de(f.k). t = de(u,k) transform (23) into 
q(z) }p{sn (de l+¢(— E(f)—A2f—de Z}. (24 
c 
where E(¢) | dn? u du, 
0 
and the elliptic functions are referred to modulus k. It is convenient to 
make the further substitution { = ir, when equation (24) gives for points 
x on the crack, - en , ‘ , . 
: q(x) sip{ E(r, k’)—A?7+ iz}, (25 
where x = dn(r, k’). (26 


[t follows from (6) and (25) that displacements at points along the crack 

are given by the relation 
(l—o)p.7 ? ‘ 7 
, SUP { E(r, k’)—A®r}. 2 

be 
The elliptic function in (26) and the incomplete elliptic integral in (27 
have been tabulated (8,9), so that the shape of the cracks can easily be 
obtained with little numerical computation. The energy of both cracks 
is given by the formula 


= 
Ip sm da { E(r, k’)—A?rhk'? sn(z, k’)en(z, k’) dr, 
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which after some calculation becomes 


ar( | a) p*, 


L=' (1—)2—4h’2}, (28) 


isa useful check we may verify that as the distance between the cracks 
tends to zero, the above expression becomes identical with the energy of 
single crack of length 2 units. 

[he theory of the Griffith crack considers the stress distribution when 
ie crack is opened by the application of an average tension p to the 
surface of the body, and the surface of the crack is assumed free from 
tress. This condition may be obtained by superposing on the stress 
inction (23) a second function giving a uniform tension perpendicular 

the cracks. The presence of the two cracks then lowers the potential 


nergy of the body by an amount L. But the cracks have a surface energy 
l 4(.1—k)T, (29) 


ere 7’ is the surface tension of the material. Thus the resultant loss 


f potential energy due to the presence of the cracks is 


The condition (@é/ek)(L—U) = 0 holds when the two cracks are just 
the point of extending inwards. This leads to the result that the 


naterial will rupture if the tension exceeds the critical value p,, where 


K'k'h { 4Tu ) a 
p ; . 3} /}. b (30) 
ID K'k?} J \x(1—co)} 
The corresponding critical value for a single crack of length (1—;) is 
given by 


] { 47 un | 


’ (31) 
v(1—k) a \a(1 o)} 


Pe 
The value of the ratio p,/p, is a measure of the amount by which the 
hody is weakened by the presence of the second crack. The variation of 


this value with the ratio 


distance between cracks 


length of crack 


s shown in Fig. 1, from which it can be seen that if the cracks are 
separated by a distance comparable with the length of each crack, then 
the effect of the second crack is very small. 


[he maximum shearing stress may be shown from equations (8) and 
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(22) to be given by the relation 


r (1A 2A*)y(a?+y?)!a?-+ (y +B)2}a2-+ (y—B)"} 
. | y ( > 


r—1P+y athe tye —kyP ty 


where p 











oO 
a ] 


20  ©»~—— 30 ~ 4-0 


——— f ——_ » 


Fic. 1. The variation of p,/p’, with r. 


In order to give some idea of the distribution of stress in the neighbour 
hood of the cracks, the maximum shearing stress t was calculated fot 
several values of x and y for the particular case k = 0-1, when the cracks 
are separated by a distance equal to two-ninths of the length of each 
crack. The results are shown in Table I, and the variation of 7 with 4 
and y is shown graphically in Fig. 2, where on account of the symmetr' 
of the problem only one-half of the configuration is shown. For purposes 
of comparison corresponding curves for a single crack are shown in Fig. 5 
these being drawn from results given in Sneddon’s paper (2). It can be 
seen that the presence of the second crack does increase appreciably the 
stresses near the ends of the crack, particularly the inner end, but the 
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TABLE [| 


Variation of r p with « and y. 
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Fic. 4. The isochromatic lines in the vicinity of the double crack. 


general character of the curves is preserved. A convenient method of 
showing the variation of 7 consists of plotting the contours of equal 
maximum shearing stress, i.e. constructing the family of curves 7/p = 4, 
where « is a parameter. These curves, which are the isochromatic lines 
of photo-elasticity, are shown in Fig. 4. 

The shape of the crack is shown in Fig. 5 which also shows, for purposes 
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f comparison, the elliptical shape of a single crack of the same length 
inder the same internal pressure. It can be seen that the presence of the 
second crack makes only a slight change in shape. 


The analysis of this section could be extended by the methods of §3 








Fic. 5. The shape of the cracks. 


to obtain the stress distribution near a double crack in an aeolotropic 
iterial with two directions of elastic symmetry. The methods of this 
paper could also be used to find the stress distribution due to the indenta- 


tion of the plane boundary of an isotropic or aeolotropic material by a 
single punch of any shape, or by a flat-ended double punch. 


[ wish to thank Dr. A. E. Green for acquainting me with the problems 
fthis paper, and for many helpful suggestions. I also thank my wife for 


preparing the diagrams. 


REFERENCES 
\. A. U.iFFItH, Phil. Trans. A, 221 (1920), 163. 
I. N. SNeppoN, Proc. Roy. Soc. A, 187 (1946), 229. 
H. M. WESTERGAARD, J. Appl. Mech. 6 (1939), A 49. 
. IN. SNeppon and H. A. Exxior, Quart. of Appl. Math. 4 (1946), 262. 
. MuscHELIsviuti, Z. angew. Math. Mech. 13 (1933), 264. 
6. A. E. Green, Proc. Cambridge Phil. Soc. 41 (1945), 224. 


7. Duranp, Aerodynamic Theory (J. Springer, 1935), vol. ii, p. 83. 


we wD Ne 


8. L. M. Mitnr-Tuomson, Di elliptischen Funktionen von Jacobi (J. Springer, 
193] 

9. A.M. LeGENDRE, Tables of Complete and Incomplete Elliptic Integrals (Cambridge, 
1934 











THE FORMATION OF CLOSED WAKES 
IN FLUID MOTIONS 


By D. N. pe G. ALLEN (Jmperial College, London) 
teceived 7 January 1948] 


SUMMARY 
Southwell and Vaisey (1) (1946) have applied relaxation methods to obtain 

variety of solutions to problems of two-dimensional (steady) motion in an inviscid 
fluid. These are concerned with the determination of ‘free’ streamlines, and includ 
one type of solution which is believed to be new, in that examples have been found 
of a flow which exhibits a junction of two ‘free’ streamlines. At such a junction th 
wake boundary must be cuspidal, and a remark of Sir Geoffrey Taylor about tl 
existence of such cusps has stimulated this investigation. It was undertaken wit] 
the aim of finding a solution of cuspidal type by ‘orthodox’ mathematical analysis 
i.e. by an extension of the classical method (using conformal transformation) whic! 
was employed by Helmholtz, Kirchhoff, and Rayleigh. 


1. THE method of determination of ‘free’ stream boundaries, in the two- 
dimensional steady motion of an inviscid fluid, by the use of conformal 
transformation is described by Lamb (2) (1932), who gives a number of 
illustrative examples. We write 

z= a+, w= d+, (1 
where ¢ is the velocity-potential and % the stream-function of the flow. 
Then ia er . 

dw/dz = ¢d/ex—i.éd/ey = —u+iv, 

uw. v being the components of velocity at a point in the axial x-, y-directions. 
Let q denote the resultant velocity at a point and @ its inclination to the 
axis of x. Then 


u q cos 8, v qsin 6, 

dz l a 
and we have — - — = —¢! (, say; 

dw u—w q 7 
and log € = log(1/q)+70. (2 
Now log ¢ and w are both functions of the complex variable z. If then the 


complete boundary of the flow can be specified both in the planes of w and 
of log ¢, and if a relation can be found whereby the regions so bounded in 
these planes can be conformally transformed into each other, we shall have 
obtained an equation of the form 


F (log f, w) 0. (0 
On a free streamline the real part of log ¢, i.e. log(1/q), has a known constant 
value, since the pressure is constant and therefore also the velocity, q = 4: 


[Quart. Journ. Mech. and Applied Math., Vol. II, Pt. 1 (1949)] 
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savy. Also on a free streamline the imaginary part of w, i.e. i, has a known 
onstant value, and the velocity-potential satisfies the relation 


0d /es a, (4) 


vhere s is arc-leneth measured along the free streamline. Integration of 


4) along a free stream boundary, with a suitably chosen origin for s. gives 
db —iie. (5) 


Thus (3) reduces along a free streamline to a relation between @ and s: it 
s the intrinsic equation to that boundary, which is therefore at once 
lefined. 

The appropriate form of equation (3) is found by introducing an inter- 
nediate function ¢ between log ¢ and w: the corresponding regions in the 
planes of log € and w are in turn transformed into a half-plane in the plane 


ft. bv relations of the form 


f (log t. $i 0). q(w. t) 0. (6) 


nd ¢ can then be eliminated between the equations (6) to obtain (3). Or, 


sufficiently, we may regard (6) as giving the transformation between log ¢ 


nd win parametric form 


2. It is convenient to take units of leneth and time so that the value 


fthe velocity along a free boundary ‘open’ to the atmosphere is 1. Along 


the boundary of a ‘closed’ wake. of the type described by Southwell and 
Vaisey (1), the velocity g has some other constant value, a. less than 1, 
idg = 0 at a stagnation point. In examples previously treated by this 
method, by Helmholtz (3) (1868). Kirchhoff (4) (1869), and Rayleigh (5) 
1876), all free streamlines are ‘open’ to the atmosphere, and so a single 
mstant value can be given to g along all parts of the stream boundary 
hich are free. In the example treated here the boundary of the stream 
mprises both an ‘open’ and a ‘closed’ free streamline, so that different 
nstant values must be given to qg along different parts of the stream 


boundary which are free 


3. Complete specification of the fluid boundary, both in the planes of 
gC and of w, is possible when every part of the bounding streamlines is 
either free or fixed and straight: a fixed but curved boundary introduces 
mplications. In the plane of log ¢ a fixed and straight boundary (since 
‘is then constant) becomes a line parallel to the real axis, whilst a free 
undary (since log(1/q) is then constant) becomes a line parallel to the 
maginary axis. In the plane of w. the stream-function % is constant along 
each bounding streamline. which therefore becomes a line parallel to the 


eal axis. 
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4. The example to which we apply this method of solution is illustrated 
in Fig. 1. The axes Ox and Oy represent two fixed straight boundaries 
fluid flows in a stream of finite breadth from the direction J far upstream 


and in contact with the fixed boundary Oy. The direction of flow is turned 


y l 





>_> 


Fic. 1. Plane of z x+y. 


through a right angle, and far downstream the fluid flows in the direction / 
in contact with the fixed boundary Ov. A pocket of air (or of ‘dead’ fluid 
is entrapped between the fluid stream and the fixed boundaries at the 
angle O so that one of the two bounding streamlines, yA Bx, is in three 
parts—yA and Bz being fixed and straight, and AB being a free boundary 
of a closed wake; the other bounding streamline, JJ’, is a free boundary 
open to the atmosphere. Far upstream towards J and far downstream 
towards I’ the uniform velocity of flow is given by g = 1 which is also 
the velocity at all points of 77’. On the closed boundary A B, the constant 
velocity is given by g = a < 1. We take d to denote the uniform width 
of the flow both far up- and downstream, so that if the stream-function / 
is given the value 0 along JJ’, then along yA Bx we have & = d. 
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5. The boundary of the region of the flow in the plane of log ¢ is shown 


» Fig. 2: the fixed and straight boundaries JA and BI’ become parts of 
the lines @ 7/2 and 6 0. while the free boundaries AB and JI’ 
ecome parts Ol the lines log(1 q) log(1l/a) and log(1 q) 0. The first 


|0 

















B log (1/4) 
| A 
IG. 2 flog ¢ — log(I @ 
6 ,O a 2 
| 


3. Plane of ¢. 


onformal transformation to be found is that which transforms the rect- 
ingular region JA BI’ in the plane of log¢ into the half-plane above the 
real axis in the plane of the intermediate variable ¢ (Fig. 3). The rectangu- 
lar boundary in the plane of log ¢ is transformed into the real axis in the 
plane of ¢. The values ¢ l and t 1 may be assigned to the points 
l’ and I respectively, and then f¢ —(1/k) at B and A, where k has a 
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value (< 1) dependent on that of a. The transformation of the rectangle 
into the half-plane is given by the Schwarz-Christoffel formula in the form 


t 
P dt 


J fa—#8)0—Re)} 


0 


log ¢ x +B. (7) 


The constants « and f are determined by the conditions that, at J, 


log ¢ —iz/2 when t 1, and, at J’, log = 0 when t 1; therefore, 
from (7), 
var /2 rk +B | 
and 0 vA +B J’ ") 
: 
where kK ae (9 


/{(1—t?)(1—K42)) 


. 


0 
4K is the real period of Jacobi's elliptic function sn(w, &) and the imaginary 
period of this function is 27K’, where 
‘ dt 


Kk’ . 5° (10) 
J{(P®—1)1—k*#?)} 


Solving equations (8) for a and £ and substituting in equation (7) we find 
that t 
4Ki, , ‘ 4 dt 
log S kK ss ( 9 ~ 7,942\) 
Vid t*)(1 - k?#?)} 
0 
Pp 


7 


or t sn 


log £ K). (11) 
\ 

6. The second conformal transformation to be found is that which 
transforms the plane of ¢ into the plane of w. The region of flow in the 
plane of w becomes an infinite strip between the lines 0 and J=d 
and is shown in Fig. 4. The velocity-potential is conveniently given the 
value 0 on the line of symmetry CD and hence, using equation (5), we see 
that the origins from which the arc-distances s are measured along the 
two free streamlines are the points C and D (Fig. 1). The infinite strip 
is a rectangle with two vertices at J and two vertices at J’ so that the 
transformation between w and t is again given by the Schwarz-( ‘hristoffel 
formula: 
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The constants y and 6 are determined by the conditions that at C, w = id 
©, and at D, w = 01 


vhen t 0: therefore, fr 
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Fic. 4. Piane of w dh +- ius. 


7. The complete transformation between log and w is given by the 
equations (11) and (13) taken in conjunction. To find the equation of the 
closed free boundary AB we substitute in these equations the conditions 


log(1/a)+76 | 
as | id j 


btained from equation (5) and § 4. Thus, along AB. 


_ 1Ke 
K 4 te ma th A) 
a 7 


1—t\ ad 
al; re 


3) that ¢ is real along AB. and, therefore, in (15) it 


K log + 


vr a 


t sn( 


~-K+i2K'’—4K06/z) 
sn(K+iK’+4K0/7) 


1 dn(4K6/z) 
k on(4K6/7) 





















(14) 
(15) 
(16) 
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(18) 
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Equations (16) and (18) form the intrinsic equation to the closed boundary 


AB in parametric form; elimination of t and differentiation yields 


ds ; Skdk n( ==") 


-———— (19) 
dé 7a 7 


\ 


8. In a similar way the equations which hold on the open boundary J]’ 
may be found. There we have 


og f= 10 | 
(20) 
and w s | 
in place of (14), so that 
t sn(K +4K 6/7) (21 


d 1|—+t 
and s log ) (22) 
T \1--¢ 


replace (18) and (16) and 

ds a SdK ] 

de a sn(4K 6/7) 
replaces (19). 


9. The boundaries of the flow (in the plane of z) have been computed 
for two particular values of a. The first case is that in which a = 0 and 
the closed boundary AB reduces to a stagnation point at O (Fig. 1). We 
have then that 

h 0, K = x/2, RK’ == 00, (24 


and the equation (23) to the open boundary JJ’ becomes 


ds dé (4d /zr)cosee 26; 
therefore dx dé (2d /z)eosec 8 
and dy dé (2d 7 )Sec A. 


and, on integration, 
x = d—(2d/z)logtan(—6@/2) | ne 
- (25) 
y d—(2d/z)logtan(6/2+-7/4) | 
Fig. 5 shows the shape of the boundaries in this flow of a fluid stream 
round a right-angled corner. The result presumably is not new, but no 


previous reference has been found. 


10. The second computed example is one in which a has a typical 


non-zero value. To lighten the arithmetic k is assigned the value 1/\- 
so that 


a 0-45593 | 


kK K’ 1-85407 | 


(26) 
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dar \long AB we have, from (19), 
da SKdk 4k0 
sn| J cos 
dé Ta 7 
Q , (27) 
dy SKdk 4K6\ . 
nd ; sn }sin 6 
an dé TA" 7 
Vil . 


nd along JJ’, from (23) 


da Skd cos ¢ 
yA dé 7* sn(4K6/z) | 
(28) 
dy Skd sin @ 
dé 72 sn(4K6/z) ' 


<«—_~< 











Wi 
>> 
> xX 
Fic. 5. 
On the free boundaries values of « and y are found by numerical integra- : 
n of (27) and (28). In Fig. 1 the flow is shown reflected in the axes 
Oy; in this way it becomes a flow. with no fixed boundaries, of two 
ial and opposite streams each of which, on meeting, bifurcate into 
its which pair with one another, departing in opposite directions at 
| | 
° ght angles to the original direction, and including a four-cusped wake 
t the junction. 
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RECURRENCE RELATIONS FOR CROSS-PRODUCTS 
OF BESSEL FUNCTIONS 
By E.T. GOODWIN (Mathematics Division, National Physical Laboratory 
[Received 10 February 1948; revised 21 July 1948] 

THE solutions of some electromagnetic problems contain cross-products of 
Bessel functions, frequently of complex argument. It is most convenient 
to compute these from the recurrence relations which they can be shown 
to satisfy and which are detailed in this note. The notation used is that 
of Watson’s standard treatise. 

if Pn =I, (®)¥,(y)—I,A(y)¥ (x) 


a 


In = I, (@)Y 2 (Y)—Sr(y)Yn(*) 


: So (la 
rn a J, (x)d nl¥) —JS,(y)¥ n(X) 
Sn = J (x) ¥,(y)— J) (y)Y 7, (x) J 
where dashes denote derivatives; then 
2n 2n 
Pn 1 Pn 1 —— 4 “Tn 
a y 
n n+] 
Gn+1 ln “Pn — — Pni+t 
A Y 
n n ] (1 h 
ln 1 q, . Pn a : Pn+i 
Y 
1 n= 
>) 2p 1 2P yn 1 Fe 
xy 
1 
Pn Sn Gn ln — 
T°xy 


The last relation is analogous to the Wronskian and affords a powerful 


check on the values of p, g, 7, s as they are calculated. 


Similarly, if fy, = J,lw¥,(y)—J, (WY, (x) ) 


In = J,(")¥n aay) —Fnaaly)¥, (2) +> (2a 
h,, J, s(@)¥,,(Y)—S, (YY naa (%) ] 
then ° 
; tn- 2n 2n } 
Tn 1 tn 1 I Ba Int In+- h,, 
vy = y 
2in+1), 
Yn 1 h,, — fn 1 
y (2b 
2(n+1) 
hiss rga > ? ; =Sn+t 
z | 
var { 
Tai 11 Jn h,, — > | 
TXY J 
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The first of these relations may be expressed in either of the simpler forms 





; 2n 2n 2n 2n 
j 1 » * 1 g In 1 h a 2 h n—1° 
x y y x 
ratory Similar relations for the functions J, K are as follows: 
ucts of Pn A K,,(y) [,(y)K,,(«) 
enlent Y, I \4 )K(Y) I’, (y)K,, (2) (3a) 
show ry, = T,(x)K,(y)—L,(y) Ky (2) 
is that Sy F (a Ky) I’, (y) Ki (x) J 
then 
2n 2n } 
P ‘= Py l q, “T ln 
1 i] 
nt U7 ] 
Gn+1 y Pn— Pri 
z i] 
n n+] a (35) 
Vy 1 1 Pn ail r Py 1 | 
n- 
S ly 1 — ) 
o/ 1 2P 1 xy" n 
. | 
P Sn q | 
XY J 


In this case the cross-products g, h analogous to (2a) which are of interest 
ntain positive signs. Thus, if 
ps I (x)K,(y)—L,(y)K,(&) 
In = L,(%)Knily)+haQbn(@) +: (4a) 
h,, Dn a(@)K, (YW) +L, (y) Ky a1(*) 


ve 
en 
in? 2n 2n | 
J n+i—J 1 = Ju a er h n 
XY a y | 
f 2(n+1), 
g ] / | * 1 
y 
5) . (45) 
2(n+1), 
h ’ ] — tz +] | 
Hr | 
l | 
- r. 1 g,f 7 
9 , ry J 


For the spherical Bessel functions, if j,,, y, are defined by 


ne (= ) J, 43(2), Y,(x) = Js) +4(7), 
f \Bz 24 




















74 CROSS-PRODUCTS OF BESSEL FUNCTIONS 


and if > 
Pr = Inl®)YnY—In(WY nl | 
In = Jnl YrlY—In(WYnl@) | ; 
” ° ae (Da) 
rn Inl®@)YnlY) Ie Y yi, (2) 
8, = jrle)y(y)—Frlwynle) J 
then 
l l ] 
Pnii— Pn-1 (=n 1) q, . rn r -Pn 
2 y xy 
n n 2 
( — ) 
In 1 y , P; y Py i 
n n-+2 
Vat T Tn Pn- = — Py 1 Pe (5b) 
y Sn 
n+] n n(n-+-1) 
S —_—~— ) - ) —_ ) 
m In!" gpa Pet ay -’ 
l 
Pn Sy dn rn > 9 
wry 


cross-products arising naturally are not 
those of j,,(x), y,,(x) but those of 27, (: 
this case, if 


In scattering problems the 


v), xy,(x) and their derivatives. In 


Pn [a7,,(2) | lyy,(y)| 
Ty, [ aj, (a) | yy, (Yy) |’ 


[Wn(Y¥)| [ey (x) 
LWn(y)|'[xy,, (x 


(x)| 
(x)] 
[winy)| [ry,(x)]’ 
(a) |’ 


° ‘ (> (6a) 
rn [ aj, (x) | lyy,(y)| 
8), [27),.(%)|'Lyy,(y) |’ ly), (2 /)\’ xy, (a } 
then 
7 ae 
pak eae ion yin a y - rye” 
l l 
Ln 1 ' (” I ) Ln aie y Pn 1 
I I . (65) 
" 1 Gs (n 1)( Pn ~Pn+1 
y a 
n n--] n(n--1) 
Ss -_ ) ) ) 
? an j P, 1 In l Pn 1 ry / 
Pn Sn Tn Vn | 


Finally it may be noted that similar expressions can no doubt be 
produced for the cross-products of any functions which themselves satisfy 
similar recurrence relations. 

The work described above has been carried out as part of the research 
programme of the National Physical Laboratory, and this paper is pub- 
lished by permission of the Director of the Laboratory. 
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SUPERSONIC FLOW PAST SLENDER POINTED BODIES+ 


By G. N. WARD 
Department of Mathematics. The Unive rsity, Manchester) 
Received 6 January 1948] 
SUMMARY 

By assuming irrotational inviscid flow of a perfect gas moving past slender pointed 
ies at supersonic speeds, the limiting values of the aerodynamic force coefficients 

s the thickn tends to zero are determined for bodies of general cross-section. 
The expression for the drag coefficient at zero incidence is the same as that for a 
ly of revolution when expressed in terms of the cross-sectional area at any point. 


lateral force coefficients and the extra drag coefficient at incidence depend only 


the shape of the base section 

[he general results are applied to the problem of a body of revolution carrying 

gs of smal pect ratio, and the interference between the body and the wings 

letermined ; this leads to a particularly simple expression for the lift in terms 
lift of the wing alone and that of the body alone which, it is suggested, may 

ly approximately to any wing system on a body of revolution. 


1. Introduction 
[HE isentropic flow of a perfect inviscid compressible fluid past slender 
winted bodies of general cross-section at supersonic speeds will be treated 
an approximate method based on the well-known linearized equation 
motion. The body, either pointed at both ends or with a flat base, 
vill be taken to be of unit length (no loss of generality is incurred, since if 
e wish to treat a body of length 7, we have only to scale up all the 
lengths involved by a factor /) and of maximum thickness ¢. Since we are 
onsidering only slender bodies, ¢ must be small compared with unity; 
he angle which any tangent plane to the body boundary makes with the 
indisturbed stream direction must be small and O(t), and the rate of 
hange of this angle along the direction of the body must also be small 
nd O(t). One further condition on the shape of the body is required in 
general: this concerns the radius of curvature of any section of the 
oundary of the body by a plane perpendicular to the stream direction. 
Ifdis the maximum diameter of the section, then the curvature must be 
Y(1/d) for all points where the section is convex outwards: there is no 
restriction at points where the section is concave outwards. This last 
ondition is not always necessary for bodies at zero incidence. but is 
ways required if we wish to calculate the flow at incidence within a 
nown approximation. These conditions ensure that the extra velocity 
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due to the disturbance caused in the uniform stream by the presence of the 
body is everywhere small compared with the velocity of the main stream, 

Some interesting body shapes are excluded by the above curvature 
condition, and it may be relaxed under certain conditions, but in these 
cases the degree of approximation of the solutions becomes difficult to 
estimate. Examples of such bodies are treated near the end of the paper. 


2. The solution of the equation of motion 

Consider the isentropic flow of a perfect inviscid compressible fluid, 
We shall take Ox, Oy, Os to be rectangular Cartesian coordinate axes 
forming a right-handed system, the undisturbed stream at infinity to be 
parallel to the s-axis flowing in the direction of increasing s with velocity 
U,, and c, to be the velocity of sound in the undisturbed stream. It may 
be shown that the general velocity components parallel to the coordinate 
axes are given by 

u,. y 2 u pit U, o(1+ 3), (1) 
Cx cy 

where the velocity potential ¢ satisfies the well-known equation of com- 
pressible flow 


c T24 — (1+) sat (; \S x2 +(4) <e49(1 *) oh ¢ oC? 
\ C 


1 > 9 ‘ ‘ 9 ai a 
™ os os* xc} eax cy cy” c CX CXCS 
/ > ¢ 9 
-2(14.§ ve pale AON ps dpe de 
Cy Cycs Ox CY Cxcy 


¢ being the local caged of sound given in 


anata (OF) +(9)] 


and V? being the Laplacian operator. 

As a first approximation we assume that the disturbance velocity UV¢ 
is everywhere small compared with l’, so that the squares and higher 
powers of the velocity components may be neglected. Equation (2) then 
becomes 2 

V2d wee? (4) 


9? 
os” 


M Uc, being the stream Mach number. This is the linearized equation 
for the velocity potential. We shall be considering supersonic flow, s0 
that M 1, and we shall write B? — M?—1. 

By transforming (4) to cylindrical polar coordinates r, @, s such that 
x = rcos#, y = rsin#@, we obtain 
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Consider now a slender pointed body placed in the stream with its 
nointed nose at the origin. The disturbance caused by the body will be 
confined to the region behind the Mach cone s Br from the nose and 
so the disturbance velocities and the potential vanish for s < 0. Thus 


the Heaviside operational form of (5) for this problem is 


o*dh - | Ch ] Ch 
} = ‘a 


|— —. — Bp, (6) 
OY’ ¢ OT FoF 
We require a solution of (6) representing waves travelling outwards and 
downstream from the body only: hence the general solution of (6) for 
flow past finite bodies in s > 0 is, assuming that the series obtained by 

integrating converges outside the body under consideration, 
¢ = A, K,(Bpr)+ > A,, K,,( Bpr)cos(né+-B,,), (7) 

n=1 

vhere the A, and 8, are arbitrary functions of p which must be chosen 
to satisfy the boundary conditions on the body; K,, are modified Bessel 


functions of the second kind, which for large r have the asymptotic value 
7 1 
kK Bpr) ~ (- . Pea, (8) 
2 Bpr} 


showing that (7) represents waves travelling outwards from the axis, and 
long the surfaces s— Br = constant at infinity. 
Near the body, where 7 is a small quantity and O(t), the dominant 


terms of (7) are included in 


| 2 n 
dp log i Bpr ; \A, } > (n (5 A, r-” cos(né +B); (9) 
n=] \ P) 
y being Euler’s constant here. The error is a factor 1+ O(t). 

If our slender body lies just inside a cylinder r = R = constant, where 
R is O(t), then it will become apparent later that the series obtained by 
nterpreting (7) and (9) certainly converge for r R, and ¢, may be 

ntinued suitably outside the body for r R if they do not converge 

1 this recion 

By writing z vty e'’ (9) may be expressed as the real part of 

w = dytiby = aglogz+b,+ Da,2-, (10) 
n=1 
where a,, bo. a, ete., are functions of s given by 





(log t Bbp—y)A,. a 1eiBn(n 1) . 
oer 7 7 I 
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and us, is a real function of s, 7, 6 which is not a true stream function for 
the motion, but is denoted by ys, because it has the properties of a stream 
function in a two-dimensional incompressible flow. Thus dy is a harmonic 
function of r, 6 approximately near the body, and the results of two- 
dimensional incompressible flow may be used to determine w—b, from the 
boundary conditions on the body, and hence 6b, from a», by means of 
the second equation of (11). The expansion for w given in (10) certainly 
converges for |z! > R since there can be no singularities outside the body 
but may have to be continued analytically for z| < R. Usually, however 
it is the analytic continuation of w—b, which is obtained from the boundary 
conditions, and dp, a,, etc., have to be obtained by expanding the function 
in the form (10). 

In order to determine the accuracy of (10) or (9) as a solution of the 
equation of motion, it is advantageous to transform (2) and (3) to the 
independent variables z, 7, s, where z = x+ iy, 7 = x—iy. The equation 
of motion becomes 


nt 4 é 24 i =a) = (1 re. “ys > ry ( P 4 =(: d od 14 ( \: ‘ ao 
2 os* Cs cs* CSI OS : C2 Gz” 






— 
t 
~ 
i 
~ 
ay 


A 

z 
~ 

z 
~ 


y Ca Cd\?2 Od Od 
where c ei om an. Mt 2 +(4) +448 | ; (13 


Oz OZ 


It is now necessary to assume some of the results which will be proved 
later: when we come to determine the coefficients in the expansion (10 
it will be found that a, and b, are O(t?) and that the a,, (n 1, 2, 3,...) are 
O(t"*2). By using these results, it will be seen from (10) that near the 
body where the series converges 


hi 24 
C@d C"-O 5 
~, — are Olt? logi),. 


> 


Cs Cs S Cc 


e2 
ete. are O(t), and - “bo ete. are O(1). 


a1 
C dy C Do 
C2 Cs Cz* 


Thus near the body we may write (12) in the form 


Ps et 400) ( : f) (ee) a ( e) ae Op oh P| 


> = 9 9 F - =( 
CZCz2 cae Os 2 


Cc 


Gaz OZ C2 C204 


Gz Oz cz oz 
_QO(tlog*t). (14 
If we now assume that 


d db, dy, S(w } v)+4,, (15) 


where ¢, and its derivatives are of sufficiently smaller order than ¢o and 
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ts corresponding derivatives so that ¢, may be neglected in the right-hand 
ide of (14). then substituting in (14) we obtain 





Vl { | +s(se) os wpe i, (16) 
les\ez ¢ 2\ez}] oz? ' 2\az] ez? | 
ri2.2.8) say. 
From this we have, as a particular integral, by solving in series, 
Lf | F035) ded24 > © [ [ [ [ fe. z.9) (dz)2(dz)2-+... (17) 
» e 16 . s* e * . * 


Q 
Vi2\ 3 * lew\2 low Cw\? ,_ 
(wir) is | dz +- | — an) 4..... (18) 
{ s 2 ; C 202. OZ } ; 
m which it will be seen that 
; is O(t4 log*t 
24 te are A] O02): 
és, 07h, /és, ete., are Olt log*t): (19) 
Dy C2 0", CZ2C8 etc., are O(# log t), and 
Dy oz". ete are O(t? log t), 
the terms omitted in (18) being of lower order than those shown. 


Now our conditions on the shape of the body ensure that ¢, and its 


lerivatives are of the same order of magnitude on the body as they are 


near it, so that the order of magnitude of ¢, and its derivatives on and 


ear the body are those given above. Hence ¢, gives the velocities on 


ind near the body within a factor 1+ O(¢? log t). 


3. The boundary condition on the body 


Let the contour C be the cross-section of the boundary of the body in 


plane s = const., and let C’ represent the projection on the same plane 
the cross-section at s+-ds (see Fig. 1). If v is the outward normal and 
ris the tangent to C' at any point, and dv is the distance between C' and C” 
easured along the normal, then the boundary condition of zero velocity 
rmal to the body is 
" dy ch 
14 ) (20) 
ds Cs 
\ 7 y > 
\ow dvds is OU by definition. so that é¢/cv must be O(t) on the body, 


nd if the curvature of C satisfies the condition of $1, then @¢/é7 will also 


Q(t) on the body. Thus |éu is O(t) on the body and near to it, and 
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SO Mp, a,, etc., in the expression (10) will be O(t?), O(#*), ete. ed/és wil 
then be O(t? log t) and (20) becomes 





F i a © 4 Ont. (21 
Cv cv ov ds ° 7” ° 

ma Vv 

7 
y> | 
C 
fod 
Fic. 1. 


Thus, by using the results (19), it will be seen that we can determine 


od, from 


Ody dy 
— ? (22 
CV ds 
and we shall have 
Cb od 
= —/1 O(#logt), 
CX CX 
Ch Ch P 
— 91 O(log t), 23 
cy cy 
Ch od . 
= —24 O(t*log*t). 
Cs Cs 
4. Determination of the coefficients a, and b, 
‘ ; : r Cho > ° ‘ ; . 
Consider the integral dr. By using Gauss’s theorem. since ¢ 
CV 
( 
satisfies 
cd od 
Po) “0 0. 
Gn? ey 
the contour C' may be replaced by a circle of small radius 7, R. centre 
the origin of 2 and y, and the expression (10) may be used for do, sine 
it converges for r R. Thus we get 
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Also, from the boundary condition (22), if the area enclosed by C is 


enoted by S(S) . J . 
. dy = ~ 
“7° 4. dr S'(s). (25) 
CV Z ds 
tf ( 
Then, comparing (24) and (25), we have the result 
ws ¢ 
a S'(s). (26) 
) 


From (11), by using the product theorem to interpret the operational 


form. we have for b 


bh - S'(s)log log(s a)S"(a) do 


We sometimes require the derivative of b) with respect to s when 
leulating the pressure: this is 


3 4 m 
}! S”(s)log : S”(O)log s log(s—a)S (c)do\. (28) 


5. The drag force on the body 


We shall calculate the drag force from the rate of change of momentum 


the s-direction through a cylinder r = 7,, where 7, is small but greater 
in R, with plane ends at s = 0 and s 1. Let S, denote the end at 
0, S, denote the curved surface r = 7,, and S, denote the end at s = 1 


itside the body (see Fig. 2). Then the drag force D, measured parallel 


the main stream direction, is given by 


F P | db d 
D= | (py+p, U2)d8,— | pU%(14 —*\: dS, 
Cs] cl 


. 


| $\2 
/P pl (1+) | dS,—Pe S(1), (29) 
lars 


here Pp is the base pressure if S(1) + 0. 


Conservation of mass flow through the cylinder gives 


| p:U dS,— | pU ds, pl (1+! ‘| dS, = 0. (30) 
‘ Cl 7 . cs ; 
Therefore, by mult iplving (30) by U and subtracting from (29), we have 


r : . 7 ,CO CO i 9 b db Y Y 
2D Pp, as, ol e = aS | 'P | pl - (1. c )}aSs PR S(1). 


Ss Cs 


. 
S S S 


(31) 
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From Bernoulli’s equation, by using the results (23), we get 


P =< % 7 9 obo _ “ta / Cho “4 O(t loz), (39 
p, l Cs cr m\ OG; - 
and from the condition of isentropy we get 
p/ Py 1+-O(# log t). (33 
ee ee 
fe % 92 Is. \ 
bee ' l 
. 7 
. # / 
a nnn ena DanS a, 
Fic. 2 


By substituting these last two results in (31), we find 


D ~~ Po Po gg hese ee do “| dS, 
Lo, U? 7 or Os ” . \\ er a } 


ry 


ate S(1 )+ O(t® log?t). (34 
bP 


The quantity (p,—pz)/3p,; U? is usually called the base pressure coefi- 


cient, and is denoted by Cz. 


The integral over S, can be transformed into integrals taken round the 
: 3 2 


to 


7 


boundaries of S,; by using Green’s theorem and we can write (34) as 
1 27 
D r £ [edo ed r c 
oe sae nl | | (<9 7 r, dds +. | $y 2 r,d0 
tp, l ' \ a (or O8 J, r F. Or |r r1,8=1 
0 0 


“| ] $y Pod | +0, S(1)+O(t log). (35) 


‘Ss 1 


We may now substitute the expansion for 4, in the first two integrals 


r) 1 


O(t® log?t). (36) 


of (35) and we get 


1 
) e 
l = = 47 | ayb,ds—2zx/(a,b,) at } by Pdr 
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By putting in the values of a, and &, found in § 4, (36) becomes 


1 1 1 
) 23 l S’(1) [¢ l ” 
l . log we (sys (c)dods — 4) | log -S (c)do— 
U? 2a os s oC 2a a ] oO 
0 0 0 
: € 4 . 9 _ 
| do mad +C), 5 S(1)+O(tlog*t). (37) 
CV 
C /s=1 


This appears to be the simplest form into which the expression for the 
drag force can be put for the completely general case. Two somewhat 
general special cases arise for which the result assumes a simpler form. 


These are (i), when S(1) and the body is pointed at both ends, and 


when the body is of general cylindrical form near s = 1 so that 
(1) = 0, and the generators of the cylinder are parallel to the main 
stream (so that we may think of the body as at zero incidence) which 
makes Chp_/ 1 dv/ds 0. For both of these cases, omitting the base 


lrag term tor case (ii), 


: @ 


aan | J i ; ' 
log S”(s)S" (a) dads -+- O(t8 log?t). (38) 
27T s o 


This is the same result as that obtained by von Karman (1) and Lighthill 
2) for the special case of a body of revolution at zero incidence and 
having S'(1) \). 

The extra terms in (37) occur when S’(1) + 0 or when the body is at 
incidence: the extra drag due to incidence comes entirely from the third 


ntegral in (37) 


6. The lateral forces on the body 
The lateral forees X, Y on the body in the directions of positive x and y 


respectively may be obtained by considering the rates of change of 


momentum in the w and y directions through the cylinder r = 7, as 
fore. We have then 
uy r > OF Ch " F cal c Ch . 
X (pl <— pcos) d8, | pl (1 li ‘) dS,, 
3 cr ca : J CS] Cx 
S Sa (39) 
rd -_ dh Cd. " . i ms c ch ’ 
) (pl 2° psin @)dS, | pl (1+ g . dS... 
J cr cy = - Cs Cy 
So S3 


It will be found to be convenient to treat the complex combination of 
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X and Y, F = X+iY. On substituting the approximate potential and 
retaining only the terms O(t*), (39) becomes 


Ree [ byte 8S _ (2 Po oer OO) th ag [ Of ag 4 
4p, U? J \ or Oz \ Os 02 OZ | i J o& ‘ 
Ss Ss 
| O(t* log*t). (40) 
Now 
Qn ; 
i‘... ' Cho .¢ Po y oo i ' ; r) 
C os dS, = (- x ie Cy dSs ion \ | $y dz ar ($0), rn 11 ao) 
S; Ss : rs 0 )s=1 
, 1 27 
a " f (0d, 0 
= {7 | d,dz -L. | r, e’ déds, (41) 
Cs 
C s=1 0 0 —_ 
by applying Stokes’s theorem for two dimensions; hence (40) becomes 
A 1 27 
F iF * £ 1. 0b, 0 OW OW 
GC —2a | bydz) —_ | | ae en Pe oe r, d0ds +- 
4p, U? J J J\ Of & Gz O),<,, 
C e=3 0 0 


| O(t® log@*t). (42) 
On substituting the expression (10) for w and do, or otherwise by 
remembering that F must be independent of r, and letting r, > ~, we 
find that the double integral vanishes, and we are left with 
F vm | 5 lec? 9 
7a = | 2 | b,dz + O(t° log*t). (43) 
~ c s=1 
This expression for F can be put into a somewhat more convenient form 
as follows. Since ¢, = w—idy. (43) may be written 


* = = | 20 w dz—2 | bode) + O(t* log?t). (44) 
2P1 \ oy FS pen 
Now w is an analytic function of z in a cut z-plane, and C encloses all 
the singularities of w, none of which lies on C. Thus the expansion (10) 
converges on some are of C' and we can choose another contour C, on 
which the expansion (10) is valid and which has at least one common 
point z, with C. We make the cut in the z-plane pass through 2) and 
start and end our contours (, C, at 2 (on opposite sides of the cut). 
Then we have 


m -719'(] , 
w dz = | (“ i | by 24 Sh. 


Rats. .t One 45) 
TIZ9)-+ 270. 
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i F F Cy ; 
Also we have ub, dz [hy z | z——dr. (46) 
OT 
( Cc 


Now since ¢, and ys, are conjugate functions 


ly, Ob, dv ” 
T ol ds 
; * dv a 
that ty |e | = dr S'(s): (48) 
| ds 
ence (46) becomes 
- - | gee 
| us, dz S'(s)zo | 2 dr. (49) 
( 


The last integral has a simple geometrical interpretation: it is the rate 
f change (with respect to s) of the moment of area of C. If we let the 
ordinates of the centre of area of C be given by x 


> Yq 80 that 


z (s) Ly t+Yy. 


then 2 dr 1 S( s)z,(8) | S’(s)z,(8) f S(s)z,(s). (50) 


By substituting all these results in (44) we obtain finally 
F 
l l 


1 


t77(a,),—, +28’ (1)z,(1)+28(1)2,(1)+ O(# log?t). (51) 


The moments of the lateral forces about the 2 and y axes may be 
ilculated from this expression. If we write 
F'(s) 
oP | Z 


on 
to 


t7(a1),+ 2S"(s)z,(8) + 2S(s)z9(), (: 


then the lateral force per unit length is given by F’(s) and we have for 


the complex moment 


1 1 
M M,.+iM, | sF’(s)ds oF (1)—2 F(s) ds. (53) 
0 0 
On substituting (52) in (53), we have, finally, 
l 
Mv. ‘ , . ws ’ , y 
nic Arri(a,),.,—477 | a,ds+2i| S’(1)z,(1)-4 S(1)z,(1)—S(1)z,(1) | 


| O(tlog’t). (54) 
These results have been given in their most general form for a slender 
body satisfying the conditions given in §1. We shall now go on to show 


how they may be applied in special cases. 
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7. Bodies of revolution with a flat base 
Consider a body of revolution at incidence «, the incidence being applied 
vs and the appropriate potential w is 


S(s) i. ; 


in the (x, s)-plane. Then z 


g 


w—b, = =a log(z %)— 


- a w “q 


(55) 
neglecting a factor 1+ O(a?) (since we have taken the sections to be circles 
whereas they are really ellipses); w—b, is composed of a term representing 


radial flow outwards from z, and a term due to the cross-flow past the 


circular section, centre z, This may be expanded to give 
w—b, = = log z— “9 i s) 4 (56 
On the body z—z, = R(s)e", say, so that 
w—b, : (s) “log R(s)+i9|+aR(s)e-?. (57) 
Therefore 
do = 8 (8) og sa si s | log(s—oa)S"(c)do+aR(s)cosé, (58) 
0 


and from the boundary condition, or otherwise directly from (55), 


Cho _ R’(s)—a cost, (59) 
Cv 
so that 
[44 Pe<r) 
i Ov . 
r [(s‘a), BR 
Tg 2 : 
ii; a) 2)8"(o) do} R'(1)—a2 R(I)eos*s 
0 a . i 0 
< R(1)db 
2 
[S — log — 5 S AN) wf log(1—o)S "(a )dao— S(1)a? (60) 
0 


On putting this expression into the formula (37) for the drag, we have 


x 3 1 
D ] ‘ ¥ ] ” ” S’ we 
== — | | log S"(a)S"(s) dads — (1) [ tox . S"(a)do+ 
> P l es 27 . E S Co 7 J l Co 
0 0 
] 2 . . . 2 1) 
+ — [8’(1)} log ——— + S(1)a?+ O(t®log*t). (61) 


2n BR(1) 
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Thus the drag coefficient is given by 
D 
Cp = = Cp, +?, (62) 
I ve I 
 $p, U2S(1) - 
vhere Cp) is the drag coefficient at zero incidence. This result was 
btained by Lighthill (3). 
The lateral force F is given from (51) and (56) 
F 1.’ Y }.. ~ ] 5 Q’ ] Pad g 1) O(t lowt 
ag = — 42, S(1) 22, 8'(1) +22, 8'(1) +22, S(1) + O(F log*t) 
2P 
22’ S(1) = 2aS(1). (63) 
Thus there is only a lifting force in the direction of positive x, and the 
t coefficient is given by 
Lift(X ‘ 
C, el 21+ O(# log?t), (64) 
I py l =5( | ) 


result due originally to Tsien (4). 
From (62) and (64) we see that the direction of the extra force due to 
neidence bisects the angle between the normals to the stream direction 


nd the body axis in the plane of incidence. 


8. Bodies with a flat base of elliptical section for which S’(1) = 0 


Let the base of the body be an ellipse with major and minor axes of 


engths 2a, 2b respectively, and let the 








yjor axis make an angle f with the v-axis. | 
The centre of the ellipse is at 
Sa ¥ (Z,(S) xs) po 
see Fig. 3). Since S’( 1) 0, the potential y ~e 
b, is that appropriate to an elliptic aN. 
vlinder moving in an infinite fluid at rest J 
infinity with velocity . in the z- /b 
lirection. / 
y | § Cho . 1 
Now “ Do dr is the well-known 
ha CV 





° , : ; Plane s a 
expression for the kinetic energy of the 
’ . . ; Fic. 3. 
compressible fluid motion due to the 


hovement of the cylinder C (for density = 1). For the elliptic cylinder 
lefined above, the kinetic energy of the motion is 


377 r7(a? sin?8 1 62 cos?8) 











88 G. N. WARD 


Thus the drag coefficient is 
D 


” ‘ a. 9 b 9 4 9 

Cp = < Coo+( sin?p + — cos*B}+- O(#® log?t), (65 
4p, l 27ab ib a 

where Cp, is the drag coefficient for « 0, from (38). The complex 


potential (w),_, is 


a . (a+b)? 
\W), <4 = box \ (C?—c?) a i 79 3 AF 
v(G°—e*)- 4 
where c? = (a2?—b?)e?8, € — z Zp, 
; 2 eae 
haf{(a+b)?—c*}— +... (66 
C 


Therefore, 
F ; . oe ae 5 Incr? 
- 5 = mra{(a-+b)*—(a*—b?)e?'P} —27ab+ O(t log*t) 
1, U2 

2P1 
27a{a? sin?B +b? cos?B—i(a2—b?)sin B cos B}-+ O(f log’t). (67 
Thus we have both the components: a lift force in the plane of incidence 
giving a lift coefficient 


xX 


?; — — — 
Ip, U®xab 


re —_ Sinn? 
L= 2a(' sin?B + cost +. O(# log?t), (68 


b a 


and a force perpendicular to the plane of incidence, giving a force coefficient 
A 2n(1—") sin Boos 8 O(# log?t). (69 
bp, U*nab boa 

These results are valid when ab = O(1) from our condition on the 
curvature of the cross-section. However, by taking « sufficiently small 
when b/a is small, and if the body has two planes of symmetry, we cai 
keep the velocities at the extremities of the major axis small enough to 


keep the approximations valid, and we have 


oC," Sa b . m 
L) - 2 sin?B cos?8 , ete. (0) 
CX ]x=0 b . a 
When +a is small, it is more convenient to take a different lift coefficient 
C’,. defined by C7, = Lift/}p, U2a?. ete., in which case we may write 
aC; -— b? : iia - 
— E 27 | sin*8 + —. cos*8)+ O(#? logs), (71 
a a= 
aCe b\ . mn - =< 
a a | ;} sin B cos B+ O(t? log’t), (72 
CX Jy=9 ac 


a- 


, ° > b? 9 9 9 
"p = Cpo+ =(sinB ~ cos'f) x? + O(# log”t). 
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Again, it will be seen from (65) and (68) that the direction of the extra 
foree at incidence in the plane of incidence bisects the angle between the 


normals to the stream direction and the body axis in the plane of incidence. 


9, The drag force at incidence 
The result noticed at the ends of §§ 7, 8 for special bodies may be shown 
to be true generally if we define the position of zero incidence to be the 
position of zero lateral force. Let do, be the potential in the plane s = 1 

vhen the body is at zero incidence, so that from (43) we have 

20 hop dz (), (74) 

: 
ind let dp, be the extra potential in the plane s | due to incidence, so 


that the lateral force is given by 


Sh 21 | do, dz 2 | do, dy—2t | bo, 42. (75) 
ori . a oy 

If » is the angle made by the tangent to C at any point with the x-axis, 
nd if we take the incidence to be applied in the (2, s)-plane (as we may 


do without any loss of generality) and to be «a, then the boundary conditions 


give , 
CO . 
= x SIN 7 (76) 
CV 
on the contour ¢ 
From (37), the drag force is given by 
D D ' Ch 4 Cd F Cd 
0 J ai? ree 2... Ton J. 6 2 
— doo di bo, —— a1 do, —— dr+ O(t log*t), 
+P} " J O1 J Cv P Cv 
( ( oy 


(77) 
where Dp. is the drag foree at zero incidence. 
Consider now the values of the three integrals in (77). We have. for 


the first integral 


doo | dr \ doy SIN 7) dr X doo dy == @) (78) 
¢ C C 


by using (76) and as a result of (74). 


For the second integral, we have 


U 
{ J 
, CO ¢ CMD 
d 00 7_ J ~“Yo0 J rol 
Pot (l (401 Doo d1 > 
Cl l Cv 


(40 ; Pon — poo Eat) ry dé 
a Or |r=my, 


(79) 


‘rom (78) and by applying Green’s theorem, since do, and ¢», both satisfy 
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Laplace’s equation in two dimensions. We may now let 7, > %., and since 
the series for do) and ¢), commence with multiples of log r and 1/r respec. 
tively, we see that this integral vanishes. 

Finally, by using (76) and (75), 


do, Por dr x doy sin uf] dr — oa } bor dy atin =~» (80) 
, ev ; / p, U? 
C Co Cc 
Thus, we have D = Dy+4aX + O(t8 log*t), (81) 


which may be written in terms of the force coefficients 


Up Cpo T k(eC, C va 0 ” 1 O(t4 log*t). (82) 


10. Flat laminar wings of small aspect ratio 

The condition on the maximum curvature of any section of our bodies 
prevents the study of some very interesting problems, and we now con- 
sider the effect on our approximations of removing this condition in 
certain special cases. Such a special case is that of a flat laminar wing 
of small aspect ratio with highly swept back leading edges ; we will consider 
this as an illustration. 

The flat laminar wing may be taken to be the limiting form of a slender 
body of elliptical cross-section as the eccentricity of the ellipses tends to 
unity. As we have seen at the end of § 8, if the eccentricity is not actually 
unity we can calculate the expressions for C, and Cp for sufficiently small 
incidences. Now letting the eccentricity tend to unity, we see that the 
expressions for the force coefficients tend to definite values, or more 
precisely, we see that (@C',/éa), 9, (C’p)y—9, and (@2C'p/éa?), _, tend to definite 
values. If we reverse the order of the above two limiting processes, we 
find that for finite incidence the velocities become infinite at the points 
where the curvature becomes infinite. However, the solutions for the 
potential ¢,) and its derivatives are still an approximation to the true 
potential and its derivatives except in the immediate vicinities of the 
points of infinite curvature, and so the pressure in particular will be 
inaccurate only for comparatively small regions on the surface of the body. 
Since the pressure singularities are integrable, it seems to be a reasonable 
assumption that the integrated body pressures will give a true approxima- 
tion to the correct body forces for small incidences. 

In what follows, we shall use the above assumption to calculate the 
body forces when there is no restriction on the curvature of sections, 
using the formulae obtained previously, but omitting the O-terms as these 
have lost some of their meaning. 

Consider now a flat laminar wing, pointed at the upstream end and 
with highly swept back leading edges in order to conform with the 
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ondition for a slender body, lying in the plane a = —as, so that its inci- 
lence isa. Let a be the semi-span at any section s = constant and let the 
rigin of 2, ( vy tty, z+-as) be taken at the centre of the local span. 
[hen the appropriate incompressible flow in two dimensions is that for 
: flat plate with its edges at 2, ia. and we have 

w = af,|(z2+a2)—z,}. (83) 


The flow given from this expression for w is appropriate so long as the 


edges of the wing are ‘leading edges’. When the edges are ‘trailing edges’ 


this form for w is no longer correct, since a vortex sheet will be formed 
lownstream from the trailing edges: the strength of this sheet must be 
letermined from the condition of Kutta and Joukowski, that the flow 
eaves the trailing edges smoothly, and from the condition of constant 
orticity along each vortex line. Thus, in principle, we must add other 
terms to (83), representing the motion due to vortices, in order to cancel 
the singularities in the velocities at y, -_a. If the incidence is small, 
then these vortices will lie approximately in the plane of the wing (actually 
ve shall find that they do so exactly) and so we shall have (€¢9/€y,),,-9 = 9 
upproximately for |y, a. Hence the vortex lines, following the stream- 
nes, will extend downstream in straight lines. These conditions are 
sufficient to determine the strength of the vortex sheet everywhere. 

If the maximum span of the wing is 2b, and the wing is symmetrical, 
then we have to determine the flow for a plate of width 2a with two vortex 


sheets extending from Yy ato yy, _b, 

Now if q dw/dz,, then we must have Re(q) constant for x, 0, 
j;, <a, and Im(q) = 0 for 2, = 0, \y,| > a: hence dq/dz, must be a pure 
maginary for |y, a, and a pure real for |y,| >a ona, = 0. Also we 
should expect that dq/dz, will have poles at z, Lib, and we must have 
ly dz, = O(1/z?) as z, >00. Thus a possible form for dq/dz, is 

dq A 
a , (34) 
dz, (z7+?),/(23+a?) 
where A is a constant to be determined from the boundary conditions. 
By integrating, we obtain q satisfying the conditions Re(q) —a on the 
wing and q 0 at infinity, viz. 


cosh-"[(b/a),/f{(z3+-a2)/(23+-6?)}] 


-l}. (85 
cosh-1(b/a) >) 


rhe strength of the vortex sheet given by (85) is constant and equal to 


77 Xx 


(86 
cosh 1(h a)’ 
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which depends upon the value of a. The strength of the vortex sheet mus 
be independent of a@ since it must be independent of s, and so we must 
superimpose a certain distribution of flows of the above form in order to 
satisfy this condition. 

[If the strength of the vortex sheet is f(y,) at any point y,, then th 
appropriate form for dq/dz,, from (85) and (86), is 


b 
dq @ {By (6°—a*))_ f(B)dB (87 
dz, J eB\ #2+f J./(2-+<a?)’ 
which gives for the velocity Re(q) on the wing 
b b 
. 2 * 6B) dB 
Re(q) = | (cos! \ f(a) ap = | TENE , 
e(q) ( aa f(B) dB } 7(@—a) (8 


We must determine f(8) as a function of B and b only, so that (88) is 
independent of a. It is found that the distribution given by 
a a (89 
7 4 (b°—B?) 
satisfies this condition and the boundary condition on the wing: this 
corresponds to an elliptic lift distribution for the wing, and we have 
finally. from (87) and (89), 


(90 


dq _ 2« [ € (By (B°>—a*)| B dB 


dz, om J OB\ 22+ B? | ,/{(b?—B?)(z2+-a2)} 


This integral can be evaluated to give 


—_ (91 
dz, (27+-b?)! 
and we obtain finally 
w x{,/ (z27-+b?)—z,}. (92 


Thus we find that the expression for w downstream from the section 
of maximum span is the same as that for the flow due to a flat plate of 
width 2b, on comparing this result with (83). 

We may now put z, z+as in (83) and (92) and expand in inverse 
powers of z. The coefficient a, of 1/z is, for the sections upstream from 
the maximum span, from (83), 


ay sa" (da/ds > 0), (93) 
and, for the sections downstream from the maximum span, from (92), 


a, = $0? (a < bd). (94) 
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Thus if b is the maximum semi-span (whether it occurs at s = 1 with 
aids > 0. or for some s |), by using (52) we find for the lift coefficient 

Lift . 
Uy 1, 2p2 > 27a, (99) 
5p, U “b* 
referring the lift coefficient to the area of the wing planform W,,, 
, Lift 
CU, —- == 4rAa, (96) 
op; U*S,, 
here A is the aspect ratio ( 4h? S,, ). 


The drag coefficient is, from (82), 
, Drag 
( = LA a?, (97 
D [ T2q, : ) 
It will be noticed that there is no lift on sections downstream from the 
iximum span. This is due to the fact that the downwash caused by the 
rtex sheet is just the velocity required to satisfy the boundary condition 
the wing. In fact, the downwash is al’ over both the wing and the 
ortex sheet itself, so that the vortex sheet remains coplanar with the wing 


verywhere, including the portion downstream from the wing. The shape of 


the wing planform downstream from the maximum span is thus shown to be 


mmaterial, so that we need not impose upon it the condition of slenderness. 
The results of this section were obtained originally by Jones (5). 


1. Winged bodies of revolution 

The appropriate complex potential w for the case of a body of revolu- 
nm with wings of small aspect ratio may be obtained from the results 
the previous section by a simple conformal transformation. 


Consider a pointed body of revolution of radius R at any section, at inci- 


ence « to the stream so that z vs, and carrying wings symmetrically 
wed on either side of the body in the plane x vs. Let the edges 


{the wings be at a distance a from the body axis at any section (if there 
reno wings at some given section, then a R). Randa are both func- 
ns of s, but the indication of the functional dependence will be omitted. 
\t zero incidence, when the axis of the body is in the direction of the 


in stream. the flow is that for the body alone. since the wings lie in 


‘tream surfaces. At incidence «a we have to add to the complex potential 


r zero incidence terms due to the cross-flow over the body and wings. 
r} : > ‘ . . nal 
‘hese extra terms may be obtained from expressions of the form (83) by 
eplacing z, by z,+ R?/z, and a by a— R?,a in the square root: the origin 


*, IS again taken at z vs. Thus we have for w—b, 


2\2 2\2 
uv b, : SN (Ss log -] X | (2: — a + (a— . ) | “| (98) 
LTT A | ' zy a } 
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This form for w is appropriate when the edges of the wings are ‘leading 
edges’. The definition of a leading edge in this case is not quite as simple 
as it was for the flat wing alone, since the streamlines will be curved jf 
R is not constant. From (98) we find that ¢¢,/éy, in the plane of the wings 
and for y, > a is given by 


, ats 
( a = Sa (99 
CY3) 2,=0.y:>a 2ry, 


this quantity is the slope of the projection of the streamlines near to the 
plane x, 0 on to this plane, and so the equations of the streamlines 
near the plane of the wings are approximately 

yz R?+ constant. (100 


Equation (100) can be used to find the value of s for which the edge 
of the wing is tangent to a streamline, and downstream from this section 
the edge is a ‘trailing edge’. As for the case of the wing alone, there will 
be a vortex sheet extending downstream from the trailing edge, and the 
strength of this vortex sheet will be determined from the conditions given 
in §10. The determination of the distribution of vorticity at any sectio1 
for a general variation of R with s is very tedious and, in general, the 
expression for w seems to involve hyperelliptic integrals, so we will be 
content with the solution of the problem when R# is constant downstream 
from the section of maximum span. If b is the distance of the wing edge 
from the centre of the body at maximum span, then the expression fo 
w may be obtained from (92) by the same conformal transformation as 


Now toe-S 3] 


The conformal transformation by means of which (101) is obtained from 


before. and we have 


w Xx 





since S’(s) is now zero. 


(92) does not alter the geometry of the vortex system, which still remains 
coplanar with the wings. Thus the boundary condition on the wing 1s 
still satisfied by the downwash from the vortex sheet and again the plan- 
form of the wings is immaterial downstream from the maximum spad 
provided that a < 6 everywhere. There is, of course, no lift on these 
sections. These results are true only if R = constant; if R varies, then 
in general there will be a lift force. 

In general the dimensions of the body and wings, i.e. R and a, must be 
continuous functions of s together with their first derivatives, otherwise 
we get infinite d?w/ds? and discontinuous pressures everywhere at the 


section of discontinuity, and the analysis given here breaks down. This 
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condition has been shown to be unnecessary downstream from the maxi- 
mum span, and it can be shown that, in addition, it is not necessary on 
the leading edges when a R, and da/ds may have a discontinuity at 
this point. By differentiating w, given by (98), it can be shown that dw/ds 
is continuous and d?w/ds? remains finite when a R if da/ds is discon- 
tinuous there. Thus the leading edges of the wings can come from the 
body at a finite angle without the necessity for introducing further 


:pproximations 



































vortex sheets 


Fic. 4 (b). 


We are now in a position to determine the lift and drag coefficients for 
vinged bodies of revolution of two different types, (i) a body of revolution 
f general variation of cross-sectional radius, with wings whose edges are 
everywhere ‘leading edges’ upstream from the base (see Fig. 4a), and 
i) a body of revolution with a uniform section downstream from the 
section of maximum wingspan (see Fig. 45). 

If 2b is the maximum span and R, is the radius of the base in each of the 
two cases, then the lift coefficients and the extra drag coefficients due to 
neidence have the same form for both. The drag coefficient at zero 
incidence depends upon the actual body shape, of course. From (98) and 
101), by putting z, 2—z, and expanding in the form (10), we find for 


t} 


he coefficient of 1/z at s l 


1 4 ~ ¥ 
OT, YS) from (98), (102) 


from (101). 





(103) 
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By substituting these values of a, in (51) we find for the lateral force in 
both cases a lift force given by 


p. 4 b4 —b? R? + Rt 
7 = Ina (104) 
oP; U~ ” 
If we take the lift coefficient relative to the base area of the body, we 
have . 
nas Yai ) h4 b*? Ri Ri (] 
= 2a > 05 
L b?R? )} 


The corresponding drag coefficient is, from (82), 
b4—b?Ri+ Ri, 
7 

b* RF 


where Cpp is the drag coefficient at zero incidence. 


Cp Cop+ (106) 


It may be more convenient sometimes (particularly in case (i) when 
?, = 0) to refer the lift and drag coefficients to b?, in which case we have 


R2 R* 
C; featl ——24.—2 107) 
* | h2 ° i) \ 
4 as R2 R4 
and Cp = Coot “(1 — a+) (108) 


The lift coefficient given by (107) may be compared with that for a wing 
of small aspect ratio given by (95), from which it will be seen that, when 
R,/b is small, they are very nearly equal, and since R,/b < 1, the lift 
coefficient for the winged body is always less than that for a flat wing 
of the same maximum span, showing that the body is less effective in 
producing lift than a wing of the same projected planform. 

We can obtain some interesting results for the magnitude of the inter- 
ference effects between the wing and the body from a consideration of 
case (ii) when the wing root lies wholly on the cylindrical portion of the 
body. The difference between the lift L,- on the combined system and 
the lift Lz due to the body alone is given by 


Lyw—Lp = aan Ri)” (109) 
jp, 0 P 
Now the lift Lj on the wing alone (of maximum span 2 (b— R,)) is given 
from (95) by 
‘ Ly 


ro 


a 277a(b— R,)?, (110) 
2P1 

so that we have for the extra lift due to the presence of the wing in 

terms of the lift of the wing by itself 


i. ( + 7) baw (111) 
) 
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For a very small wing, for which b— R, is small, we have approximately 

Law—Lg = 4Ly. (112) 

From the potential (55) it is easy to show that the upwash near the wing 

due to the presence of the body is aU, so that a small wing is effectively 

at incidence 2«; thus we see from (112) that half of the extra lift comes 
from the wing and half from the body, approximately. 

It is suggested that the result (111) will provide a sufficiently exact rule 
for estimating the effective lift for a wing of any shape attached to a body 
of revolution. Verification of this hypothesis must await an extended 
theory or experimental evidence. 
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SUMMARY 

A method is developed for solving approximately the complete equation for 
isentropic axisymmetrical flow of a gas past bodies of revolution. The potential 
is expanded in powers of ¢ and log ¢ (¢ being the thickness ratio of the body), and 
the coefficients in this series are then expanded in powers of r and log r near the 
body, sufficient terms only being found to give the potential correct to order ¢4 on 
the body. Arbitrary constants are found by the boundary conditions. From the 
result the pressure and drag coefficients are found correct to order ¢4. 


Introduction 


THIS paper is a development of the work of Th. von Karman and 
N. B. Moore (1) and M. J. Lighthill (2), who considered a linearized theory 
of supersonic flow past bodies of revolution. 

In this paper we consider the general equation for isentropic, irrota- 
tional, axially symmetrical flow of a gas, neglecting viscosity and con- 
ductivity, past a thin body of revolution placed in a uniform, supersonic 
stream, with its axis lying in the direction of the undisturbed stream. 

The equation of the body is taken in the form r = ¢ R(x), where t is the 
thickness of the body, its length being taken as unity. The potential at 
any point will therefore be a function of t, r and x. The potential is 
assumed to be expanded as a power series in ¢ and logt, whose coefficients 
are functions of x and r. This expansion is substituted in the differential 
equation and by equating powers of ¢ and logt, equations to determine 
the coefficients are obtained. These coefficients are then expanded in 
series of powers of r and logr whose coefficients are functions of x. These 
new coefficients are determined by making use of the differential equation, 
the boundary condition at the surface of the body, and the condition at 
infinity. 

The first coefficient in the series of powers of ¢ and log ¢ is the linearized 
potential. This vanishes at infinity and by making use of its integral form 
we derive the expansion of this coefficient in powers of r and log, thereby 
using the conditions at infinity to fix the series. In deriving these series 
terms which are of order ¢*log*t or higher order} on the body are every- 
where neglected. From the results obtained the pressure distribution on 

} Terms containing the sixth or higher powers of t are neglected. 
(Quart. Journ. Mech. and Applied Math., Vol. II, Pt. 1 (1949)] 
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the body and the drag coefficient are each found as a sum of five terms 
in 2, Plogt, t*, tlogt and ¢#log*t. The pressure coefficient is given by 
equation (78). 

In actual fact the flow is not isentropic and it is necessary to consider 
to what order change in entropy at the shock-wave might affect the 
slution found. In other work [Lighthill (4), Broderick (5)], it has 
been shown that, for the case of axially symmetrical flow past a cone, the 
pressure difference everywhere is of order ¢?logt; the entropy change at 
the shock is of order ¢'*, and there is no difference in the drag and pressure 
coefficients between the complete solution and the isentropic solution if 
terms in f®log*t are neglected. It appears probable that this holds in 
seneral for bodies of revolution. 

Finally, a numerical example is worked out and curves drawn showing 
the relation between the pressure coefficients at the surfaces of cones with 
semi-vertical angles of 5°, 10° and 15°, and the Mach number; a comparison 
is made between the results found from this paper and the numerical 
results recently computed in America (6). 

[am indebted to Mr. M. J. Lighthill for suggesting the above problem 
tome and for much helpful advice and criticism in its solution and in the 


preparation of this paper. 


I, Solution of the equation for isentropic, irrotational flow of a gas 
past a body of revolution 
If V is the velocity of the uniform stream, and v(1 +3) and yo are 
OX cr 
the velocities at the point x, r when the body is placed in the stream, 
where x and r are cylindrical coordinates, then the equation of motion is 


eat} 
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vig — Pols 4 28) 8 
5 = Pi folr4 ): 


; ~.9 |? 
O21 r oxer ox cor C r| 


where c = ,/(yp/p) is the local velocity of sound, given by 


= a! v2(2re + (4) +(2) (2) 
Z | Ox Ox or J 


and pp», pp are the pressure and density in the undisturbed stream. 
Denoting V2/c2? by M2, we obtain from (1) and (2) 
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The equation of the body of revolution will now be taken as r = tR(x), 
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where its length is unity and its maximum thickness ¢, so that by varying 
t we are able to consider a system of co-axial bodies of revolution. “We 
will assume that R(x) = 0 at 2 = 0 and that R’(z) is finite at x = 0. go 
that the bodies of revolution are pointed at the origin. We also assume 
that R(«)R’(x) is an analytic function in the region 0 < x < 1. 

Since the potential is a function of t, r and 2, it is reasonable to assume 

an expansion for ¢ in the form 
$ = Pb +Pdstl byt Obs +t oot 
+tlogtd, +#6 log gf+ 
+#® log*t dg+.... (4) 
(It is important to note that an accent does not denote differentiation when 
it occurs with a term ¢,. All derivatives of terms ¢,., 4;, etc., will be written 
in full.) 

It is not apparent here why the terms in logt are included in (4), but 
the reason will appear later by examining the boundary conditions at 
the surface of the body. It is found impossible to satisfy these conditions 
unless the terms in logt are introduced. Substituting the above expression 
for ¢ in (3) we obtain 


ur — rary 4+ 26 oes 2t4 log t = +s 4 of)" 4 


“— 
(2) + o(e } | LV%4 a+ PV2d,+14 log t V7d4+t4V2d,+EV2d,+ 
+1°V24¢,+t8 log t V2d¢+#° log*t V265+ O(t7)] 
= ee x 
Ox Ox ~ OX Cx 


40° a a 


SO SF+ 
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x |e — 3 3° “F9 + tog t 
Oac* 
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= +é8 log! +48 log 3 ots ~+O(t (°)|+ 
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x jeaet+e 57°F 3 1 log 12 P41 OP 1 OE | + 
oxor oxor oxor exer 4 
+e( a) Moo, O(t?). (5) 
cr or 


+ R(x)R’(x) may be expanded in a Taylor series near each point in the interval. 
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Equating coefficients of powers of ¢ and logt in this equation we obtain 
the following equations: 


For ?: M-V%¢, _— 0, 6) 
fp: M-*V%4,—% — 0, (7) 
ar 


Hlogt: M-*V? hb, — 


fa: M-2y24, “$4 (y— 1) 2V%bs 4.992 & O21 9 Ops Ops (9) 
Ox 


Ox Ox? or exér’ 
; as a4, ad. ap 

5. M-2V2¢.——+5 y—1)—2 V2d. + (y— 1) 2 V2 

t ds Fa (> 7 bst(y ae oot 


1 92 Obs | obs Abs | obs bs Obs Mb. (49) 


+4 9 ; a 7 ie 
Ox Ox? Cx Ox? or oxer or é@xér’ 


Hts — (y—1)Bayrg + (y—1) P2944 


7 (By) v4.4 258 it 


Ox O02? 
Ob, Ob, Obs 


6: M-V24,— 


y—l1 
—- l ) Pavey, 7 oe ( ) 
Ox 2 


9 © ds c "3 1 96 py ¢ Ob. Oh. 
bs ex 








+2 


= +. 2 — 





Oa Oa? Cx Ox? ea? Cr Ox Oxor 
9 Pe c "bs . 9 Obs Pde 9 Obs @ "bs Od,\* 07h, 11) 
. + 2—— —~ +2 +") => ( 
Or oxor or oxor or exer or or 
; ae) C it ad; 
logt: M-2V94,— Pe — (y—1) Payrg, + (y—1) Se yg, + 
Cac? Ox 
9 Ops Od, | L 9 ba é P24 9 ops Od, +2 9 ob é Hb, (12) 
Ox Ox® ~ Ox Ox? Or oxor Or Oxér’ 
Slop: m-2y24" “Ps _ 0, (13) 
Ox? 


Karman and Moore (1), and M. J. Lighthill (2), have shown that the 
solution of equation (6), applicable to the case of flow round a body of 


revolution, is 2—ar 





_ { —_fody _ 14 
rs V{(a—y)?— 02} - 
K (apr) f (x) (15) 


xr 
in the Heaviside notation, with p- written for | dx and « for ,/(M?—1); 
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f(x) is an arbitrary function which, however, is zero for x < 0 and for 
« > 1, and Ky, is the modified Bessel function of order zero. Expanding 
the Bessel function we have 


dbs — > [b,(apr)?"—a, (apr)?" log(apr)] f(x), (16) 
1=0 


7 
which is a series involving powers of r and logr times a series in powers 
of r. The first few terms of (16) can be written 


K+ A logr+ Br?+ Cr? logr-+..., (17) 
where A = f(z), (18) 


K = (log 4ap+y)f(x) = f(x)log 4a— [ f’(y)log(a—y) dy,, (19) 
i) 


by the formula logx = —logp—y and the Heaviside representation of 
a Faltung integral, 


C = Jop*f (x) = Jo2f"(x), (20 
and B = }o?p*(log gap+y—1)f(x) = 4a? [K"(x)—f"(x)], (21) 


where an accent denotes differentiation with respect to 2. 

Subsequent terms in the series (17) could be similarly found, or alterna- 
tively by substituting (17) in equation (6). However, the next term in 
(17) is O(r*logr) which on the body is O(é#*logt), and thus, in the term 
t?¢, in the expansion (4) for ¢, terms O(¢®logt) have been neglected. It 
will be seen later that f(a) is determined by the boundary condition and 
hence that the series (17) is determined uniquely. The series (17) is valid 
in the region 6 < x < 1 and by obtaining it as above from the integral 
form (14), we ensure that the boundary condition at infinity is satisfied. 
The expression (14) gives the value of 4, everywhere inside the Mach 
cone x = ar while ¢, = 0 in the region in front of this cone. 

Since ¢, and ¢j satisfy the same equation as ¢, we may assume that 


¢, = K,+A,logr+ B,r?+... ) 


(22) 
and d, = k+alogr+br?+.... J 


It will later become apparent from the boundary condition at r = tR(2) 
that ¢, = 0. Using this fact and equation (10) we see that 4, has also 
the form 
k, +a, logr+6,r?+..., (23) 
and ¢, will also be seen to vanish identically by considering the boundary 
condition. 
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Now equation (9) may be written 


44 1 Oba 2b _ pos (y—1) ey 2 2 4 ogee The 
or? fr OF 5 aw CL Cx” cor oxor 
F(x,r), say. (24) 
When 7 is small and 0 < x 1, we have 
Odo Z ' ’ 
Pe - (2B+C)r+2Cr log r+ O(r log r), 
o7 i (25) 
O*dho A’ i " ” , : 
and +.(2B’+0’)r+2C'r log r+ O(r? log r). 
oxcor rT 
Hence 
24, AA’ — ' nos 
by Oe t 4 A(2B’ +0") 4+A"(2B+0)+ 
or oxo? in 
+2(A0'+A’C)log r+ O(r?log?r), (26) 
and 
Od O*h5 ( r, ‘y ae ( 2] : SK" 4”) is O 2] ’ 97 
ogi ‘K'+A’ logr+O(r? logr) K"+A*" logr+O(r*logr)}. (27) 


Hence from equation (24) we may write, when r is small and 0 < x <1, 


: r ; . i ' 
F(x.1r) --. Ylog*r+ Rlog r+ S-+ O(r? log?r), (28) 
a 
where P 2M2AA’' and we do not need the values of the other coefficients. 


We solve equation (24) by Heaviside operators and, since ¢, = 0 and 


$,/ox = 0 when x = 0, the equation becomes 
24. «1 ad 
oO C@D 
Pa Pa 22 Ne « 9 
> =—a*p*hy F(x,r). (29) 
cr” , oF 


If F(a,r) is replaced by zero then two independent solutions of (29) are 
I,(apr) and K,(apr), where J, and K, are the modified Bessel functions 
of zero order. Now denoting J,(apr) by U(r) and K,(apr) by m(r), and 
putting l’m —Im’ w, the general solution of (29) is given by 


T r 








d, = U(r) oe F dr—mi(r) | U(r) F dr, (30) 
J w w 
where w l/r, and hence we may write 
dy = K, (apr) sI,(aps) F(a, s) ds—I,(apr) | sK,(aps)F(x,s)ds, (31) 
a r 


where a@ is independent of r. The limits have been chosen so that ¢, 
represents a disturbance travelling in such a manner that 2x increases with 
', which is easily seen by considering the asymptotic expansions of J, and 
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Ky. Also ¢, is zero in the region x < ar. To find the form of 4, when r js 
small we have 











I,(apr) = 1+ O(r?), 
> 2 (32) 
and K,(apr) = (ioe = —7) +00 log r). | 
apr 
Hence, neglecting terms of order r*log*r, we have from (31) 
= —— ds—|log —-— y] C— 
= (eS IE : (oe 7 
oe 
a { Koper? ae { Koons) [F(w, 3) 5} ds, (33) 
8 s 
r 0 
where C is a function of a. 
Now 
— fe (aps Fc 3-3 ds— | K,(aps) = ds 
8 8 
re) z/a(coshu+1) | 
= — | a | ee cosh nn ) ds— 
8 
0 0 
cosh~4{(z/ar)—1] 2/a(coshu+1) P h 
- | du es ds, (34) 
8 
and 
cosh~4[(x/ar)—1] x/a(cosh u+1) p h 
: | du i} P(z—as cosh) g, 
s 
cosh~1{(z/ar)—1] zja(coshut+1) , a Pe a i 
. s” cosh”u 
— — | du | ae - here aoe — ds 
0 r i 
cosh~*(x/ar) —1] a , P(x) h" 
~_— du| P(x) log i nk 
| “| (*) ar(c = u+1) a ee n!n(cosh u+1)" 


+P Sree (35) 


where P(x) denotes the nth derivative of P(x). 
When r is small 


9 
cosh = — 1) ~ Sai seniors + O(r?), 
ar x 


ar 
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and therefore 


cosh~*{(x/ar)—1] 
cosh"u 
(1+ cosh u)” 
@ 


2x cosh” u—(1+cosh u)" 
t-te OO 





Also 


cosh~3[(x/ar)—1] 


* x 2 
— P'(x)ar coshu du = — Paar, |{(=— 7 _ 7 
or 


0 
= —xP'(x)+O0(r), (37) 
and 
cosh~*{ (x, ar)—1]) 
deb | P(x)log 


. 


0 


x 
vr(1-+ cosh w) 
cosh ~1[(x/ar) —1] 
» or 
= — P(x)log = [log = — on 0) | + Pe al | (u—log 2) du+ 


0 


+ | {log(1+cosh w)—w—log 2} du— 


0 
in | {log(1+cosh w)—u—log 2} il 
cosh "{(z/ar)—1] 
. oy . . o 
= — Ptxjlog * [log — +. 0-4) | + Pte)| dost — "jog = + 
x7 | xr x ar 2 or 


+ O(r? log r) log 2log = + “log 2-+0(r2 )+ dn? + O(r | 
ar 


il 
= Pee) dor? — 3 log? = + O(r) F (38) 
! xr 
Finally, from (35), (36), (37) and (38) we deduce that 
cosh (x) xr) —1] x x«(cosh u+1) P 
a | i | (x—as cosh U) as 
x 8s 
0 r 
9 22 
= Ple)| in? 4 log? == |—xP'(w)-+ r)+ 
r 
Hn) \n ox ny —(1 
4 5( yn = ™ z 4 {= cosh"u— ( +-cosh w)"  dul+0(r). 
= nn xr (1+ cosh u)” 


0 


(39) 
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Therefore, neglecting terms in r* log*r, we write from (33) and (39) 


9 2 : 
& = — (oe aa 7) log r—loga)P— (oe = -7¢ — 


0 z/a(cosh u+1) 





P(a—«s cosh u 
— | du i} s| Fle—as cosh u, s) —- | — ; ds+- 
s* 
0 
; out a 2x ™ 
+ P(x)| 42?—4 log? — — 4 log*r + log — log r | —aP’(x)+ 
xX x 
8) 
P(x)” 2a 
—])r+1—, - |log — —log r+ 
T > ) n!n ? ” 
n=1 
ie 8) 


+ Ss i +O(r), (40) 


(1+ cosh u)” 


and hence ¢, has the form ° 
x, log*r+-B, r? log*r+ O(r* log?r) + 
+y,logr +6, r?logr + O(r4 log r) + 
+e, +n,7r +O(r*) (41) 
(since the terms in r?logr and r must vanish as ¢, satisfies equation 
(24) when F(a,7) is given by (28)). 
From (40) and (41) we deduce that 


y = oP, (42) 
Sin = (n)/4*\~n 
y, = C—Ploga+ P(y+log tap)+ P log = — > (—1)" no Ee 
: n=1 , 
- C—Ploga, (43) 


from the interpretation of the operational forms and since 


- = I(n)/ ~\ yn 
P(x)log x— | P'(x—y)log y dy (eae (x)a" 


n!n 
0 n=1 
and 
o x/a(coshu—1) 
e, = (C—Ploga)(log 4ap+-y) — [ au | 8| F(w—as cosh u, 8)— 
0 0 


P(x—as cosh u)) 


or 
3 | ds+ Ptx)(int—4 log? == |—eP '(a)+ 
s* rs 
[ Pwlogte—y) ay} 


e 


0 


+ Pesto x— 


n\n (1+ cosh u)” 


_ coe} 
+> (ye Se | ee te. (44) 
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Eliminating C and a from (43) and (44) we find that 


. 


24 
é y,(x)log 4a vi(y)log(x y) dy+ P(a A(t 7 — 4 log? — “| —2P (a7)— 
x 
0 
r «(cosh u +1) 


— } du s F(x— as cosh u, s) — 


0 0 


s* 


P(a—«as cosh 1U ) 
. ds+ 
7) 


P(x)log x { P'(y)log(a—y) ayjiog che 


0 
: , cosh u 
|. du | P'(x—y)log y dy+ P(x)log- ses _— 
4 7 1+-cosh u 
0 rcosh u/(1+ cosh u) 


P( x Jlog vee |. (45) 
1+ cosh u 1+-cosh u 
Hence in the solution (41), a, and e, have been determined and later y, 
will be determined by the boundary condition at the surface of the body, 
thus making the solution unique. The terms neglected in (41) are 
O(#log*t) on the body, and the corresponding terms in the expansion for 
dare O(t* log?t). 

The next step is to solve equation (11) for ¢,. Using the fact that 4, 
nd ¢, vanish identically, (11) may be written 


7 O*dhe Od, 9 Od, ‘ , y—1 Od, : ‘ 
MV", “5 (y—1)—V*¢,44 (y—1)- ~ V7. +! aa — V?b_+ 
cs On OX ys On 
- oh ee “er 
a= (ea “V2ho- 2 Pe , 9 Ma Oba Os Vee 
2 \ or Cx ¢ a ; Ox Ox? Ox Ox* 
1.9 os Ody Oho 1.92 O*h, -. 9 obs C be, (06 } 2¢ “be (46) 
Or Ox OxOr’ Or OxOr- 8 Or Oxor or} er? 
™ 5 Obs ’ : 
Since V?¢, = M?—*, (46) may be written 
Cx? 
od, 1 a 20 > Od» ‘ Od od. 
a4 5 (\y 1) M* 2 V"*b4+(y—1) M4 4 —P 
or r o7 Cac? ox Ox Cx 


ax ea? Cx Ox? 


— on “Ps I Y PP she rf Ody € Obs , Ob ca)+ 


ar(ibe be 4 he i) 


or oxor Or Cxor 
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which, 
Now we have iby 
2a Od, ) or® 
V*d, = =z + O(log*r), Pat = O(log?r), | 
= 
Ob, 204 aa v1 2 | 
dxer 7. ~ + O(rlog r), | 
py em 2 | ) 
= = A’ logr+ K’+O(r? logr), (48) 
a = A + (2B+O)r+2Crlogr+O(*log?), | . 
| 
” where 
- —4 + (2B+30)+2Clogr-+O(r*log r). | and vi 
, To 
Substituting these values in equation (47) we obtain ; 
6 = 
Cy 1 Oe alte 
or ' f & ea? 
Then 
= (y—1)M*{A' logr+ K'+ O(r* log r)} = 14. O(log?r)} + 
+(y—1)M*O0(log*r) 4 M*O(log?r) + ont 
yA ep, 
" 7 
=m + O(log r) la" logr+K"+ O(r* log r)}-+ 2M?O0(log*r)+ , 
ent 


+24" logr+K’+ O(log) 4 + O(log) [4+ O(rlogn}+ | mh 


45 9 m2) 2% 108? M8" 424-4 O(rlog*n} |= + O(rlogn)| + The 
\ . J\r J it ca 
(A? othe 

+ M*\— + 2(2B+C)A+4AC log r+ O(r* log*r)} x 
- = (2B+C)A+4AC log O(r? log? ) = 
A tern 
¥ | 5+ 2B+30)+2C log r+ O(r? log n+ infi 


2a; log r 


+202 + O(r log 7) [cer +734 O(rlog? r) (49) 








(48 


sr) 4 


(49 
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hich, on collecting terms, becomes 


be 1 Og = 2 Obs 














oe |r ér ear? 
= 4 “|” ~1)M?2a, K’ 4 2K" 42M2AA'K' + 
2M2Ay’,4 ita ieahlacensdailaiiais 
+28 ly 1) M?2a, A’-+¥—* M8A2A" +2M2AA*44MPAG + 


+4M?A'a,—2M*A*C| + O(log*r), 
P,Q. , R,logr 


ee ook hoe 4 ioe. 50 
r3 r2 : r2 . \ = r) ss 
where P, A? M* = — M'f?, (51) 


and we do not need the values of the other coefficients. 


To solve equation (50) put 


$, = 2+ 8,+ O(r?)++y, log r+ O(r? log r) +8, log’r+ O(r? log?r) + 





v2 
+e, log*r+O(r? logr). (52) 
F de 2a 
Then Ls = + O(r)+ 72+ O(rlogr r)+ 
Or a 
or , - 2 
- 25, logr | O(rlog*r) +. 382!°8" + O(rlogsr), q (53) 
r r e 
and 
C2 6 9 Vo 25. oT 2 9 6 9 7 
Ve 7 - _Y2_ do a rr 2 4 oer 3¢€, log? aE. O(log’). 
er r 7? 2 r2 r2 r2 J 
0 l . 2h, 4a, , 20, , 6e, = 
Hence obs 4 Oe 25 Pe a hp ae — 2 <2 +. O(log*r), (54) 


or? r or or ptt 7 
and, comparing (50) and (54), we obtain 
> ¢ o > - 

405 P,, 26,—a*a0% = Qe, Ge, = R,. (55) 


These equations determine a, 5, and e, uniquely. By proceeding further 
itcan be shown that y, and f, are completely arbitrary, and that all the 
other coefficients are determined in terms of y, and f, and known 
quantities. Later we shall see that y, is determined by the boundary 
condition at the surface of the body and f, could then be obtained in 
terms of y, as «, was obtained in terms of y,, by using the condition at 
infinity. The only coefficient in ¢, which we require is a, which is given by 

P, M? f* 


Xo —_—_ = 


~- 56 
2= 3 j (56) 
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Terms neglected in (52) are then O(log*t) on the body and the corte. 
sponding terms in ¢ are O(#* log*t). 
Equation (12) for d, may be written 


or2 ' r Or ‘ Cac? Ox Ox ox? 


no 4° ¢ , ‘ 
o* 1 é 2 9 OP; C 
Oe , 1 Ope “bs _ = (y—1)M? bay ag' 4 (y—] 1)M4- by hs. 
o C C 0 C 0 
4. 2M? — a} be +32 M?2— Ope @ aS 9M2_-12 pe O Oba 1 9 M2— "be Oby (57) 
Ox ar ér Oxor Oxor Or 
From the values of ¢, and ¢, already found it is clear that the right-hand 
side of (57), when expanded in powers of r, has the same form as the 
right-hand side of (24), and hence 
dg = a3 log*r+ O(r? log’r)+£,+ O(r?)+y,logr+O(r2logr), (58) 
where, exactly as in the case of ¢4, ys is determined by the boundary 
condition at the surface of the body, and £, by the condition at infinity, 
and all the other coefficients are then fixed uniquely. Since none of these 
coefficients are required in the following work we shall not determine 
them. The lowest order term in ¢ due to ¢% is O(t®log’t) on the body 
which we neglect. 

Finally, the solution of equation (13) is evidently of the form 

¢, = K,4A,logr+ B,r?+C,r? log r+ O(r4 log r), (59) 
where A, is determined by the boundary condition at the body, and K, 
in terms of A, by the condition at infinity. 

This completes the useful information we can find from equations (6) 
to (13). To determine the coefficients in the previous series which are 
still arbitrary it is necessary to consider the boundary condition at the 
surface of the body r = tR(x). This boundary condition is 


é¢/er . 
$ = tR' (x), (60) 
1+ éd/éx 
when r = tR(2). 
Substituting the expression (4) for ¢ in (60) the boundary condition 


becomes 
| phe 4 pOPa 4 jog 104 4 Ps 5 ps OOs 4 oOPe 4 
or cr or cor or or 
oy : ‘ e ] : 
+ log t oe + t§ log*t ai + . = 
or or 


[+ . j2.0P2 the 1 9°31 lope a of prOha , 152 Os 6° be 


Cx Ox Ox Ox Ox 


+t log t oe + #8 log*t O$¢ +. | tR’(x), (61) 
Oa Oa 
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which, on substituting the values for the potentials, becomes 
P| A 


ar 
oa 
+Hlogt}" + O(rlogr) pa} 20 FOBT 
: 


| (2B+-C)r+2Crlogr- O(r*logr)| +0[=1+ O(rlogn)} + 
y 


4” +O(r log*n) + 
r 


r 


\r J \ 


) 
1 ¢8 ei 204 O(r log?r) 


or 
+ 1%4 O(rlogr)| +46) __ 2% +t 40(° =)h 4. 


\ 4 
r il 
+ log|= 24 O(r log r) }+. |= 
+(l -#{K’+ A’ logr+ O(r? log r)}+ {Aj log r+ Kj + O(r? log r)}+ 
+t log t{k’ +a’ log r+ O(r? log r)}+t{aj log’r+-y; log r+ €,+ O(r? log?r)}- 
t®{a; log r+k,+ O(r? log r)} +19(23-+ O(log*r) — 
y2 
| t8 log t O(log?r) +-t® log*t O(log r)+...] = tR'(x). (62) 


Now, putting r = tR(x) and equating coefficients of powers of t and 
logt in the equation (62), we obtain the following equations: 


Coefficient of t: A RR’, } 
t2: Ay 0, giving ¢, = 0, 
R 
aa — 2a,log Ry, 2a 
t (2B+-C)R+2CRlog R+- —_ + Fe — 
AK’ AA'logR - + (63) 
= . 0, giving 7, 
R R 
( 2 ‘ 
Blogt: 2CR+" Oe 0, giving a, 
Rk’ R R 
{4- ay 0 giving , = 0. 
R tan) ds j 


lt is unnecessary to write down the remaining coefficients. However, we 
note that the quantities y,, y, and A, are determined by equating to zero 
the coefficients of ¢°, log t and ¢° log’t respectively, and thus the boundary 
condition at the surface of the body is satisfied if we neglect terms in 
loos, 

Une further point has to be considered here. It seems possible, a priori, 
that by continuing the expansion for ¢ further, a term such as f°¢, might 
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be of order ¢* on the body, in which case an approximate solution for 
¢, would have to be determined. To investigate this, assume that a longer 
expansion for ¢ has been substituted in (3) and consider the coefficient 
of #8. Clearly the important terms are 





Ob, 9 be Od, Oh, oe log r 
()'< ere 7 or or ore ° re fy (64) 
a2 
Therefore M-*y24,— “Fe = (6), (65) 


and hence ¢, = O(r-*logr), and #846, = O(t®logt) on the body. 
If we consider the coefficient of ¢!° in (5) we deduce that 


M=V74,.— —- - o(5) 


and hence that ¢#%},, = O(t*) on the body. Similarly 


M~V*4,.— é ots = = of°§ "| 





76 


and therefore ¢#*¢,, = O(f®logt) on the body. From these the following 
results can be proved by induction: 


Pap-2 — O(r-*P +2), Pay = y - O[r (r-*?+? log r) (p = 2), 66 
) 
¢, = O(logr), d, = O(log?r). - 


Again, it is easy to show that ¢, = O(1/r?) and 4), = O(r-*logr). These 
results suggest that ¢{ and ¢,,-», are of the same order and this can 
be proved by induction. Finally, a complete statement on the orders of 
magnitude on the body of all the terms ¢?"log?t 6) is contained in the 
following equations: 


2" log?t dh?) = O(t"+”+2log?+4) (for n+p even and p < n—2), 
= O(t"+?+1 log?t) (for n+p odd and p < n—1), (664) 
2" log”-14d"-) — O(t?" logt), | 
t2” log"-2¢ d&™-2) — O(t?" log’t). | 


Since p is always less than or equal to n—1 we immediately see that no 
term #2" log?’téY) can contribute a term of order ¢* or lower order to ¢ 
when n+p > 3, and hence the only terms affecting ¢ to order f*, are $5, 
d4, $, and 44, all of which have been considered. Hence every term in 
the expansion for ¢ of order ¢* or lower order has been found, and all 
terms of order ¢®log*t or higher order have been neglected. 
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II. Evaluation of the drag coefficient of bodies of revolution 
The drag coefficient C, may be defined as the total drag divided by 
D . Q A 
the product of 3p) V* and the maximum cross-sectional area of the body. 

Therefore 
1 
, 4 [ p—Do : , : = 
Up ; C f097t?RR’ dx 8 | “ re dx. (67) 
dPo V* 


0 0 


y find the value of the pressure coefficient (p Po) Do V2 we have from 


uation (2), 
p -_ 4 l | 29 oP ‘ er a: ‘Od ‘| (68) 
p " Po Z { Cx Cx or } 
D Po | y72 Po 2% (2 Pas op |. (69) 
Pe By | pol ax * \ex) Var) J 
since 7 5 69) be ymes 
—| ; | | | » Po (od | od , 1 od *) (70) 
2 9 j.{—} -j—} }, i 
l 1 hence 


> | ) d Od : d yKy-)) 
ly ) Polo § + ( ‘| + (2 | } (71) 
yy p | Oxf (Or \ ¢ r | 


h on € xpansion become 


Po L Po pala! i j | (2) '\ + 
vo Z My \ Cd \or | 


nd therefore 


M | ( db \2 db f d 2 db . d 3 
| 149 (: \"  ( 14 (2 (73) 
t | \oa — Ox \or or | Cx 


On the body é¢/éa O( logt), so that we may write from equation 


Pe __ nih 5 af d)*_ (28), yn (68), 
Py V2 O2 C2 \or} — Cx \cor 


4 


co? 


: Cd A . ‘ 
ye... () + O(t®logt). (74) 
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On substituting the expression for 4, (74) becomes 


4 ‘ , ‘ ft ‘ 
pT Po — — 9f 72°F2 4 Jog 1 OP4 1 ps Oba | 1626 1 QUSlog tlog2r)\ + 
5 Do V2 Ox Ox Cx Ox P B&B 
‘ J ¢ , f y , 2 
eS yess #4 mete of. Pa -+¢6 Pe -- oOrts logt log*r)| ae 
Cx Cx Cx Cx 
— {2 P24 hog 12b4 4 22h 4 p Ps 4 O(log 1 4\o7 
| or. or cr cr \ r | 


do 
—- me Ee al O(tlog tog r) x 


Cx 


ad " o \)2 
x fp2Obo s 22a pou 4 Fa w As 4. 0(tlogt-2") + 


| or co} cr co? r 
M2 . ad 44 oo r\)2 
4 M } ,2° Po | 0 ‘log? | é (75) 
4| or r 
and, substituting the series for the potentials in (75), we find that 


ff — 25 K'+ A’ logr+ B’r?+ C’r? log r+ O(r4 log r) 


2t* log th’ +a’ log r+ O(r? log r)! 
oar.’ 2 , , J | 2\) ) ef Xo | ( 3y ! 
2t4{ a, log*r+y; log r+-e,+ O(r? log?r)}—2t \72 7 (log r+ 


Vt K’?+-2K’'A’' logr+ A” log*r+ O(r? log?r)!— 


ae 9A(2B+-C)+4AC log r- O(r? log?r) | — 


r- 


Aa A [y; , 2a, log r 
2 log t : Ollog*r)| 946) ? 3 “ | O(log*r)| — 


a ? 


- \7 


ou} 4 lf 2% (- 
2t LQ log ~ ae = EF 
lr . ah 


2 2 1 
M2t6!(.A' log» KA = Ohana t° A f 
| | 4 ? 


M*t?| A’ log r+K'+ O(7? log r)| 





yo) Y1 2a, log r 


a ) . 
+ Olr log*r)| +?®logt | + O(r log r) + 
, 


\7 ; 
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M2 [(4A3 tA* 
M"t| (AZ ope 1" 2AC log r+ O(log r)| + 
aa ie 2 J 
y2lA® | o(-)\| 20 1087 1. Orlog?r)| + 
\>3 \) r ili 


ye? 


1 1\)\ {a ) 
1/2 log t! O +. O(rlogr)} + 
8" (5) [pt O@los”) 4 


18 3A2(2B+C) . 6A?Clogr 


44)° 
19 } 7 


- {2% o(})\ 76 
nt Ce ys ! (7) \ ) 


It is to be noted that the terms in square brackets in equation (76) are 


+ O(rlog)} x 


order (® log’ (or higher order) on the body and are therefore taken to 


e negligible. They have been written out because they contain squares 
nd higher powers of 1/r and therefore might possibly tend to infinity 
ry tends to zero. However, their coefficients contain powers of f(a) which 
re O(r) as r > 0 on the body and by examination we can see that the 
efficient of oo’t in (76) is of order log*r. Hence. the corresponding 

oral in (67) converges giving a term genuinely O(f* log) in the drag 
efficient. 

Therefore. from (76) we have, on the body, 

K'+2A' log R+~_|—2A’flogt 
R2 


; 
AiI2B'R 20" R? log R+-2a, log? R 2y', log R- 


—=__ 4?(K'?+-2K’A' log R+-A” log R)- 


Jey R 
IA(2B+0C) 1AC los R- 2Ay, 1, tx, A log R 
R? R> 
ta, A -_,A2 M?A* 
= M*(A' log R+-K' _ _ 
K —s Re aR | 
{| 26 R*+-2 2a' log R ba} log R- 2y1 
2Aa 1a, A4 M?A’A? 
2K'A 9A'2loo R)\+4AC+ = og _ 
sis : RR R? | 


f 2a, Ya v2 A “a O(t® log*t). (77) 
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Substituting the values for A, B and C this expression may be written 


P Pe. . | 5 +2f' log R 2K" 2f't2 log t 


SP) V® Rk? 


M2 ft. 


| —qatprm—gargy" EL 
| (hemes ee Jlog R—o2f'? log? R-+4o2 R2K" 4 
= 


-, M2f2K’ vm , 
vf K AT aay K' log R—o2K’24 


foc, log R , 
Re" 


, 2fy a ? 2D ae > o.? day f Qa5 
+p 20, log? R+-2y, log R4 2 — 4 
—t*log t} 40? f” R?+- a2ff” — “as — 2a°f"* log R—2a?f’K'+ 


9 ft, . 
=I oa" log R-+ 2k’ 4+ 4h 


22 | r R2 


+ 4ai log R4 aA , 


t* log*t| —a®f’? + 2a’ + 204 |+ O(#8 logs), (78) 





where 


r 
- 


J RR’, K = f(x)log 4a— fT (y)log(a—y) dy, 


0 
4 = MT, X2 -7 if", P(x) = 2M?ff’, 
V1 ba®f” R?+- ff’ log R—1.2f” R? log R-+-fK’ 
—to? R?K"+-20,/R?—2 x, log R, 


r 


y1(x)log da v1(y)log(a- y) dy 


+ P(a) {d?— 4 log?(22/a)!—axP’ (x) 


P(a%—«s cosh u)) 


du s F'(x—.as cosh u, s) — ds 


9° 
sg? } 


x a (79) 
, \ tov 22 
T (x)log «— I (y)log(a y) dy jlog- oes 
. x 

0 


x r 
re .Y 


| | P’(x—y)log y dy + P(x)log 


* 
0 x cosh u/(1+ cosh uw) 


cosh u 


l | cosh u 


x x cosh u 
— Le ——— ; 
1-+-cosh u 1+-eosh u 
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. . es, CDs C 24, 9 C do 7d, 
F(a. 1) {2M?+-(y—1)M*} — 2.4. SAPS? 8, 
? C2 COx= — cr oxor J 
The pressure coefficient is given by equation (78) and the drag coefficient 


s found by substituting (78) in (67). 


III. Comparison of the analytical and numerical solutions for the 
pressure coefficient of a cone 
Using numerical integration, Taylor and Maccoll (3) found solutions for 
supersonic flow past three cones, and recently tables of solutions for 
supersonic flow past a variety of cones have been published in America (6). 
[hese are the only solutions by exact computation for supersonic flow 


past bodies of revolution which have been published and it is interesting 
to compare the approximate analytical solution developed in this paper 
vith those of (6), in order to judge to what extent the present theory 
ry be appli d to general bodies of revolution. 
For a cone ¢ is the tangent of the semi-vertical angle and R(x) = 2. 


We find the following results from equation (79): 


f : ; 
K i log = On rio r 
, = M% 
X9 pM Pa 
1 i l s x ly Ma log l 
a r(1—2M?) | 
7 | . 
h r(1 2? M?)log J r—(] 2M?\(xlogx £3. f (80) 
Fiv,r) — {2M 1) M3 — eosh-2 4. 9M2=, 
nj (O° — FT) xr 
€ WV vr log?a (3 log da vlog x 
ri a* log? Sa 2M?— })log 4a+2M?-+-}} 
2M2+(y—1)M4 
ae - s. 
2a? J 
Substituting these results in (78), we deduce that 
C, = —12— 21? log Hta-+3a%! log? Hta-+ (5M2— 1)t#log ta} 
' . y+1)M4 
sofa CEA) oy 


neglecting terms in ft® log*t. 
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We call the terms in # and # log? in equation (81) the first approxima. 
tion to the pressure coefficient and the whole expression the second 
approximation. For comparison, one other form of the pressure coefficient 
for a cone is of interest, that one found by solving the linearized equatior 
for flow past a cone, satisfying the boundary condition exactly and using 
the exact form of Bernoulli’s equation. Doing this we find that 

(i sas } _7%—! pf 2? cosh" L/at) 

, yM* 2 (Jd vt?) +-t? cosh-(1/a#) 


| f? cosh-1(1/at) ) 
(1 y7t?) +t? cosh—}(1/xt) 


. ( t,/(1—a?#?) )| yKy-1) 2 
-L — . (82 
\(1—a*l?) + cosh-(1/at)) J y M? 


/ 


Comparison of the numerical solutions (6) and the analytical solutions 


07; 
6H . : 
067 . ~ Numerical 
| ~~ - 2nd. order approximation 
\ Linear equation solved exactly 
Ist.order approacimation 
05+ ‘\ 
si ccna, 
LT 
LL 
Cy 
a ee Sa —— 
10 M 2:0 3-0 4-0 50 


Fic. 1. Graph showing the relation between the pressure coefficient at the surface 
a 5° cone and the Mach number. 

for pressure coefficients at the surface has been made for cones with 

semi-vertical angles of 5°, 10°, and 15°, and the results of plotting ( 

against MW are shown in Figs. 1, 2, and 3. 

In the case of the 5° cone it is not possible on the scale of the figure to 
show any difference in C,, between the numerical results and the second 
analytical approximation up to a Mach number near 5. Indeed, the 
difference between these two values for C,, is about 0-4 per cent. ata Mach 
number of 1-04, and about 0-04 per cent. at a Mach number of 3. 











SUPERSONIC FLOW ROUND POINTED BODIES 
OXIN 0-2r — Numerical a 
secol l ~ 2nd. order Spproscenacion 
oat | Linear equation solved exactly 
Ist. order approscimation 
ua 
| ™ 





7 YY 20 30 4-0 50 





9. Gi wing the1 between the pressure coefficient at the surface of 
a 10 e and the Mach number. 
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Ma . & ving the re n between the pressure coefficient at the surface of 


and the Mach number. 
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For the 10° cone the second approximation is good up to Mach number 
about 3, while in the case of the 15° cone the agreement is not good except 
for very smal] Mach numbers. 

We can state that the second approximation is a great improvement 
on linear theory for thin cones and the same is no doubt true for all bodies 
of revolution. 
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SUMMARY 


An approximate analytical tion of the problem of supersonic flow past pointed 

es of 1 ition was recent given by the author (1). In deriving this solution 

irrot onal flow w assumed and no account was taken of the shock- 

It dor t em possil to extend the solution when the shock-wave is 

to account in general, but for the particular case of axially symmetrical 

the extension n be made. The purpose of this paper is to show 

ution the same approximate value for the pressure 

fficient that found in the earlier paper for flow past a cone. We may infer 
t this result holds for all cor x bodies of revolution. 


Introduction 
Iv a recent paper (1) the author developed a method for solving the 
non-linear partial differential equation for isentropic, irrotational flow 
past bodies of revolution, by expanding the potential in powers of ¢ and 
gt (t being the thickness ratio of the body), and thus reducing the non- 
linear equation to a series of linear equations. It was found that the 
solution obtained differs from uniform flow only inside the Mach cone and 
ence can give no approximation to the position of the shock-wave; 
furthermore, it can be shown that the series diverges on the Mach cone. 
However, it appears to give correct results near the body and the pressure 
ind drag coefficients were found to order f4. 

It would be desirable to extend this previous work in order to find a 
solution which would hold also in the neighbourhood of the shock-wave 
and which would give the position of the shock-wave and the pressure 
ind drag coefficients when the shock-wave is taken into consideration. 
One would then know definitely whether the solution already obtained 
lor the drag coefticient to order t4 was correct. However, it does not 
ppear possible to carry out this work for bodies of revolution in general. 
Hence the extension has been made in this paper for the particular case 
i axially symmetrical flow past a semi-infinite cone, as this case can be 
dealt with analytically. We find that the pressure coefficient for a cone 
obtained by using the complete equations is the same as that found in 
the earlier paper to order ¢*; a first approximation to the position of the 
shock-wave is also deduced. We may infer that the expression for the 
pressure 


oefficient given in the earlier paper is correct to order ¢* for all 
convex bodies of revolution. This paper has not been written in order 


[Quart. Journ. Mech. and Applied Math., Vol. II, Pt. 1 (1949)] 
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to determine a solution to the problem of supersonic flow past a cone, 
This problem has already been considered numerically by Taylor and 
Maccoll (2) for three particular cones, and further exhaustive numericaj 
results have been published in America (3). The first analytical approxi 
mation to the position of the shock-wave has also lately been found by 
Lighthill (4), who used another method of solution. 

Comparisons of the analytical solution and the numerical solutions are 
made in the author’s earlier paper to which the reader is referred. 

In conclusion, I wish to thank Professor Goldstein and Mr. Lighthill 
for their very helpful advice and suggestions. 

The following notation is used throughout this paper: 


V velocity of main stream. 


Ox 


¥ Cd : : . 
J (1 velocity parallel to axis of cone. 


od : ; — 
} velocity perpendicular to axis of cone. 
or 
Ch Cd 0d ed 
‘ values of —, at the shock-wave. 
€ vy C ry CXL cr 
c = local velocity of sound. 
Cy = velocity of sound in main stream. 
y 
M 
Co 
9 
x x ( M 1) 


iv Mach angle. 


n angle between shock-wave and axis of cone. 
t tangent of semi-vertical angle of cone. 
l 
Yy tan p. 
X 
| 
3, = tan 7. 
Q / 
fied 


I. Hydrodynamical equations 
The general equation for the flow of a gas with axial symmetry in 
cylindrical polar coordinates is 
C* 24 (1 o> eas 1 OF OP OP of) (1 
V2 cx} Cat ox} exer or \er cr- 
For flow past a cone the variables such as pressure and density will be 


functions of r/x only. Thus ¢é¢/éx, and éd/ér, are constant along the 
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shock-wave and we may put ¢ = af(r/x). Then using the Rankine shock- 
wave equation, Bernoulli’s equation, and (1) we deduce that f satisfies 


the equation 
M22" 6 ff) (ful PPL ty HT 


ty¥i+f—yf')—f' PF", (2) 


vhere rly 7 and dashes denote derivatives. 
The boundary condition at the cone ist d, (1 $,) rye t which 
eeoOmMmeSs r tf (fF) (] ! (2) f(t), (3) 


nd using again the shock-wave equations we deduce two conditions at 


the shock-wave which may be written 


f(B,) 0, (4) 


nd { (P,) ;(B,—«,)+ O(B8, —a,)?. (5) 


1 (1 x3)? 
We have then to solve equation (2) with boundary conditions given by 
uations (3 1), and (5). Since equation (2) is of the second order the 
undary conditions determine the relation between 8, and a. 
II. Solution of the equation of motion near the body 


Following the method developed in (1) we substitute 


(?2f,+t4f,+ tf, f,, flog Matfe, flog $ta+...+-fgst® log? dta 
(6) 
equation (2), and equating powers of t and log $ta we deduce a series 
f linear equations for the functions fy, fy, f4.-.... The solution for f, is 
{| | 2..9\\ 6a . 
A} log v (1 ey) v( v-y") |, (7) 
| yx 


xcept for the addition of an arbitrary constant. It will be seen later 


m the boundary conditions that this arbitrary constant is zero, so we 
glect it at this stage in order to save unnecessary terms in f,. The 
bitrary constant A will be fixed by the boundary conditions at the 
surface of the cone. The solution for f, is also found and f, and f, are 


en expanded in powers of y (since y is small near the cone). The results 
1{ 1 +log dya—loPy?+ O(y')| (8) 


W*(y 


Pr 
ax” 


1) , 
| log hyo O(y- log Ly x) 


BS) log by O(y*)}, (9) 


1 tis the value of y at the boundary of the cone. 
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where B is another arbitrary constant. Making use of these expansions 





we then deduce ued 
M243 
Is —— + O(log? Syx), (10) 
4y* _ 
and fa = —Cl1+log dya+ O(y?)], (11) 


where C, again, is an arbitrary constant. 

[It can then be shown (1) that no term, other than fy, f,, fg and f,,, can 
affect the value of f to order t*, so we need not consider any other term 
in the expansion (6) in order to find f correct to order t* near the body. 

From the solution (7) it appears that y = 1/a (which gives the position 
of the Mach cone) is a singularity of f. We next consider the behaviour 
of the solution near this point. Using (7) and, putting z = (1/«)—y, we 
deduce 

fa(y) fo (2) Aw2 vizif1 + O(2)}, (12) 
where the positive sign of the square root has to be taken to agree with 
equation (8). Then from a consideration of equation (2) we can see that, 
to obtain the asymptotic form of f;,,,.. near z = 0, we have only to con- 
sider the equation 


‘ n ad (y T 1)M4 : even » 9 
22fon+2(2)—Sonsa(2) = - vi > feo) fon 2-2p(2). (13) 
p=1 
From (13) we find 
” y-+-1)M4A? 
fale) = VIA + 124), 


which suggests in general 


re A (vy 1)” 1,2 Vrn 1) > 1—3n — 
fon(2) —_ ks a + O(z?-*"). (14) 
v6 y(n 1) X 
Substituting (14) in (13) we find 
n F 
NAnsy = > A,(1—4p)Ansi-p- (15) 
p=1 
In a similar manner the form, near z 0, of other terms in f’ such as 


fi, and fg, may be determined and the solution near z = 0 may be written 


x 


A " A 2/,,—1 1)” Lyf xn 1) ~ 1—in rns 
f@=> ee: a aa (=) + Ole! »)| 4} 
yr" X 


n=1 
f4\1—}n 
olen? logt -} | (16) 
ilies” ae 
and this expansion diverges at z = 0 (the Mach cone). Furthermore, the 


expansion found for f does not exist outside the Mach cone. To find the 
value of f in this region it is therefore necessary to consider a solution neat 


the shock-wave. This is done in the next section. 
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nsions Ill. Solution of the equation of motion near the shock-wave 
Substituting equations (4) and (5) in (2) we deduce that 
(1( 4 
f"(B,) —--t- O(B,—«a,). (17) 
(1] (7 1)(1+-a7)* 
constant , 
Similarly r"(B, ——_— +. O(1), 
1> ( Py Xy 
tern nd we infer that a series exists of the form 
My. ; yy 9 -) 
sitior f(y) —+ O(f, v1)" |(\Y—P) 
; . 1 (1 xi)" 
Wiou ) ; 
Y, We + 5x3 + O(B,— 9%) | (Y—B1)? + 
y+1 (1 v7)" { 
| constant : 
. | = +oyly B,)®+-..., (18) 
- Jy Oy : 
that hich may be regrouped as 
) Q 
» CO! . we 2 ‘ Y p 
f(y (B;—o4)°9| 7 + (B,— (2 ) _— (19) 
2 | P34 % 
{ | l ] " 
here g(w) 2 isang smmaneessli lb (20) 
1 (1 xy)" y 1 (1+-aj)* 
Put then 
y , y x, +é(1+w), 
3 Oy (21) 
f(w) &7q(w) &h(w) 
{ ! lequate powers of € in equation (2). 
| The term independent of € vanishes identically, and equating coefficients 
[é we dedue e that 
q’ 9(1+-w)q’” 
5 L)(1 WIIGJ + J om 9 \. (22) 
] XT | XT 
the solution is 
ch as 
16 o\e0..7 ‘ 
ritten g : Tor aa ly 1)(1+aq)*g —(1+w);, (23) 
34 Z 7) 
j ] : : 
ising the fact that g —_ when w = 0 to fix the arbitrary 
; 1 (] ri)" 
1 mstant From (23) we deduce that 
8 4 
— — /(1—3w). (24) 
3() L)(] i)? S(y+1)(1+ a7)" 
d th Since w is negative in the field of flow. q’ is regular therein and g itself is 
g g 
neal 


btained as | g'(w) dw. Other terms in equation (21) could be found in 
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a similar manner, and then equation (21) gives the value of f near the 
shock-wave. We now determine the value of f’ at a point near the shock- 
wave and inside the Mach cone. Put w (k+1), where k is an arbitrary 
positive constant (noticing w 1 on the Mach cone, from equation (21)), 
Then —z = y—a, = k(a,—f,) from equation (21). At this point from 
(24) and (21) we deduce that 


Ty) (B, X,) T O(B, = 04)? 


(25 

IV. Boundary conditions 
The boundary conditions at the shock-wave have already been satisfied. 
The two values for f’ are now equated when z = k(8,—.«,). With this 


value of z, we have from (16), if k be chosen large enough to secure 
convergence (k > 4 is sufficient), 
2 <3 TA, (y+ 1)" Mm (8 1—in ' v 
f(z) = > en Snr) ea oa + Of(B,— ay)! | 4 
— yin 1) | Y | 
7 
+ Oft?"+? log t(B;—a,)'-2"}, (26) 
and putting B,—a, = O(t’) we see immediately from (25) and (26) that 
r= 4, 
k(B 1 -(y+1)M!,, - . 
Put then (P1 au} K% Z #{1+ O(# log t)', 27) 
x \" 


where K is a constant, and substitute (27) in (25) and (26) and equate 
the results; we obtain 


y Lie. @ J 3h 
> 4, K-* = 21848 ff 4 Sh (28 
n=1 k\3° 3n 4 }] 
Now put F(u) = > A, uw"; from equation (15) we find 
n=1 
F(u)—uF'(u) F(u)| F(u) suk’ (u)|. (29 


Solving (29), we have 
, A? u?t,/(Atut+4A? u?) 
F(u) 1 \\44y 1 


) 


(30 


the arbitrary constant and radical sign being fixed by the fact that 
F(u)/u—+> A, asu>0. Now, from (30), 
a ae A‘ ; 4k? 9 
F(K-) cits. ft... 3l 
2k? | A Aj }) 


and, comparing (28) and (31), we see that 
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r the | from equation (12).’ Then, from (27), 


hoc] (y+1)2M8 
| tig we EAP ll 
tra y? 


+ O(t® log t). (33) 


This equation determines the first approximation to the position of the 


iron shock-wave when A is known. 
\ further boundary condition remains to be satisfied, that given by 
equation (3). Substituting the expansion (6) in (3), after replacing fo, f,, 
b and f,, by their power series in y and log y near the body, and equating 
25 powers of ¢ and log sat we deduce the following values for the arbitrary 
ynstants 
v nstan 4 " | | 
h this B J4 M24 MMy+)) (34) 


meee , 20? | 
; C = 2M*—1. 


1 


It can also be shown that by continuing the solution the boundary condi- 
tion at the surface of the body can be satisfied theoretically to any degree 
approximation. Other arbitrary constants in the solution for f are 


é btained by integrating equation (16), to give f, and then introducing 





) that | the result into (21) with 2 = &(B,—a,) and w (k+-1). Since f in (21) 
s O(t8), this shows that the arbitrary constants which may be added to 
‘ , and /, in (8) and (9) must vanish, otherwise there would be terms in 
' nd ## in the value of f obtained by integrating (16). Similarly the 
quate rbitrary constant which may be added to f,, in (11) is zero. Thus the 
statement made after equation (7) is justified. The solution for f near 
the body mav now be written 
" l | pe / at) hs U(y* 
| : :, , M*(y+1) 
nu | 10 Slog dya M? : . - O(y? log Sya) | + 
- 2Zaé : ov 
. | 16 | M 2 f 1 2)) 
| | 7 U(log’ sya t* log Stal (2M? 1)i1 log dya- O(y*)}| t 
O(t® loo Dy (35) 
Y. Derivation of angle of shock-wave, entropy change at shock, 
T 


strength of shock, and pressure coefficient 


| fhe angle of the shock-wave is given from (33) and 
M6 
. i—s (9 1)? #4 O(#® log t). (36) 
3 
sing the sh vave equations, we deduce that the change of entropy 
8 proportional to J/1%¢!*/.° and hence that the strength of shock is 


toportional to .W**/x*?. These results agree with those obtained by 
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Lighthill (4). Using again the shock-wave equations and Bernoulljj 


equation, we may write the pressure coefficient in the following form: 


a ) ’ a 9 » "r\9 "> 1 9 iJ "VOLT 1 , 
= Pe = 2 f—uf')+0%( ful P—f2+ Mf — af Vf +e 
2Po"~ 
(37) 
if we neglect terms in ¢*log*t on the cone. Substituting (35), with y =§ 
in (37) we find 
C, t? — 2#? log $ta+-3a%t! log? dta+ (5.M?—1)t* log dta+ 
M4(y-- 
if, ees . (38) 


x 
This equation gives the first and second approximations to the pressung 
coefficient at the surface of the cone, and it is the same as that found) 
in the author’s paper (1) where the effect of the shock-wave was com 
pletely neglected. Thus it has been demonstrated in this paper that, for 
the case of supersonic flow past a cone, the additional effect of the shock) 
wave does not alter the value of the pressure coefficient to order f. Itig 
reasonable to infer that this result holds in general for pointed bodies of 
revolution which have a convex surface. For a comparison of the analyti¢§ 
and numerical solutions the reader is referred to the earlier paper men=s 


tioned above. 


REFERENCES 
. J. B. Broperick, ‘Supersonic flow round pointed bodies of revolution’, se 
above, pp. 98-120. 
. G. I. Taytor and J. W. Maccott, ‘The air pressure on a cone moving at high 
speeds’, Proc. Roy. Soc. A 139 (1933), 278-311. 
Massachusetts Institute of Technology, Technical Report No. I. Tables @ 
Supersonic Flow around Cones, The Murray Press Company, Cambridge) 
Massachusetts, 1947. 
M. J. Licururn, ‘The position of the shock wave in certain aerodynamié 
problems’, to be published in Quart. J. Mech. and Applied Math. 





oulli’g : -/ 


(38) 
‘essure 
found 

com- 
at, 10r 


shock- 
It is 
dies of 


1a lytie 


r mene 





